August 1, 2009 6:56 WSPC/INSTRUCTION FILE hq-qgp4 



International Journal of Modern Physics E 
© World Scientific Publishing Company 



HEAVY QUARKS IN THE QUARK-GLUON PLASMA 



RALF RAPP 

Cyclotron Institute and Physics Department 
Texas A&M University, College Station, Texas 77843-3366, U.S.A. 
rapp@comp.tamu. edu 

HENDRIK VAN HEES 

Institut fiir Theoretische Physik, Justus- Liebig- Universitdt Giessen, D-35392 Giessen, Germany 
Hendrik.vanHees@theo.physik.uni-giessen.de 

Received July 20, 2009 
Revised July 20, 2009 

Heavy-flavor particles are believed to provide valuable probes of the medium produced 
in ultrarelativistic collisions of heavy nuclei. In this article we review recent progress in 
our understanding of the interactions of charm and bottom quarks in the Quark-Gluon 
Plasma (QGP). For individual heavy quarks, we focus on elastic interactions for which 
the large quark mass enables a Brownian motion treatment. This opens a unique access to 
thermalization mechanisms for heavy quarks at low momentum, and thus to their trans- 
port coefficients in the quark-gluon ffuid. Different approaches to evaluate heavy-quark 
diffusion are discussed and compared, including pcrturbative QCD, effective potential 
models utilizing input from lattice QCD and string-theoretic estimates in conformal 
field theories. Applications to heavy-quark obscrvablcs in heavy-ion collisions are real- 
ized via relativistic Langevin simulations, where we illustrate the important role of a 
realistic medium evolution to quantitatively extract the heavy-quark diffusion constant. 
In the heavy quarkonium sector, we briefiy review the current status in potential-model 
based interpretations of correlation functions computed in lattice QCD, followed by an 
evaluation of quarkonium dissociation reactions in the QGP. The discussion of the phe- 
nomenology in heavy-ion reactions focuses on thermal model frameworks paralleling the 
open hcavy-ffavor sector. We also emphasize connections to the heavy-quark diffusion 
problem in both potential models and quarkonium regeneration processes. 



1. Introduction 

The investigation of strongly interacting matter constitutes a major challenge in 
modern nuclear and particle physics. Of particular interest are phase changes be- 
tween hadronic and quark-gluon matter, similar to the one which is believed to 
have occurred in the early Universe at a few microseconds after its birth. While the 
theory of the stron g force is by now well established in terms of Quantum Chro- 
modynamics (QCD)Il213|^ ^.^^ ^j- ^-^^jQi- manifestations in the world around us 
- the confinement of quarks and gluons and the generation of hadronic masses - 
are subject of vigorous contemporary research. Both phenomena occur at cnergy- 
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momentum scales of Q'^ < 1 GeV^ where the QCD coupUng constant is rather large, 
tts > 0.3, and therefore perturbation theory is not reliable and/or applicable. In a 
hot and dense medium at sufficiently large temperature (T) and/or quark chemical 
potential (/ig), one expects the finite-size hadrons to be dissolved into a dcconfined 
Quark-Gluon Plasma (QGP) where the condensates underlying hadronic-mass gen- 
eration have melted. Numerical simulations of latticc-discrctizcd QCD (IQCD) at 
finite temperature predict the phase change from hadronic to quark-gluon matter 
to occur at a "pseudo-critical" temperature of Tc ~ 200 Me\Sl. This appears to be 
a rather small scale for a "perturbative QGP" (pQGP) of weakly interacting quarks 
and gluons to be realized, even though the computed energy density matches that 
of an ideal (non-interacting and massless) QGP within 20% or so for T > 1.2 Tc- 

In the laboratory, one hopes to create a QGP by colliding heavy atomic nuclei 
at ultrarelativistic energies, with a centcr-of-mass energy per colliding nuclcon pair 
well above the nucleon rest mass, \/s/A 3> Mn. If the energy deposition in the reac- 
tion zone is large enough, and if the interactions of the produced particles are strong 
enough, the notion of an interacting medium may be justified, despite its transient 
nature. This notion has been convincingly verified in nuclear collision experiments 
over the last ^25 years at the Super-Proton-Synchrotron (SPS) at the European 
Organization for Nuclear Research (CERNj^ and at the Relativistic Heavy-Ion Col- 
lider (RHIC) at Brookhaven National Laboratory (BNLjSl. Transverse-momentum 
{Pt) spectra of different hadron species in the low-pT regime {px ^ 2-3 GeV) reveal 
that the produced medium explodes collectively reaching expansion velocities in 
excess of half the speed of light. In the high-pT regime {pr ^ 5 GeV), which in the 
heavy-ion environment became available at RHIC for the first time, hadron spectra 
are suppressed by up to a factor of ~5 relative to p-p collisionSjindicative for a 
strong absorption of high-energy partons traversing the mediunl^J. The inclusive 
production of charm-quark bound states {J/ip mesons) is suppressed by a factor 
of 3-5 at both SPS and RHIC, indicat ive for their dissolution in the medium (pos- 
sibly related to deconfincmcnt jHl 9| 1 | large excess of electromagnetic radiation 
(photons and dileptons) is observed, indicative for medium temperatures around 
200 MeV and a "m elting" of the p-meson resonance (possibly related to hadronic 
mass de-generation)PJJ12l A more differential analysis of hadron spectra in non- 
central Au-Au collisions at RHIC reveals a large elliptic asymmetry of the collective 
flow ("elliptic flow"): the spatial asymmetry of the initial nuclear overlap zone is 
converted into an opposite asymmetry in the final hadron px spectra. Within a 
hydrodynamic modeling of the exploding fireball this observation requires a rapid 
thermalization and a very small viscosity of the interacting mcdium-'^^ ^"^ 15il6j^ 
Only then can spatial pressure gradients build up fast enough to facilitate an ef- 
fective conversion into azimuthal asymmetries in the energy-momentum tensor of 
the system. The agreement of hydrodynamic predictions with elliptic-flow data at 
RHIC led to the notion of an "almost perfect liquid" , with a ratio of viscosity to 
entropy density close to a conjectured lower bound of any quantum mechanical 
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systenJl^. The microscopic mechanisms underlying these rather remarkable trans- 
port properties are yet to be determined. In this context, heavy quarks (charm and 
bottom, Q=c and b) and their bound states (charmonia and bottomonia) are rec- 
ognized as particularly suitab le pr obes of the medium produced in ultrarelativistic 
heavy-ion collisions (URHICsp^- In the present article we will attempt to review 
the current status of the theory and phenomenology of this promise. 

Let us first focus on the sector of individual heavy quarks (open heavy flavor). 
The fact that their masses are well above the typical temperature of the system, 
niQ 3> T, has at least three important implications: 

(1) The (hard) Q Q p roduction process is essentially restricted to primordial 
N-N collisions^, i.e., re- interactions in the subsequently evolving medium 
are not expected to change the number of hea vy quarks (reminiscent of the 
"factorization theorem" of perturbative QCIJ^U); this is borne out exper- 
imentally by a scaling of cc production, Ncs, with the number of binary 
N-N collisions, A^coii, at different collision centralitie*^. 

(2) The thermal relaxation time of heavy quarks ought to be larger than for 
light quarks, parameterically by a factor ^mq/T « 5-20. With a light- 
quark and gluon thermalization time of Tq^g ~0.3-lfm/c (as indirectly 
inferred from hydrodynamic modeling at RHIC) and an estimated QGP 
lifetime of tqqp ~ 5fm/c in central Au-Au collisions, one expects {jh) 
to be on the same order as (significantly larger than) tqgp. Thus, charm 
(and especially bottom) quarks are not expected to reach thermal equilib- 
rium, but their re-interactions should impart noticeable modifications on 
the initial momentum spectrum (less pronounced for bottom). The final 
heavy-quark (HQ) spectra may therefore encode a "memory" of the inter- 
action history throughout the evolving fireball, by operating in between the 
limits of thermalization and free streaming. 

(3) The theoretical task of describing HQ interactions is amenable to a diffusion 
treatment, i.e., Brownian motion of a heavy test particle in a bath of a light- 
particle fluid. Nonrelativistically, the typical thermal momentum of a heavy 
quark is p^^ ~ ^rnqT ^ T^, and therefore much larger than the typical 
momentum transfer from the medium, T^. This allows to expand the 
Boltzmann equation in momentum transfer to arrive at a Fokkcr-Planck 
description of HQ diffusion in the QGP, which directly yields the pertinent 
transport coefficients as well. 

The above three points provide a well-defined framework to construct in-medium 
HQ interactions in QCD matter and test them against observables in URHICs 

''The (weak-decay) lifetime of the top quark of ~0.1 fm/c is too short to render it a viable probe 
in URHICs; thus, heavy quarks will exclusively refer to charm and bottom in this article. Strange 
quarks are in between the heavy- and light-quark limit, forming their own complex of valuable 
observableJEl. 
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(quantitative comparisons additionally require to account for effects of hadroniza- 
tion of the quarks, as well as reinteractions in the hadronic medium). The Fokker- 
Planck approach is readily implemented for the c ase of elastic p + Q ^ p + Q 
scattering off partons in the medium {p = q, q, g j23 | 24 | 25 | 26 | 27 | 28 | ^^le light- 
hadron sector, however, the large suppression of high-py spectra is believed to be 
largely caused by radiative energy loss of high-energy partons traversi ng the QGP, 
i.e., medium-induced gluon radiation of type q Even in the 

low-pT regime, perturbative 2,3^3 scattering processes have been suggested to 
facilitate the rapid thermalization required by phenomenology (albeit in connection 
with rather large coupling constants of as ~ O.sj^. The situation could be quite 
different in the HQ sector. In the low-momentum limit, giuon-Bremsstrahlung of 
a heavy quark is suppressed^ and the dominant momentum-transfer reaction is 
elastic scatterin^^SI^ jg ^gjj known from classical electrodynamics, the radia- 
tive energy loss of a muon is suppressed relative to an electron by a mass ratio 
{m^/me)'^. In perturbative QCD (pQCD) it is currently an open question at what 
momentum scale radiative energy loss of a heavy quark takes over from the coUi- 
sional one (which, most likely, will depend on additional parameters such as tem- 
perature, path length, etc.). In fact, this may not even be a well-defined question 
since nonperturbativc processes at moderate momentum transfers may supersede 
perturbative ones before the elastic part of the latter dominates over the radiative 
one. The relations between perturbative and nonperturbativc interactions is one of 
the key issues to be addressed in this review. 

Experimental signatures for the modifications of HQ spectra in URHICs are 
currently encoded in single-electron (e*) spectra associated with the semileptonic 
decays of charm and bottom hadrons, Ac, . . . — > evX. These measurements 

require a careful subtraction of all possible "photonic" sources of electrons, such as 
photon conversions in the detector material, Dalitz decays of tt and 77, vector-meson 
decays, and others. The modifications of the "non-photonic" electron spectra (asso- 
ciated with heavy-flavor decays) in Au-Au collisions are then quantified by the stan- 
dard nuclear modification factor, R^^, and elliptic flow coefficient, w|. The available 
RHIC data in semicentral and central Au-Au collisions at -y/ijvw = 200 GeV ex- 
hibit a substantial elliptic flow of up to — 10% and a large high-pT suppression 
down to — 0.25, respectively^*^ ^•^ "^"^ 25j Both values are quite comparable 

to those for light hadrons (the pion V2 reaches somewhat higher, to about 15%). 
Radiative energy-loss modelJ^Sl based on perturbative QCD cannot explain the 
data. These data were, in fact, instrumental^ in reco nsidering ela stic scattering 
as a significant source of parton energy loss in the QG P ' ' ' The combina- 
tion of pQCD elastic and radiative scattering does not suffice either to reproduce 
the observed suppression once a realistic bottom component is accounted for in 
the electron spec tra^. Elastic scattering based on nonperturbativc interactions, as 
proposed in Refs.^^SI^ simultaneously accounts for the elliptic flow and suppres- 



sion reasonably welP^. This has reinforced the hope that HQ observables provide 
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the promised precision tool to characterize transport properties of the "strongly 
coupled QGP" (sQGP). E.g., if a clear mass hierarchy in thermal relaxation times, 
TQ oc ^^Tq (as well as Tb/rc = mi,/mc) emerges from a quantitative analysis of 
URHIC data, it would be suggestive for a universal behavior of light- and heavy- 
quark transport in the QGP. However, there is still a substantial way to go before 
such a program can be realized, as discussed belo'wl^. 

It should not be surprising if in-medium properties of open heavy flavor, espe- 
cially at low momentum, are closely related to medium modifications of heavy 
quarkonia. The latter have a long history as "probes" of the QGP in heavy- 
ion coll isions, especially as potential indicators of the deconfinement transition, 
cf. Refs.lS^^Iini fQj- a broad up-to-d ate cover age of this topic. In particular, progress 
in finite-temperature lattice has triggered vigorous reconsideration 

of the question whether quarkonia, especially their ground states, can survive in 
the QGP significantly above the critical temperature. These developments include 
the application of potential models at finite temperature, coupled with the hope 
that heavy-quark free energies as computed in thermal IQCD can serve as a model- 
independent input for the low-energy heavy-quark interaction. If charmonium bind- 
ing indeed remains sufficiently strong in the QGP to support bound states up to 
rather high temperatures, it is conceivable that the underlying int erac tion is of a 
more general rele vance and therefore also operative in hcavy-lighlf^ and maybe 
even light-lighil^^l systems. Especially in the former case, from the point of view of 
elastic (on-shell) scattering of a heavy quark in the medium, the conditions for mo- 
mentum transfer are comparable to the heavy-heavy interaction governing quarko- 
nium properties. Since low- momentum HQ interactions determine their transport 
properties, one immediately recognizes an intimate relation between HQ transport 
and in-medium quarkonia. These connections are also being exploite d in the analysis 
of thermal IQCD computations of quarkonium correlation functiontPSl. In addition 
to the binding properties, the inelastic reaction rates of quarkonia with surround- 
ing partons or hadrons are a key ingredient for a quantitative description of their 
spectral function in QCD matter (also here "quasi-elastic" scattering of thermal 
partons with a heavy quark inside the quarkonium bound st ates may play an im- 
portant role, especially if the binding energy becomes smalP^). A good control 
over all of these aspects is mandatory to utilize quarkonium properties as diagnos- 
tic tool in heavy-ion collisions and eventually deduce more general properties of the 
medium produced in these reactions. As in the open heavy-flavor sector, this has 
to be built on a solid knowledge of the space-time history of nuclear collisions, as 
well as of the initial conditions on quarkonium spectra. The latter aspect could be 
more involved than for single heavy quarks, since (a) measurements inp-A collisions 
show that cold- nuclear-matter (CNM) effects from the incoming nuclei (e.g., the so- 
called nuclear absorption) affect the primordial charmonium number significantly 

''A rece nt r eview articlJSSl addresses similar topics but from a more elementary perspective; see 
also Ref.Hni 



August 1, 2009 6:56 WSPC/INSTRUCTION FILE hq-qgp4 



6 R. Rapp, H. van Hees 

(e.g., with up to 60% suppression for J/ij] at SPS energies when extrapolated to 
central Pb-Pb collisions); (b) the bound-state formation time introduces another 
rather long time scale (soft energy scale) which is easily of the order of (or longer 
than) the thermalization time of the medium (at least for charmonia and excited 
bottomonia at RHIC energies and higher). 

Charmonium suppression beyond the level of CNM effects has been established 
in semi-/central Pb-Pb and Au-Au collisions at SP^31! and RHICpSl, respectively. 
An intriguing finding is that the observed suppression pattern and magnitude is 
very comparable at SPS and RHIC, despite the different collision energies which 
lead to substantial variations in, e.g., light-hadron observables (most notably a fac- 
tor ~2 larger charged particle rapidity density and s tron ger collective phenomena 
at RHIC). However, this "degen eracy" wa s predicted^ as a consequence of char- 
monium regeneration mechanismiSHISlISSl ^ stronger suppression in the hotter and 
denser medium at RHIC is compensated by the coalescence of c and c quarks in 
the QGP and/or at hadronization (the cc production cross section at RHIC is 
about a factor of '--'100 larger than at SPS energies). While an "extra" source of 
charmonia increases the complexity of pertinent observables in heavy-ion reactions, 
it also provides another, rather direct, connection between the open and hidden 
heavy-flavor sectors. Obviously, the secondary yield from c-c coalescence necessar- 
ily carries imprints of the charm-quark distributions, both in its magnitude (softer 
c-quark spectra are expected to result in larger coalescence probabilities) and in 
its momentum spectra (including elliptic flow). A comprehensive theoretical and 
phenomenological analysis of open and hidden heavy flavor is thus becoming an 
increasingly pressing and challenging issue. As a final remark on quarkonia in this 
introduction, we point out that bottomonium production in heavy-ion reactions 
is less likely to receive regeneration contributions (at least at RHIC and possibly 
neither at LHC). In addition, the increase in bottomonium binding energies (com- 
pared to charmonia) render them rather sensitive probes of color screening which 
strongly influences its dissociation rateJ^SI. Bottomonia thus remain a promising ob- 
servable to realize the originally envisaged "spectral analysis of strongly interacting 
matter"151. 

Our review is organized as follows: In Sec. [2] we outline the theoretical frame- 
work of evaluating HQ diffusion in equilibrium QCD matter. We first recall basic 
steps in setting up the HQ diffusion equation (Sec. 12. 1[) which determines the time 
evolution of the HQ distribution function in terms of pertinent transport coeffi- 
cients based on elastic scattering amplitudes. This is followed by a discussion of 
several microscopic approaches to calculate the HQ friction and diffusion coeffi- 
cients in the QGP: perturbative QCD (Sec. [Ml at leading ([22111 [M21 iMiSl) and 
next-to-leading order p.2.4p as well as for three-body scattering (|2.2.5|) ; nonpertur- 
bative calculations (Sec. 12. 3p implementing resonance-like correlations in the QGP 
using HQ effective theory (Sec. 12.3. ip . in- medium T-matrices with HQ potentials 
estimated from thermal lattice QCD (Sec. I2.3T2|) . or coUisional-dissociation mech- 
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anisms of heavy mesons (Sec. 12.3. 3^ : and string-theoretic evaluations based on the 
conjectured correspondence to conformal field theories fSec. l2.4[ ). The variety of the 
proposed approaches calls for an attempt to reconcile the underlying assumptions 
and basic interactions (Sec. [23]). This is followed by a discussion of inelastic (radia- 
tive) energy-loss calculations and their relation to elastic ones fSec. l2.6( ). We briefly 
consider interactions of open heavy-flavor hadrons in hadronic matter (Sec. 12.7]) . In 
Sec. [3] we discuss applications of HQ diffusion to URHICs using relativistic Langevin 
simulations of the Fokker-Planck equation within an expanding finite-size thermal 
medium (Sec. 13. ip . A realistic description of the latter (utilizing hydrodynamics, 
transport models or suitable parameterizations thereof) is an essential prerequisite 
to enable a quantitative extraction of transport properties of the QCD medium 
(Sec. 13. 2p . Further ingredients are reliable initial conditions (possibly modified by 
nuclear effects) and the conversion of quarks to hadrons (Sec. 13.3]) . Implementations 
of different HQ diffusion coefficients in various space-time models are quantitatively 
analyzed in terms of the resulting HQ spectra at RHIC, in particular their nuclear 
modification factor and elliptic fiow (Sec. l3.4( ). Including effects of hadronization (as 
well as semileptonic electron decays), a quantitative comparison of these calcula- 
tions to single-electron spectra at RHIC is conducted (Sec. 13. 5|) . We emphasize the 
importance of a consistent (simultaneous) description of pt spectra and elliptic flow. 
Only then can these observables be converted into a meaningful (albeit preliminary) 
estimate of charm- and bottom-quark diffusion coefficients in the QGP. We finish 
the discussion on open heavy flavor with an attempt to utilize these coefficients for 
a schematic estimate of the ratio of shear viscosity to entropy density in the QGP 
(Sec. 13.6)) . In Sec. [3] we elaborate on theoretical and phenomenological analyses of 
quarkonia in medium and their production in heavy-ion collisions. We first address 
spectral properties of quarkonia in equilibrium matter (Sec. 14.1]) : Euclidean corre- 
lation functions computed in lattice QCD with good precision have been analyzed 
in terms of potential models based on screened HQ potentials (Sec. I4.1.ip . The 
interplay of color screening and parton-induced dissociation reactions has impor- 
tant consequences for the evaluation of quarkonium dissociation widths fSec. l4.lT2)) . 
In light of the charmonium equilibrium properties the current status of the phe- 
nomenology in heavy-ion collisions is discussed fSec. 14. 2| . First, quarkonium trans- 
port equations are introduced along with their main ingredients, i.e., dissociation 
widths and equilibrium numbers using relative chemical equilibrium at fixed HQ 
number fSec. I4.2T|) : this is followed by model comparisons to J/tp data at SPS and 
RHIC, scrutinizing suppression vs. regeneration mechanisms and their transverse- 
momentum dependencies (Sec. I4.2.2p . and a brief illustration of predictions for T 
production at RHIC. In Sec. [5] we recollect the main points of this article and 
conclude. 
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2. Heavy-Quark Interactions in QCD Matter 

At an energy scale of the (pseudo-) critical QCD transition temperature, the large 
charm- and bottom-quark masses imply that the HQ diffusion problem is a non- 
relativistic one (unless initial conditions bring in an additional large scale). In the 
weak-coupling regime this further implies that the dominant interactions of the 
heavy quark are elastic scattering (gluon radiation is suppressed by an extra power 
in Us and cannot be compensated by a large momentum transfer as could be the 
case for a fast quark; see, e.g., the discussion in Ref .^^^ . It turns out, however, that 
the perturbative expansion of the charm-quark diffusion coefficient, evaluated using 
thermal field theory is not well convergent even for a strong coupling constant as 
low as ttc 

= O.lESl Thus, non-perturbative methods, e.g., resummations of large 
contributions or interactions beyond perturbation theory, are necessary to improve 
the estimates of HQ diffusion. This is not surprising since transport coefficients 
usually involve the zero-momentum limit of correlation functions rendering them 
susceptible to threshold effects which may increase with the mass of the particles. A 
simple example of such kind are Coulomb-like bound states (e.g., heavy quarkonia), 
where the binding energy increases with increasing HQ mass, es cx a^mg, to be 
compared to thermal effects, e.g., at a scale ^gT for Debye screening (to leading 
order in g) or at for inelastic dissociation reactions with thermal partons. An 
interesting question in this context is whether potential models are a viable means 
to evaluate HQ interactions in the QGP. If a suitable formulation of a potential at 
finite temperature can be established, a promising opportunity arises by extracting 
these from first principle lattice computations of the HQ free energy. In the heayy- 
quarkonium sector such a program has been initiated a f ew years agc l^S1 57 | 58 | 59 | 6 0l 
with fair success, although several open questions remairl^IElIMES] applicable, 
potential models have the great benefit of allowing for nonperturbative solutions 
utilizing Schrodinger or Lippmann-Schwinger equations; the calculated scattering 
amplitudes can then be straightforwardly related to transport coefficients. A key 
issue in this discussion is the transition to the (ultra-) relativistic regime, which 
becomes inevitable in applications to experiment toward high momentum. While 
relativistic kinematics can be readily accounted for, the opening of inelastic (ra- 
diative) channels poses major problems. However, here the contact to perturbative 
calculations may be possible and provide a valuable interface to match the differ- 
ent regimes, at least paramctrically (e.g., in the limit of a small coupling constant 
and/or high temperature). This reiterates the importance of identifying the com- 
mon grounds of seemingly different calculations for HQ properties in medium. 

We start the discussion in this Section by setting up the Brownian Motion 
framework for heavy quarks in the QGP (Sec. 12. ip . The main part of this Section 
is devoted to the evaluation of the Fokker-Planck transport coefficients. We focus 
on elastic interactions, classified into (various levels of) perturbative (Sec. 12. 2| ) and 
nonperturbative approaches (Secs. [2?3l and l2.4]) . As wc will sec, there is considerable 
conceptual overlap in the calculations available in the literature, the main difference 
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being that they are carried out in different approximation schemes (Sec. [^75]) . Our 
presentation also encompasses inelastic reactions with an additional gluon in the 
final and/or initial state, i.e., radiative energy-loss calculations within perturbative 
QCD fSec. l2.6( ). This raises the issue of their relative magnitude compared to elastic 
interactions which has recently received considerable re-consideration even for light 
quarks and gluons. Finally, we address interactions of hadrons carrying charm or 
bottom in hadronic matter (Sec. 12. 7p . Bottom- up extrapolations in temperature 
(or density) in the hadronic world are useful complements to top-down ones in the 
QGP, to reveal qualitative trends of, e.g., the HQ diffusion coefficient toward Tc- 



2.1. Heavy- Quark Diffusion in the Quark-Gluon Plasma 

As emphasized in the Introduction, an attractive feature in analyzing H Q m otion 
in a QGP is the ensuing simplification to a Brownian motion frameworl*^. The 
latter is characterized by a Fokker-Planck equation where HQ interactions are con- 
veniently encoded in transport coefficients. These, in turn, are readily related to 
underlying (elastic) scattering matrix elements on light partons in the QGP which 
allow for direct comparisons of microscopic models of the HQ interaction (as elab- 
orated in subsequent sections). 

Starting point for the derivation of the Fokker-Planck equatiorP^ is the Boltz- 
mann equation for the HQ phase-space distribution, /g. 



d p d „ d 

\- — 1- F — 

dt u!p dx dp 



fQ{t,X,p)=C[fQ], (1) 

where iVp = A- p^ denotes the energy of a heavy quark with three-momentum 

p, F is the mean-field force, and C[/q] summarizes the collision integral which 
will be analyzed in more detail below. In the following, mean-field effects will be 
neglected, and by integration over the fireball volume, Eq. ([1]) simplifies to an 
equation for the momentum distribution, 

^/Q(i,p)=C[/Q], (2) 

where 

/Q(i,p) = j d^xfQ{t,x,p). (3) 

The collision integral on the right-hand side of Eq. ^ encodes the transition rate 
of heavy quarks due to collisions into and out of a small momentum cell d^p around 
the HQ momentum p, 

= J d^k[w{p + k,k)fQ{p + k)-w{p,k)fQ{p)] . (4) 

Here w{p, k) is the transition rate for collisions of a heavy quark with heat-bath 
particles with momentum transfer fe, changing the HQ momentum from p to p~k. 
Accordingly the first (gain) term in the integral describes the transition rate for 
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HQ scattering from a state with momentum p + k, into a state with momentum p, 
while the second (loss) term the scattering out of the momentum state p. 

The transition rate, w, can be expressed through the cross section of the collision 
processes in the heat bath. For elastic scattering of a heavy quark with momentum 
p on a light quark in the heat bath with momentum q, one finds 

j^fq,giq)vrd-^{p,q-> p- k,q + k), (5) 

where fq^g are the Fermi or Bose distributions for thermal light quarks or gluons, 
and 7g = 6 or 7g = 16 the respective spin-color degeneracy factors. The relative 
velocity is defined as 



V^e\ = , (6) 

where p = {ujp, p) and q = {ujq, q) are the four momenta of the incoming heavy and 
light quark, respectively. Upon expressing the invariant differential cross section, 
dcr/dri, in Eq. ^ in terms of the spin-color summed matrix element, l-^P' the 
collision term, Eq. ([4]), takes the form 

r\f [ [ ^P_ f dV 1 ^ 

^■"^^ 2upj {2TTf2ojJ {2nf2ujpJ (27r)32u;q, 7Q ^ ' (7) 

X (2^)4jW(p + q-p'- q')[fQ{p')UM) - fQ{p)UAQ)] 

with k = p — p' — q' — q. 

The key approximation is now that the relevant momentum transfers to the 
heavy quark obey |A;| <^ \p\. This enables to expand the HQ momentum distribu- 
tion function, /q, and the first argument of the transition rate, w, in the collision 
integral, Eq. (|4|), with respect to fe up to second order0 

wip + k, k)fQ{p + k, k) ~w{p, k)fQ{p) 

d I d"^ 

+ k—[w{p,k)fQ{p)] + -fc.fcj ^^^^ [w(p, k)fQ{p)] 

(8) 

(i,j=l,2,3 denote the spatial components of the 3-vectors, with standard summa- 
tion convention for repeated indices). The collision integral then simplifies to 

d 1 , , 92 ■ 

dpi 2 PiPj 

i.e., the Boltzmann equation ^ is approximated by the Fokker-Planck equation, 
pQit^P) = ^ [MP)fQ{t,p) + ^[i?.,(p)/Q(t,p)]| . (10) 

"^According to the Pawula theorerrlSH any truncation of the colhsion integral at finite order is 
only consistent with fundamental properties of Markov processes if the truncation is made at the 
2"'^-order term. 



cVq] ~ / d^fc 



[w{p,k)fQ{p)] , (9) 
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The drag and diffusion coefficients are given according to Eq. © by 

A,{jp) = i d^kw{p,k)ki , 

1 f (11) 

For an isotropic background medium, especially in the case of (local) equilibrium 
(implying that the coefficients are defined in the local rest frame of the heat bath), 
rotational symmetry enables to simplify the cocfhcicnts to 

B,, ip) = Bo {j>)4 ip) + Bi (p)i^t (p) , 

where the projection operators on the longitudinal and transverse momentum com- 
ponents read 

= Pm = S.-^- (13) 

Implementing these simplifications into the collision integral, Eq. ([7|), the scalar 
drag and diffusion coefficients in Eq. are given by integrals of the form 

^ ^^>^ 2cjpJ {2^f2uJ (2^)32c.pJ (2^)32c.q,7Q^'-^' (^4) 
In this notation, the coefficients can be written as 

- (' - ^ 
1 / _ (rt>)^ 

Note that Eq. includes the sum over gluons and light quarks (w, d, s). 

The physical meaning of the coefficients becomes clear in the non-rclativistic ap- 
proximation of constant coefficients, 7 = A{p) = const and D = Bq{p) = -Bi(p) = 
const, in which case the Fokkcr-Planck equation further simplifies to 

^/Q(i,P) = 7^[p./Q(t,P)] + DAMt,p) . (16) 

E.g., for an initial condition 

/Q(t = 0,p) = <5(^)(p-Po) , (17) 
the solution takes the form of a Gaussian distribution, 



— 3 /2 

fQit,p) - [1 - exp(-27i)]} cxp 



7 [P - Po exp(-7t)] 
2D l-exp(-27t) 



(18) 
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From the equation for the mean momentum, 

(p) = poexp(-7t) , (19) 

one sees that 7 determines the relaxation rate of the average momentum to its 
equiUbrium value, i.e., it is a drag or friction coefficient. The standard deviation of 
the momentum evolves according to 

on 

(P')-(P)' = — [l-cxp(-270] , (20) 

7 

i.e., D is the momentum-diffusion constant, describing the momentum fluctuations. 

In the limit t ^ 00, Eq. psp approaches the (non-relativistic) Boltzmann dis- 
tribution, 

Since in thermal equilibrium the heavy quarks have to obey an equilibrium distri- 
bution with the temperature, T, of the heat bath, the drag and diffusion coefficients 
should satisfy the Einstein dissipation-fluctuation relation, 

D = mQjT . (22) 



The relativistic Fokker-Planck equation will be discussed in Sec. 13. II in connection 
with its formulation in terms of stochastic Langevin equations. 

We note that the spatial diffusion coefficient, £><,, which describes the broadening 
of the spatial distribution with time, 

(x^t)) - {x{t)f ~ 6Dst , (23) 

is related to the drag and momentum-diffusion coefficient through 

T 

Ds^ ^jr. (24) 

mQ7 D 



2.2. Perturbative QCD Approaches 

In a first step to evaluate HQ diffusion in a QGP perturbation theory has been ap- 
plied, thereby approximating the medium as a weakly interacting system of quark 
and gluon quasiparticles. Such a treatment is expected to be reliable if the temper- 
ature is large enough for the typical momentum transfers, ~ T^, to be in the 
perturbative regime, > 2 GeV^ or so. This is most likely not satisfied for matter 
conditions realized at SPS and RHIC. For more realistic applications to experiment 
several amendments of the perturbative approach have been suggested which are 
discussed subsequently (focusing again on elastic HQ scattering on light partons). 
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Fig. 1. Feynman diagrams for leading-order perturbative HQ scattering off liglit partons. 
2.2.1. Schematic Leading Order 

The init ial estimates of equilibration times and energy loss of heavy quarks in the 
QQpl23] jjg^vg started from the leading-order (LO) perturbative diagrams involving 
the minimum o f tw o strong-interaction vertices, as displayed in Fig. [TJ Pertinent 
matrix elementJ^ figuring into Eq. (jl4p in the vacuum have been computed in 
Ref.E51. 

The dominant contribution arises from giuon t-channel exchange, i.e., the 
3'''^ and 4**^ diagram in Fig. [T] For forward scattering, the gluon propagator develops 
the well-known infrared singularity which has been regularized by introducing a 
Debye-screening mass, 

G{t) = i ^ , fiD^gT , (25) 

t t — 

where g = v^SttoI denotes the strong coupling constant. Even for a value as large 
as as = 0.4, and at a temperature of T = 300 MeV (typical for the early stages in 
heavy-ion collisions at RHIC), the thermal relaxation time, Toq = 1/7, for charm 
(bottom) quarks turns out around ^15(40) fm/c (and therefore much larger tha n a 
typical QGP lifetime of ~5fm/c at RHIC), see, e.g., right panel of Fig.[7](in Rem 
the corrections due to quantum-equilibrium distributions (Bose/Fermi) have been 
investigated and found to be small). Note that with the above gluon propagator, the 
pertinent total HQ-parton cross section is parametrically given by aqp cx a'^/^'jj, 
i.e., it essentially increases only linearly in (p = 5, q,g). 



2.2.2. Leading Order with Hard Thermal Loop Resummation 

In Ref.l^, the schematic introduction of the Dcbye mass into the t-channcl gluon- 
exchange propagator has been extended by a LO hard-thermal loop (HTL) calcula- 
tion of the charm-quark drag and diffusion coefficients in the QGP. In this approach, 
the screening of the gluon propagator in the t-channel diagrams (Fig. [T|) is realized 
by inserting the HTL gluon propagator for the region of small momentum exchange. 
In Coulomb gauge, with q = |q|, this propagator is given by 

G^^(tj, q) = 2 , „ + — 2 I TT ' (26) 

q^ + Hoo q^ - ijj^ + Ht 
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Fig. 2. (Color online) HQ transport coefficients in HTL improved perturbation theorjESl. Left 
panel; spatial diffusion coefficient at p = as a function of an independently varied Debye mass, 
mo = fJ'D} figuring into the t-channel gluon exchange propagator, for different quark- flavor content 
of the (Q)GP. Right panel: momentum dependence of the drag coefficient, rjoip) = j4(p), for three 
values of fi£) in t-channcl gluon exchange; the lower curve, with dp/dt oc v, resembles a calculation 
in the non-relativistic limit (A/ = rric = 1.4 GeV in the x-axis label denotes the charm-quark mass). 



where the i,j G {1,2,3} denote the spatial components oi fi^iy £ {0,1,2,3}. The 
HTL self-energies read 



In 



(- 




— m 


1 


\<1 


-uj) 







In ( — I — ITT 

q — LU 



(27) 



For small energy transfers, ui, and a slowly moving heavy quark, v 1, only the 
time component of the propagator contributes to the squared matrix elements which 
in this limit reduces to the Debye-screened Coulomb- like propagator, Eq. ([25]) . Fig. [2] 
shows the spatial diffusion coefficient and the momentum dependence of the drag 
coefficient resulting from this calculation. Compared to the screening description 
with a constant Debye mass, the drag coefficient shows a slight increase for an 
intermediate range of momenta (cf., e.g., the pQCD curves in Fig. I12p . 



2.2.3. Leading Order with Running Coupling 

As indicated in the Introduction, the current data situation at RHIC does not allow 

for an understanding of the electron data in terms of LO pQCD with reasonably 

small couplin g constant (say, as < 0.4). This was a motivation for more recent 
studiePEII^ 

augmenting the LO pQCD framework in search for stronger effects. 
Two basic amendments have been introduced. First, the idea of of introduc- 

ing a reduced screening mass in the gluon propagator was made more quantitative. 
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Gr{t) oc - — . (28) 



Starting from an ansatz for the screened gluon propagator, 

1 

t - r/^^ 

the objective is to obtain an estimate for the constant r (it was denoted k, in 

]^g£g|6o|671 changed the notation to avoid conflicts in what follows below). This 

has been done in analogy to a corresponding QED calculatiorPSESl^ requiring 

that the energy loss of a high-energy quark obtained in a LO-pQCD calculation with 

the screened propagator, Eq. ([28|) . matches a calculation where for low momentum 

transfers, \t\ < \t*\, the HTL propagator, Eq. (pS)) . and for \t\ > \t* \ the perturbativc 

gluon propagator, Eq. ([25]) . i s used: \t* \ is a momentum-transfer scale between g^T^ 

and T^. The QED calculatio n"" ! "^ ! yields an energy loss which is independent of the 

matching scale while this is not the case in QCD. This problem is treated by 

introducing an infrared-regulator mass into the hard part of the energy-loss integrals 

involving the ^-channel exchange-matrix elements, chosen such that the dependence 

on is weak for \t*\ < (the validity range of the HTL approximation). This 

translates into effective values for the r coefficient in Eq. ([28]) of r ~ 0.15-0.2. 

Second, a running strong coupling constant is introduced well into the non- 

perturbativ e re gime but with an infrared-finite limit. The justification for such 

a procedurj^ is that it can account for (low-energy) physical observables (e.g., 

in e+e~annihilatiorPIJ) in an effective way. The parameterization adopted in 
Refs.lMi67|jg i^j^ggj 

on an extrapolation of Ref.^^ into the spacelike regime, 

, Att iL-^ for Q2 < 

/3o \l/2-7r-iarctan(L+/7r) for > q, 

where f3o = 11-27V//3, Nf = 3, and L± = \n{±Q^/A^). The pertinent substitution 
in the f-channel gluon-exchange matrix elements amounts to 

7-?^, (30) 

where the regulator mass is chosen as jl^ € [l/2,2]fljj, while the Debye-screening 
mass is determined self-consistently from the equation 

r.2 _(Nc^Nf 



nD-[^f+ 6 j ^^<-~^d)T' ■ (31) 

To find the optimal value for the regulator mass a similar strategy of matching 
the energy loss with a Born approximation has been employed, using the substitu- 
tion, Eq. (j30p . in the i-channel diagrams, with a HTL calculation along the same 
lines as summarized above for the calculation with non-running a^- The results 
for the drag coefficients for charm quarks under the various model assumptions de- 
scribed above are depicted in Fig. [31 Changing the screening mass from the standard 
Debye mass, fin, to that reproducing the HTL energy loss, with r = 0.15 in Eq. ([28]) . 
increases the drag coefficient by a factor of 2. In view of the large reduction in r 
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Fig. 3. (Color online) The drag coefficient as a function of HQ three-momentum in the ar nended 
pQCD scheme with reduced infrared regulator and running coupling constant (left 

paneljSMZl. 

The corresponding legend (right panel) details the different parameter choices in the calculation. 



this appears to be a rather moderate effect. This is simply due to the fact that the 
change mostly enhances forward scattering which is little effective in thermalizing 
(isotropizing) a given momentum distribution. Implementing the running-coupling 
scheme with a small screening mass yields a substantial enhancement by a factor 
of -5. 

2.2.4. Next-to- Leading Order 

The rather large values of the coupling constant employed in the calculations dis- 
cussed in the previous sections imminently raise questions on the convergence of the 
perturbative series. This problem has been addressed in a rigorous next-to-leading- 
orde r (NLO ) calculation for the HQ momentum-diffusion coefhcicnt, k = 2D, in 
Refs.^SSEH This work starts from the definition of k as the mean squared momen- 
tum transfer per unit time, which in gauge theories is given by the time-integrated 
correlator of color-electric-field operators connected by fundamental Wilson lines: 

^=^j d^Tr^ {W\t,Q)Et{t)T^HW{mE\mH) ; (32) 

W{t] 0) denotes a fundamental Wilson line running from t' ~ Q to t along the 
static trajectory of the heavy quark, are the generators of the gauge group in 
the representation of the heavy quark and du its dimension. In leading order this 
reduces to a Wightman-two-point function of fields at zero frequency, i.e., in the 
usual real-time propagator notation, 

K^^-Jj^pG> (^ = 0,p), (33) 

with Ch =4/3 the Casimir operator of the HQ representation. The integral is IR 
regulated by HTL corrections, i.e., a Debye mass, = g^T^{Nc + Nf/2)/3. In 
the left panel of Fig. |4] the NLO corrections to the LO result, Eq. ([33|) . arc de- 
picted in terms of Feynman diagrams. The double line represents the heavy quark, 
and all propagators and vertices include HTL corrections, leading to a gauge in- 
variant expression as it should be the case for an observable quantity like k. The 
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(C) {£>) 0.5 1 1.5 2 2.5 



9s 

Fig. 4. NLO calculations for HQ diffusion in the Left panel: NLO diagrams for the 

momentum-diffusion coefficient, k; the double line represents the heavy quark, all propagators are 
soft and HTL resummed, and all vertices include HTL vertices. Right panel: comparison of LO to 
NLO result for function of the strong coupling, aa. 

diagrams are evaluated in Coulomb gauge within the closed-time path (real-time) 
Kcldysh formalism of thermal quantum-field theory (TQFT). The real part of di- 
agram (A) provides a correction to the Debye mass. Diagrams (C) and (D) take 
into account real and virtual corrections by additional soft scattering or plasmon 
emission/absorption of the light or heavy scatterer, respectively. Diagram (B) rep- 
resents interference between scattering events occurring on the light scatterer's and 
on the heavy quark's side. Contrary to naive power counting, the NLO calcula- 
tion provides 0{g) corrections due to scattering with soft gluons with momentum, 
q c:i lio, and due to overlapping scattering events, dominated by t-ehannel Coulom- 
bic scatterings involving soft momentum transfers, ~ fijj oc gT. The right panel of 
Fig. |4] shows that the NLO correction to k is positive, i.e., the momentum-diffusion 
coefficient becomes larger compared to the LO calculation. The convergence is poor 
even for rather small coupling constants. A rigorous resummation scheme to cure 
this behavior is not known to date, especially to establish convergence in the typi- 
cal range of coupling con stants under conditions in relativistic heavy-ion collisions, 
as — 0.3-0.4. In Ref.^^ the investigation of NLO corrections is extended to the 
weak-coupling limit of A/" = 4 supersymmetric Yang-Mills (SYM) theory. Also in 
this case the perturbative series turns out to be poorly convergent, even for low 
couplings. 

2.2.5. Three-body elastic scattering 

Another step in the (would-be) perturbative hierarchy are three-body collisions, 
which are expected to become increasingly important at high parton density. An 
attempt to assess the effects of three-body elastic scattering for HQ diffusion has 
been conducted in Ref.[Z21, with pertinent Feynman diagrams as depicted in Fig. [5l 
Special care has to be taken in regularizing contributions from diagrams with in- 
termediate particles going on-shell; these can lead to divergent real parts in the 
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(a) 



(b) 



(c) 



Fig. 5. Different topological classes of diagrams for thr ee-b ody elastic scattering of a heavy quark 
(thick lines) off light quarks and antiquarks (thin lines|ZSl. 



scattering amplitude and represent successive tw o-b ody scatterings (rather than 
genuine three-body scattering). Therefore, in Ref.E3l 

the intermediate quark hnes 
in diagram (a) and (b) are supplemented with an in-medium collisional width, and 
only the real part of their propagator is kept in the evaluation of the diagrams. For 
three-body elastic processes involving one or two gluons, it has been assumed that 
the dominant contributions arise from diagrams with similar topology as diagram 
(b) in Fig. [5] for Qqq scattering; all other contributions are neglected. To compare 
with two-body gluo-radiative inelastic scattering the LO diagrams have been used 
to evaluate matrix elements for Qq Qqg, Qq Qqg and Qg Qgg processes. 
Within this scheme, at temperatures T = 200-300 MeV, three-body elastic scatter- 
ing processes are estimated to contribute to the c- and 6-quark friction coefficients 
with a magnitude comparable to two-body elastic scattering. Again, this raises the 
question of how to control the perturbative series for HQ diffusion. As a by-product, 
the friction coefficient for radiative scattering, Qp — > Qpg, was estimated to exceed 
the one from elastic two-body scattering for HQ momenta p > 12 GeV (for both 
charm and bottom). 



2.3. N on- Perturbative Interactions 

The evidence for the formation of a strongly coupled QGP (sQGP) at RHIC has mo- 
tivated vigorous theoretical studies of the possible origin of the interaction strength 
(see, e.g., 

Rcf.l24! 

for a recent review). In particular, several lattice QCD computa- 
tions of hadronic correlation functions at finite temperature have found indications 
that hadronic resonances (or bound states) survive up to temperatu res of twice 
the critical one or more (for both a gluon plasma (GP) and a QGpj3ll 42 | 75 | 76 |^ 
cf. also Sec. 14.1.11 of this article. Pertinent spectral functions (extracted from Eu- 
clidean correlators using probabilistic methods, i.e., the maximum entropy method) 
exhibit resonance peaks in both QQ and qq channels. The consequence s of hadronic 
resonances in the QGP fo r HQ tran sport have been elaborated in Refs.^^SEMTTllS] 
The starting point in 

Refs.ESEHIIII 

is the postulate that heavy-light quark [Qq) res- 
onances, i.e., "D" and "i?" mesons, persist in the QGP. In Refs.-^ this has been 
realized within an effective resonance model for Q-q scattering (Sec. I2.3.ip while in 
Ref.l^HQ fragmentation into mesons and their subsequent momentum broadening 
was considered (Sec. l2.3!3|) . The phenonicnological success of these models (cf., e.g.. 
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the right panel of Fig. I26p caUed for a more mic rosc opic evaluation of the heavy- 
light quark correlations. This was realized in Ref. 135] where in- medium heavy- light 
quark T-matrices were computed with interaction potentials estimated from HQ 
free energies in lattice QCD (cf. Sec. I2.3.2p . This approach is the direct analog to 
the potential models used in the heavy quarkonium context (cf. Sec. I4.1.ip . We 
finish this section with a brief discussion of a recent suggestion to extract informa- 
tion on HQ diffusion more directly from thermal lattice QCD (Sec. I2.3.4p . which 
would constitute a valuable benchmark for both perturbative and non-perturbative 
calculations. 



2.3.1. Effective Qq-resonance model 

The heavy-light quark resonance modePSl has been set up by combining HQ ef- 
fective theory (HQET) with chiral symmetry in the light-quark sector, q = (u,d), 
based on the Lagrangian, 



^Dcq + ^l, - iGs ( qn^c - qi'^^c + h. 



Gv ( qY'%^J^c - q^^Y'^,l±- 



(34) 



written in the charm sector (an equivalent one in the bottom sector follows via the 
replacements c — *■ 5 and D B for the HQ and resonance fields, respectively; v: 
HQ four- velocity). The pertinent free Lagrangians read 

■^c,q = c(i^ - mc)c + qi^q, 

$ and $g denote the pseudoscalar and scalar meson fields (corresponding to D and 
Dq mesons) which are assumed to be degenerate chiral partners (mass mg) as a 
consequence of chiral restoration in the QGP. The same reasoning applies to the 
vector and axialvector states (mass my), $* and (corresponding to D* and 
-D*). HQ spin symmetry furthermore asserts the degeneracy of spin-0 and -1 states 
with identical angular momentum, implying mg = my and the equality of the 
coupling constants, Gs = Gy. In the strange-quark sector only the pseudoscalar 
{Ds) and vector (-D*) resonance states are considered (i.e., chiral symmetry is not 
imposed). 

The boson-resonance propagators are dressed with heavy-light quark self- 
energies at the one-loop level (cf. left panel of Fig. [S]). To leading order in HQET, 
in accordance with spin symmetry, the self-energies for the vector /axialvector res- 
onances are given by 



n_D*,^ii. = {Vf^v^ - g^^)HD(s), 



(36) 
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Fig. 6. Left panel: one-lo op d iagram representing the D-meson self-energy in the QGP within the 
effective resonance modei25J. Right and middle panels: clastic Qq and Qq scattering diagrams for 
s- and u-channel resonance exchange, respectively. 

where s ~ denotes the meson's four momentum, and is the self-energy of the 
pseudoscalar/scalar resonances. Its imaginary part reads 

Imnz,(.) = -^ii^^e(. - ml) , (37) 

while the real part is calculated from a twice-subtracted dispersion relation with 
the wave-function and mass counter terms adjusted such that the following rcnor- 
malization conditions hold, 

a.4"')(s)|,=o = , Ren(;"')(s)|,=,„.^ = . (38) 
As an alternative regularization scheme, dipole form factors, 

have been supplemented to simulate finite-size vertices of the resonance model, 

lrnnf{s)=ln-,no{s)FWq\) , (40) 

with |q| = (s — m^)/(2-y/s). In this scheme, the real part is calculated from an 
unsubtracted dispersion relation, while the bare resonance mass is adjusted to obey 
the second renormalization condition in Eq. (|38p . 

With charm- and bottom-quark masses of TOc = 1.5 GcV and mf, = 4.5 GeV, 
the physical resonance masses are adjusted to niD = 2 GcV and ms = 5 GeV, 
respectively. This is in appr oximat e accordance with earlier T-matrix models of 
heavy-light quark interactionJ^SESl^ Likewise, the coupling constant. G, is adjusted 
such that the resonance widths vary as Tub = 0.4 . . . 0.75 GeV. The resulting 
heavy-light quark scattering matrix elements (cf. middle and right panels of Fig. [6]) 
have been injected into Eq. (fTS]) to calculate HQ drag and diffusion coefficients. In 
the left and right panel of Fig. [7] we compare the total HQ elastic scattering cross 
sections and resulting thermal relaxation times, Tcq = ^/A{p = 0), of the resonance 
model with LO pQCD (cf. the diagrams in Fig.[T]). Although the total cross sections 
are not very different in magnitude, the thcrmalization times decrease by around 
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Fig. 7. (Color online) Left panel: total HQ scattering cross sections off light partons in LO pQCD 
(blue lines) and within the effective resonance model (red lines). Right panel: thermalization times, 
T = 1/A(p = 0), for charm and bottom quarks in LO pQCD with as = 0.4 and Debye-screening 
mass ixo = gT, compared to the results from the resonance+pQCD model, as a function of QGP 
temperature. 



a factor of ~3-4 when adding resonant scattering, for all temperatures T = 1-2 Tc. 
The main reason for this behavior is that s-channel Qq scattering is isotropic in the 
rest frame of the resonance, while the pQCD cross section is largely forward-peaked 
(t-channel gluon exchange) , and thus produces a much less efficient transport cross 
section (which encodes an extra angular weight). The charm-quark equilibration 
times in the resonance+pQCD model, r^^ — 2-10 fm/c, are comparable to the 
expected QGP lifetime at RHIC of around tqgp — 5 fm/c. Thus, at least for charm 
quarks, substantial modifications of their pt spectra towards local equilibrium in 
the flowing medium can be expected. 

The consistency of the Fokker-Planck approach can be checked with the 
dissipation-fluctuation relation, Eq. ([77]) . at p = 0, cf. left panel of Fig. [H For the 
forward-peaked pQCD-matrix elements, the relation is fulfilled within 3%, while 
with the isotropic resonance scattering deviations reach up to 11% in the renor- 
malization scheme and up to 26% in the formfactor-cutoff scheme at the highest 
temperatures considered (T = 400 MeV). Note however, that for a typical ther- 
mal evolution at RHIC, average fireball temperatures above T = 250 MeV are only 
present within the first fm/tf^ below this temperature, the deviations are less than 
5% for all cases. The right panel of Fig. [8] illustrates that (for identical resonance 
widths) the formfactor regularization scheme leads to somewhat larger (smaller) 
friction coefficients at low (high) momentum than the renormalization scheme. 

In Fig. [9] the momentum dependence of the drag and transverse diffusion co- 
efficients is depicted using either resonance-scattering or pQCD-matrix elements. 
Resonance scattering becomes relatively less efficient for higher HQ momenta since 
the center-of-mass energy in collisions with thermal light antiquarks increasingly 
exceeds the resonance pole. The variations of the coefficients with the strong cou- 
pling constant in the pQCD scattering-matrix elements or the resonance-coupling 
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Fig. 8. Left panel: check of the dissipation-fluctuation relation at p = for c and b quarks with and 
without resonance interactions and in the renormahzation and form-factor cutoff schemes. Right 
panel: momentum dependence of the drag coefficient in the renormahzation and the formfactor- 
cutoff scheme, where the coupling constants have been chosen in both schemes as to obtain a 
resonance width T = 0.4 GeV. 
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Fig. 9. Drag (left panel) and transverse-diffusion coefficient (right panel) for pQCD and resonance 
model with varying interactions strengths as a function of the HQ momentum at a temperature 
of T = 200 MeV. 



constant in the effective resonance-scattering model are rather moderate. This is 
due to compensating effects of an increase of the matrix elements with or G*, 
on the one hand, and the accordingly increased Debye-screening mass for pQCD 
scattering or the broadening of the resonances widths, on the other hand. 



2.3.2. In-Medium T-matrix with IQCD-based Potentials 

The idea of utilizing HQ free energies computed in lattice QCD to extract a dri ving 
kernel for heavy-light quark interactions in the QGP has been carried out in Rcf.^^, 
with the specific goal of evaluating HQ diffusion. Since the latter is, in principle, de- 
termined by low-energy HQ interactions, the potential-model framework appears to 
be suitable for this task. Moreover, with a potential extracted from IQCD, the cal- 
culation could be essentially parameter-free. Currently, however, such an approach 
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bears significant uncertainty, both from principle and practical points of v iew, e.g., 
whether a well-defined potential description can be constructed in 

mcdiuninilM] 

and, if so, how to extract this information from, say , the H Q free energy. In the 
vacuum, both questions have been answered positivel y I I, thus validating the 30 
year-old phenomenological approaches using Cornell potentials for heavy quarko- 
nia, which provide a very successful spectro scopj^ . The potential approach has 
been extended to heavy-light mesons in Refs.l^^l^. 

A Brueckner-likc in-medium T-matrix approach for heavy-light quark scattering 
in the QGP has been applied in Ref.S^, diagrammatically represented in Fig. [TUl 
The underlying (static) two-body potential has been identified with the internal 
energy 

Ui (r, T) = F, (r, T) + TS, (r, T) = F, (r, T) - T^^^f^ , (41) 

extracted from two IQCD computations of the color-singlet HQ free energy above 
Tc, for quenchecP^ and two-flavoJ^ QCD (pertinent parameterizations arc given 
in Refs.^^ andSSl, hereafter referred to as [Wo] and [SZ], respectively). This choice 
(rather than, e.g., the free energy) provides an upper limit for the interaction 
strengtlilmSHl. 

To use Eq. (j4ip as a potential in a T-matrix calculation, the 
internal energy has to be subtracted such that it vanishes for r — > oo, 

yi(r,T) = C/i(r,T)-C/i(r^oo,r) , (42) 



which is dictated by the convergence of the T-matrix integral in momentum space. 
It is suggestive to interpret the asymptotic value C/f° = Ui{r ^ oo,T) as an in- 
medium HQ-mass, 

"^Q(T)-mO, + ic/- , (43) 

where itiq denotes the bare HQ mass (e.g., m° ~ 1.25 GeV^^for the bare c-quark 
mass). However, close to Tc, the values for Ut ^(T) e xtracted from IQCD calculations 
develop a rather pronounced peak structurJ^SlMl^ which renders a mass interpre- 
tation problematic. Progress in understanding these properties is closely connected 
with the proper identification of the potential. First IQCD estimates of the in- 
medium HQ mass (extracted by relating zero-mode contributions to quarkonium 
correlators to the HQ su sce ptibility) indicate a moderate increase when approaching 
Tc from abovcHS] in Ref.H^ constant (average) in-medium charm- and bottom-quark 
masses of nic — 1.5 GeV and nit = 4.5 GeV, respectively, have been employed. 

The T-matrix approach is readily generalized to color-configurations other than 
the singlet channel of the qQ pair. The complete set of color states for Qq (sin- 
glet and octet) and Qq (anti-triplet and sextet) pairs has been taken into account 
assuming Casimir scaling as in pQCD, 

Vs^^lvi, V-, = ^Vu Ve = -\Vi. (44) 
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Fig. 10. (Color online) Diagrammatic representation of the Brueckner many-body calculation for 
the coupled system of the T-matrix based on the IQCD static internal potential energy as the 
interaction kernel and the HQ self-energy. 



which is, in fact, supported by finitc-T To augment the static (color- 

electri c) p otentials with a minimal rclativistic (magnetic) correction for moving 
quarkJS^ ^j^e so-called Breit correction as known from electrodynamicfP^ has been 
implemented via the substitution 

Va^Vail~ai-a2), (45) 

where ai_2 are quasiparticle velocity operators. 

The above constructed heavy-light potentials can now be resummed in a two- 
body scattering equation. In accordance with the static nature of the potentials, it is 
appropriate to use a three-dimensional reduction of the full four-dimensional Bethe- 
Salpeter equation. This leads to the well-known ladder series which is resummed 
by the Lippmann-Schwinger (LS) integral equation for the T-matrix In the q-Q 
center-of-mass (CM) frame it takes the form 

TaiE; q', q) =Va{q\ q) - j 7^K(g', k) G,q{E; k) 



xTa{E;k,q) [1 - f,{col) - fgioj'^)] 



(46) 



The driving kernel (potential) can now be identified with the Fourier transform of 
the coordinate-space potential extracted from IQCD, 



Va{q\q) = j d?rVa{r) exp[i(q - q') ' ^1 



(47) 



in a given color channel, a G {1, 3, 6, 8}. The concrete form of the intermediate q-Q 
(or q-Q) propagator, GqQ{E, k), depends on the reduc tion scheme of the underly- 
ing Bethe-Salpeter equa tion. It has been verifiecPIESl that, e.g., the Thompson^ 
and Blancenbecler-SugaA^ scheme lead to very similar results in the present con- 
text (as was found for nucleon-nucleon scattering). In the former, the two-particle 
propagator is given by 

G,q{E, k) = ^ . „ . .J^ , a . .^o. > (48) 



i?-(c.^ + iS|)-(. 
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where E and k denote the CM energy and relative momentum of the qQ pair, 
respectively. The quasi-particle widths are chosen as Tf^ = 2^'^ = 200 MeV, 
and the light quark masses as constant at — 0.25 GcV, with on-shell energies 



^r=^J<Q + k'- (49) 
The latter figure into the Pauli blocking factor with Fcrmi-Dirac distributions, 

= exp(J/T) + l 

(at the considered temperatures their impact is negligible). The solution of the 
T-matrix Eq. (|46|) is facilitated by a an expansion into partial waves, I, 



K(g',q) =47r^(2? + l) K,i(g',g) P/[cosZ(q,q')] 



Ta{E; q', q) = 4^ ^(2/ + 1) Ta^E] q' , q) P/[cos Z(q, q')] 
I 



(51) 



which yields a one-dimensional LS equation, 

Ta.i{E-q\q) = K,i(g',g) + ^ J dkeVaAq'^k) Gq^iE-k) 

xTa,i{E-k,q) [1-Jf{o^2)^If{^1)] , 



(52) 



for the partial- wave components, Ta,;, of the T-matrix. Eq. ()52|) can be solved 
numerically by discret izati on and subsequent matrix-inversion with the algorithm 
of Haftel and TabakiJEH. 

The resulting 5- wave (Z = 0) T-matrices indeed show 
resonance structures in a QGP in the channels where the potential is attractive, i.e., 
in the meson (color-singlet) and diquark (color-antitriplet) channels. The pertinent 
peaks in the imaginary part of the T-matrix develop close to the Q-q threshold, and 
melt with increasing temperature at around 1.7 Tc and lATc, respectively (cf. left 
panel of Fig. [TT|) . In the repulsive channels, as well as for P- waves, the T-matrices 
carry much reduced (non-resonant) strength. The increasing strength in the meson 
and diquark channels (the latter relevant for baryon binding) when approaching Tc 
from above is suggestive for "pre-hadronic" correlations toward the hadronization 
transition. 

The next step is to use the T-matrices to compute the light-quark contribution 
to HQ self-energies, i.e., the last diagram in the second line of the Brueckncr scheme 
illustrated in Fig. 1101 In a given color-channel, a, the T-matrix induced self-energy 
is given by 

S?(c.,p) = ^ J^^[M^,) + f^ic.+u;,)]UE;p,k) , (53) 

where dsi = 4(2Z + 1)A^/ denotes the spin-isospin and angular momentum degen- 
eracy of all Qq (or Qq) configurations (assuming spin and light-flavor symmetry) 
and da the color degeneracy of channel a; the factor 1/6 averages over the incoming 
HQ color-spin degrees of freedom. The resulting charm-quark selfenergies (summed 
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Fig. 11. (Color online) Results of a Brucckner-type approach for c-quarks in a base d on 

a potential corresponding to the internal energy extracted from quenched lattice QCeEHEQI. Left 
panel; imaginary part of the in-medium T matrices for S-wave c-q scattering in color-singlet and 
-triplet channels at two different temperatures; right panel: real and imaginary parts of c-quark 
selfenergies based on the T-matricos in the loft panel. 



over all light quarks and antiquaries) arc displayed in the right panel of Fig. 1111 One 
finds rather small corrections to the HQ mass (presented by the real part of S), 
but the imaginary parts are substantial. F^ = — 2ImI]c ~ 100-300 McV for temper- 
atures T = 1.1-1.8 Tc (with the largest values attained close to Tc). These values 
were the motivation for the choice of input widths in the propagator. Eq. (|48|) . 
of the T-matrix equation, thus providing a rough self-consistency. The in-medium 
mass corrections, on the other hand, are associated with the gluonic contribution 
to the HQ self-energy (corresponding to the first term in the lower line Fig. [TO)) . 
which have not been calculated explicitely in Ref. SSl^ represent a rough (aver- 
age) representation of the asymptotic values of the HQ potential, J/i^'' (as discussed 
above). 

The final step is to implement the T-matrix elements into a calculation of HQ 
drag and diffusion coefficients via Eq. (jl5p : one finds 

= ^{s -ml+ mlf{s -ml + mlf 

XiV/^rfa(|Ta,/=o(s)|'+3|ra.i = l(s)cOs(0cm)n • ^^^^ 

a 

The resulting friction coefficients are summarized in Fig. [T^] as a function of mo- 
mentum for three temperatures and for two potential extractions from IQCD | 44 | 90 | 
Generally, the Qq T-matrix based coefficients are largest at low HQ momentum, as 
to be expected from the resonance formation close to threshold. The values exceed 
the LO-pQCD coefficients at small temperatures and for both potentials by a factor 
of ~3-5. At higher temperatures the enhancement reduces considerably, to a factor 
of less than 2 in the [SZ] potential and to essentially equal strength for the [Wo] 
potential^. In fact, the coefficients computed with the [SZ] potential have a slightly 



"^Recall that the LO-pQCD calculations employ a rather large coupling (as =0.4), and are domi- 
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Fig. 12. (Color online) The drag coefficients at different temperatures, using the parameterization 
of the HQ potential from [Wo] (left panel) and [SZ] (right panel) compared to LO pQCD with 
cis = 0.4 and = gT. 



(30%) increasing trend with T, while the [Wo] potential leads to a decreasing trend. 
This difference needs to be scrutinized by future systematic comparisons of IQCD 
input potentials. Compared to the resonance model (cf., e.g., left panel of Fig. [9]), 
the T-matrix calculations yield quantitatively similar results at temperatures not 
too far above Tc but become smaller at higher T due to resonance melting (which 
is presumably a more realistic feature of a non-perturbative interaction strength). 



2.3.3. Collisional Dissociation of Heavy Mesons in the QGP 

In Ref a so-called reaction-operator (GLV) approach has been applied to resum 
multiple elastic scatterings of a fast Qq pair. The quenching of heavy quarks in 
a QGP is calculated by solving coupled rate equations for the fragmentation of c 
and b quarks into D and B mesons and their dissociation in the QGP. The main 
mechanisms for HQ energy loss are collisional broadening of the meson's transverse 
momentum and the distortion of its intrinsic light-cone wave function. The latter is 
modeled in a Cornell-type potential ansatJ^ (note that this bears some similarity 
to the T-matrix approach discussed in Sec. 12. 3. 2|) . This in-medium HQ/heavy- meson 
fragmentation/dissociation mechanism leads to comparable liigh-pt suppression for 
B and D mesons, which is quite contrary to perturbative calculations for both 
collisional and radiative energy loss (where the suppression of b quarks is signif- 
icantly less pronounced than for c quarks). This feature largely results from the 
much smaller formation times of B-mesons compared to D- mesons, leading to a 
faster fragmentation-dissociation cycle for b quarks/i? mesons. 



nated by scattering off gluons in the heat bath; thus a minimal merging of the gluon sector with 
the T-matrix calculations would consist of adding the gluonic part of LO-pQCD; this procedure is 
adopted below whenever combined results are shown or utilized. In principle, a non-perturbative 
treatment should also be applied to HQ-gluon scattering. 
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2.3.4. Estimates of HQ Diffusion in Lattice QCD 

It has recently been suggested that, unUke in the case of other transport coefficients 
(e.g., the shear viscosity), the HQ diffusion coefficient might be amenable to a 
determination within IQCD, based o n an analytic continuation of the color-electric- 
field correlator along a Polyakov loo]Jl^ (see also Ref.^^^for earlier related work). 
The starting point of these considerations is the spectral function of the HQ current 
correlator, 

p^^'icu) = J dicxp(ic^t) J d^a;^! [j^{t,x), J^{0,0)\^ , (55) 

where Jq denotes the HQ current operator in the Heisenberg picture. The spatial 
diffusion coefficient, Dg, can be extracted from this spectral function by the pole 
position at cij = —iDgk^, where k is the HQ momentum. The condition for a pole 
leads to the Kubo relation 

As = :^limE^, (56) 
1=1 

where xoo is the conserved-charge susceptibility 

X°° = M d'x{j'Qit,x)j'QiO,0)) . (57) 



T „ 

For a heavy quark, the spectral function, Eq. is expected to develop a sharp 
peak around w = which can be described by a Lorentzian function close to this 
point. The width of this function is given by the drag coefficient, which obeys the 
fluctuation-dissipation relations, discussed in Sec. 12.11 Using HQ effective theory 
techniques it is shown that in the static limit the momentum-diffusion coefficient, 
K = 2D, is given by a correlator involving color-electric fields and Jq operators 
whose Euclidean analogue can be mapped to an expectation value involving Wilson 
lines and color-electric fields, similar to Eq. ([32]) . This purely gluonic correlation 
function can in principle be evaluated in IQCD. 

Anot her lattice-based approach to assess HQ diffusion has been suggested in 
Ref.™ in terms of (discretizcd) classical gauge theory. The limitation of this ap- 
proach is set by the thermal (hard) scale ~7rT where quantum theory suppresses 
excitations. Since the HQ thermalization rate is expected to be governed by the 
electric screening scale ^gT, this limitation may not be severe for small and mod- 
erate coupling {gT ^ T), and thus allow for valuable insights. First, it has been 
verified that, in the weak coupling limit, the discretized (nonperturbative) classical 
computation indeed agrees well with pQCD. Upon increasing the coupling strength, 
the classical lattice results for the HQ momentum diffusion coefficient increasingly 
exceed the LO perturbativc result, by about an order of magnitude for a moderate 
coupling strength corresponding to as — 0.2. Next, the NLO term (with slightly 
increased strength to account for HTL effects), as calculated in Ref.^^, has been 
added to the LO calculation which extends the agreement of pQCD with the classi- 
cal lattice results to larger (but still weak) coupling. Forofs ~ 0.2 the increase over 
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LO amounts to a factor ^-^2, which means that the classical lattice result remains 
substantially larger (by a factor of ^^5) than the NLO value. Besides reconfirming 
the poor convergence of pQCD, this also suggests that the perturbative series is 
not alternating but that higher order terms keep increasing the value of the HQ 
momentum diffusion coefficient. Semi-quantitatively, such an enhancement is in the 
ball park of the factor ^3-4 found in the effective resonance model (Sec. 12.3. ip or 
T-matrix approach fSec. 12.3?^ . 

2.4. String Theoretical Evaluations of Heavy-Quark Diffusion 

The conjectured correspondence between certain classes of string theories, formu- 
lated in five-dimensional Anti-de-Sitter space (AdSs), and gauge theories with con- 
formal invariance (conformal field theory, CFT) has opened interesting possibilities 
to address nonperturbative aspects of QCD. The so-called AdS/CFT correspon- 
dence implies a "duality" of a weakly coupled gravity to a strongly coupled super- 
symmetric (and conformal) gauge theory, specifically N — A SU(Afc) super- Yang- 
Mills (SYM) theory. This connection has been exploited to formulate the problem of 
HQ diffusion at finite temperature and extract m \ "exact" no nperturbative result 
for HQ transport coefRcients in the SYM plasmiMMUnZ]. 

The translation to 

QCD matter is beset with several caveat a I, e.g., the particle content of the SYM 
medium is quite different compared to the QGP. While this may be corrected for 
by a suitable rescaling of the temperature by matching, e.g., the energy densitie^, 
a more problematic difference is the absence of a scale (other than temperature) 
in conformal SYM. Thus, the latter does not possess a breaking of scale invari- 
ance, a running coupling constant, confinement nor spontaneous chiral symmetry 
breaking, and consequently no notion of a critical temperature, either. Thus SYM 
is quite different from QCD in the zero- and low-temperature regimes. However, 
at sufficiently high T, where the QCD medium deconfines its fundamental charges, 
the resemblance to SYM might be much closer. E.g., the pressure in SYM in the 
strong coupling limit amounts to about 75% of the Stefan-Boltzmann limit, close 
to what is foimd in thermal lattice QCD for a wide range above Tc. In addition, the 
finding of an extremely small shear viscosity in s trongly coupled SYM, rj/ s = l/Air 
(conjectured to be a universal lower boundjl^, and the apparently low-viscosity 
QCD medium deduced from the success of hydrodynamic models at RHIC, is an- 
other good reason to further pursue exact nonperturbative calculations in SYM for 
quantities that are relevant for RHIC phenomenology. If nonperturbative effects in 
the strongly coupled QGP at moderate temperatures, T = 1-2 Tc, are ultimately 
connected to the presence of the phase change(s) (and thus inherently to the crit- 
ical temperature as a relevant scale), the CFT-QCD connection would not be a 
rigorous one. But even in this case, the nonperturbative computation of transport 

"This procedure works quite well when comparing quantities in quenched and unquenched lattice 
QCD computations, e.g. for the critical temperature. 
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coefBcicnts of a strongly coupled system at a given reference temperature "not too 
close" to Tc should provide useful insights. 

The first step in computing HQ diffusion for CFT is the introduction of a heavy 
quark into the conformal field theory. This can be achieved by either introducing 
a heavy charge via breaking the gauge group from + 1 to Nc (which, strictly 
speaking, generates {2Nc + 1) "Higgsed" 'W bosons or by adding a finite- 
mass Af = 2 hypcrmultiplct with charges in the fundamental representation as 
a "probe" of the CFT medium. In either case, the pertinent object on the 4- 
dimensional b oundary of the 5-dimcnsional AdS space represents a fundamental 
charge. In Refs.^^l ^^^ [ the HQ drag has been evaluated by computing its momen- 
tum degradation, dp/dt = — 7P, through the force on the trailing string, resulting 
in a friction (or drag) coefficient, 

where A = 5sym^c denotes the 't Hooft coupling constant. Alternatively, in Ref.lin71 
the problem was formulated focusing on the diffusion term. For time scales longer 
than the thermal relaxation time of the medium, but short compared to the HQ 
relaxation time, the fluctuation term in the Langcvin equation (j63p dominates over 
the drag term. The evaluation of the noise (or force) correlator is then carried out 
via the fluctuations of the string, resulting in a noise coefficient which is directly 
related to the diffusion coefficient (cf. Eqs. ([S5|) and ([71]) below). Furthermore, the 
latter can be related to the friction coefficient using the Einstein relation, Eq. (|22p : it 
turns out that the result is identical to Eq. (j58p . which also verifies that a Langevin 
process consistent with the fiuct uation-dissipation theorem applies in the SYM the- 
ory (see, however, Refs.m ^l^^^[ where the applicability of the Langevin framework 
in AdS/CFT for high- momentum quarks is discussed). The square-root dependence 
of 7Ads/CFT on the coupling constant A clearly characterizes its nonperturbativc 
nature; in this sense it is parametrically large for comparatively small coupling 
constants. The temperature dependence is rather "conventional" , as to be expected 
since there are no additional scales in the problem (the HQ mass in the denominator 
implies the standard suppression of the HQ relaxation rate by ^T/mg). 

The next question is how to convert the result into a (semi-) quantitative es- 
timate for the QCD plasma. Naively, one may just insert the values of the strong 
coupling constant, gs, and QGP temperature, T, for gsYM and Tsym, respectively. 
A more suitable identification probably consists of matching p hysic al quantities 
which leads to somewhat different parameter values. E.g., in Ref.™ 

comparable 

temperatures were identified by matching the energy densities (e) of the QGP and 
SYM-plasma. Since the latter has a factor ~3 larger particle content (degeneracy 
factor), one has a smaller temperature at the same e, Tsym — For the cou- 

pling constant, one can exploit the fact that in AdS/CFT the potentia l bet ween a 
heavy charge an d anticha rge is essentially of Coulomb-type, both at zer cff^ and fi- 
nite t emp er at ur (j ^ ^ I \ 2M tj^g latter case, the potential goes to zero at some finite 
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range, characteristic for Debye-scrcening behavior. This range can be used to iden- 
tify the length scale in comparison to typical screening radii of heavy-quark fre e en- 
ergies as computed in thermal lattice QCD (although some ambiguity remains j^^. 
Matching the magnitude of the potentials at the screening radius then allows for 
a matching of the coupling constants. This leads to significantly smaller values for 
A (by a factor of 3-6) than the naive identification with Us = 0.5. In connection 
with the redefined temperature, the improved AdS/CFT-based estimate for the HQ 
friction coefficient in QCD amounts to 7 ~ 0.3-0.9 c/fm at T 250 MeV, which is 
significantly smaller than the "naive" estimate of ^2 c/fm. 



2.5. Comparison of Elastic Diffusion Approaches 

In view of the recent proliferation of seemingly different approaches to evaluate HQ 
transport coefficients in the QGP it becomes mandatory to as k to what extent they 
are related and encode similar microscopic mechanism j ^^^[ It turns out that all 
of the approaches discussed above incorporate a color-Coulomb-type interaction. 
This is rather obvious for the T-matrix approach, where the input potentials from 
lattice QCD clearly exhibit the Coulomb part at sufficiently small distance (in- 
cluding effects of color screening). The one-gluon exchange in pQCD (which is the 
dominant contribution to HQ rescattering, recall the two right diagrams in Fig.[T]), 
also recovers the Coulomb potential in the static limit (color screening enters via 
the Debye mass in the spacelike gluon-exchange propagator). The collisional dis- 
sociation mechanism involves the Cornell potential for the D- and i?-meson wave 
functions and thus incorporates a Coulomb interaction as well; the emphasis in 
this approach is on formation-time effects essentially caused by the different (free) 
binding energies of D and B mesons. In addition, the confining part of the Cornell 
potential may play a role (as in the T-matrix approach) . Finally, in conformal field 
theory (AdS/CFT), the absence of any scale promotes the Coulomb potential to the 
unique form of a potential, V{r) cx 1/r (this is the only way of generating a quantity 
with units of energy). On the other hand, scale-breaking effects are present in the 
QCD-based approaches in terms of a running coupling constant (pQCD), while the 
Cornell and IQCD-based potentials additionally feature linear terms oc or wh ere 
the string tension introduces a further (nonperturbative) scale. In fact, in Rem 
it has been argued, based on an analysis of IQCD results for the heavy-quark free 
(and internal) energy, that "remnants" of the confining force play a prominent role 
for temperatures not too far above Tc (e.g., for heavy quarkonium binding). 

If one assumes the prevalence of the Coulomb interaction, the obvious first 
question is with what strength (coupling constant) it figures into the different ap- 
proaches, which should be fairly straightforward to determine. A more involved 
issue is to scrutinize the underlying approximation schemes and their applicability. 
E.g., perturbative approaches with large (running) coupling constants have poor 
(if any) control over higher-order corrections. As usual in such situations, diagrams 
with large contributions should be identified and resummed (which is, of course, a 
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Fig. 13. (Color online) Charm-quark friction coefficients, 7, in the QGP. Left panel: three- 
momentum dependence at three temperatures (color code) for: LO-pQCD with fixed as = 0.4 
and un = gT (dash-dotted lines), heavy-light quark T-matrix plus LO-pQCD for gluons (solid 
lines and pQCD with running as and reduced infrared regulator (dashed lines jSfiJfiZI. Right 
panel: temperature dependence of 7 for LO-pQCD, T-matrix plus LO-pQCD (gluons only), pQCD 
with running as, and from AdS/CFT correspondence matched to QClJ '^^'^ l with C = 1.5-2.6 ^^^ 1 



non-trivial task, e.g., maintaining gauge invariance); it would be illuminating to ex- 
tract a static gluon-exchangc (Coulomb) potential for a given set of parameters. The 
T-matrix approach performs a resummation of the ladder series of a static (color- 
electric) potential; magnetic interactions are implemented in a simplified manner 
using the Breit current-current interaction from electrodynamics. It has been veri- 
fied that for large center-of-mass energies, the qQ T-matrix recovers the result for 
perturbative scattering. However, a number of effects are neglected and need to 
be scrutinized, including the interactions with gluons beyond pQCD, retardation, 
extra gluon or particle/antiparticle emission (e.g., in a coupled channel treatment) 
and the validity (and/or accuracy) of a potential approach at finite temperature 
(this issue will reappear in the context of heavy quarkonia in Sec. |4]). In the coUi- 
sional dissociation approach, it would be interesting to explore medium effects in the 
employed potential (i.e., on the mesonic wave function). Ideally, by improving on 
specific assumptions in a given approach, an agreement would emerge establishing 
a common result. Explicit connections with the AdS/CFT results are more difficult 
to identify. Maybe it is possible to push the T-matrix approach into a regime of 
"large" coupling, or study the existence and properties of {D and B) bound states 
in the string theory setting. In Fig. [13] we summarize the drag coefficients as func- 
tion of momentum (for three temperatures, left panel) and temperature (for p ~ 0, 
right panel) resulting from the approaches discussed above, i.e.. 



(i) leading-order pQCD calculations with fixed as = 0.4 and Debye-screcning mass, 
Hd = gT, in the t-channel gluon-exchange contributions to the matrix elements 
for elastic gQ and qQ scattering, 

(ii) in-medium T-matrix calculations using IQCD-based qQ potentials, augmented 
by the leading-order pQCD matrix elements for elastic gQ scattering^, 
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Fig. 14. (Color online) Spatial diffusion coefScient, Ds = T/{'ymQ), for c (left) and b quarks (right) 
in a QGP for: LO-pQCD with fixed et .. = 0.4 (dashed lines), effective resonance model + LO- 
pQCD (bands for Tu,B = 0.4-0.75 GeV)2 3, T-m atrix approach + LO-pQCD (gluons only p51 and 
pQCD with running as (dash-dotted linejSSEZl. The AdS/CFT result corresponds to 2-kTDs = 
27r/C ~ 1.5 — 4 (not shown in the plots). 



(iii) pQCD calculations with running as and reduc ed scr eening mas j^^ l ^ , and 

(iv) the AdS/CFT correspondence matched to QCEffSSlwith 7qcd = CT'^/mq for 
C = 1.5-2.^1151 . 



At all temperatures, the T-matrix approach, (ii), produces significantly more HQ 
interaction strength than LO pQCD, (i), while for T> 0.2 GeV the thermaliza- 
tion rate for the T-matrix is a factor of ^-^2-4 less than for AdS/CFT, (iv), or for 
LO-pQCD with running coupling and reduced infrared regulator, (iii). Close to 
Tc~ 180 MeV, however, the three approaches (ii), (iii) and (iv) are not much dif- 
ferent and share overlap around 7~ 0.2c/fm. The spread in the numerical results 
reiterates the necessity for systematic checks as indicated above. 

Finally, one can convert the drag coefficients into estimates of other HQ trans- 
port coefficients of the QGP. Within the Fokker-Planck approach the spatial dif- 
fusion coefficient, Ds, is directly related to the drag coefficient, 7, as given by 
Eq. ([24]) . Fig. [14] shows the dimensionless quantity 2ttTDs for charm (left panel) 
and bottom quarks (right panel) as a function of temperature for LO pQCD, LO 
pQCD with running coupling and reduced infrared regulator, effective resonance 
model and T-matrix approach. The former three are fairly constant as a function 
of temperature while the T-matrix approach exhibits a significant increase with 
temperature, indicating maximal interaction strength close to Tc. This originates 
from the increasing potential strength (decrease in color-screening) with decreasing 
temperature, enhancing resonance correlations at lower temperature. It is tempting 
to interpret this feature as a precursor phenomenon of hadronization. However, its 
robustness needs to be checked with a broader range of lattice potentials. We recall 
that the internal-energy based potentials probably provide an upper estimate for 
the strength of the interaction. It is interesting to note that for all approaches the 
results for b quarks coincide with the ones for c quarks within ^20-30%. The largest 
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deviation is seen in the T-matrix approach, where the (spatial) diffusion coefhcient 
is smaller for b quarks than for c quarks (B-meson resonances survive until higher 
temperatures than D resonances). This is qualitatively similar to what has been 
found for the coUisional dissociation mechanism, where the relative enhancement of 
the 6-quark energy loss (compared to charm) is due to smaller i?-meson formation 
times. Since the latter arc related to larger i?-meson binding energies, the dynam- 
ical origin of the smaller Ds for b quarks appears to be of similar origin as in the 
T-matrix approach. The relative magnitudes of the various approaches reflect what 
we discussed before for the drag coefficient. 

2.6. CoUisional vs. Radiative Energy Loss 

For slowly moving heavy quarks in the QGP, the parametrically dominant inter- 
action is elastic scattering. However, at high px, radiative scattering is believed to 
eventually become the prevailing energy-loss mechanism. It is currently not known 
at which px this transition occurs. Therefore, it is important to assess the relative 
importance between elastic and inelastic scattering processes in the medium, even 
at the level of perturbative scattering only. Toward this purpose, we first recollect 
basic results on the gluon-Bremsstrahlung mechanism for light-parton, and then 
HQ, energy loss in the QGP, followed by a direct comparison to coUisional energy 
loss for heavy quarks. 

A seminal perturbative treatment of gluo-radi ative ene rgy loss (E-loss) of high- 
energy partons in the QGP has been given in Refs.EMIIl(BDMPS). The medium 
is modeled as static scattering centers which implies that the E-loss is purely ra- 
diative. The key finding is that the E-loss due to multiple in-medium scattering 
of a high-energy parton grows as L^, where L is the path length of the parton 
traversing the medium. The static scattering centers, at positions described 
by screened Coulomb potentials, 

V^{Q) = 2 ^ 2 exp{-iqxi) . (59) 

The range of the potentials is assumed to be small compared to the mean free path, 
A, of the scattered parton, i.e., 1/ fio <C A. In this case successive scatterings can be 
considered as independent, thus enabling an eikonal approximation for the elastic 
scattering on static centers, i.e., a classical propagation of the particle with energy 
E ^ ^D, undergoing independent kicks, thereby radiating Bremsstrahlung gluons. 
In analogy to the QED case an important ingredient is the coherent resummation of 
the multiple-scattering Bremsstrahlung amplitudes ( "Landau-Pomeranchuk-Migdal 
effect") which can be formulated as a diffusion equation for the effective scatter- 
ing amplitudes (or pertinent currents) . The total radiative E-loss of a high-energy 
parton traversing a medium of path length L is then given by 



(60) 
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where q is the diffusion coefBcicnt for transverse-momentum broadening in scat- 
tering off the static scattering centers, (q^) = qL. Perturbative calculations of the 
transport coefficient result in a value of about g ~ 1 GeV/fm^ at typical energy 
densities of e ~ 10 GeV/fm'"' (translating into T ~ 250 MeV) relevant for the QGP 
at It turns out, however, that the description of high-momentum pion 

suppression at RHIC in the BDMPS formalism requires an approximate ten-fold 
increase of the perturbative value for Recent calculations of perturbative E- 
loss including both elastic and radiative contributions within a thermal-field-theory 
framework indicate that coUisional E-loss may be significant even for high-p^ light 



partonJ2QI. This would imply a reduction of the value required for q from RHIC 
phenomenology. 

An early calculatiorP^ of radiative charm-quark E-loss, — d£'/da-, in the QGP 
has found that it dominates over the elastic one down to rather small momenta, 
p < 2 Ge\^21] In Ref.EZl it has been pointed out that the application of radiative 
E-loss to heavy quarks leads to the appearance of the so-called "dead cone", i.e., 
the suppression of forward gluon radiation for 9 < toq/£', where 9 denotes the 
direction of motion of the gluon with energy i5, relative to the direction of the HQ 
momenturrPSl. It has been predicted that the reduced HQ E-loss leads to a heavy- 
to-light hadron ratio above one in the high-pT regime accessibl e at RHIC . Within 
the BDMPS model, extended to heavy quarks, it has been arguecP^, however, that 
medium-induced gluon radiation tends to fill the dead cone. As will be discussed in 
Sec. 13.51 a similar value for q as in the light-hadron sector is necessary to come near 
the observed suppression of high-pT electrons from HQ decays in terms of radiative 
E-loss alonJSSl. 

The BDMPS formalism for light partons has been generalized to resum an ex- 
pansion of gluo-radiative parton E-loss in the G P with opacity, n = i/A, employ- 
ing a so-called reaction operator approaclff^ (GLV). A reaction operator i?„ is 
constructed that relates the n**^ power in a opacity-inclusive radiation probability 
distribution to classes of diagrams of order n—1. This results in a recursion relation 
for the radiation probability distribution, corresponding to a certain resummation 
to all orders in opacity, which can be implemented in Monte-Carlo simulations for 
jet quenching. The GLV reaction-operator method for light -parton radiative E-loss 
in the QGP has been extended to heavy quarks in (DGLV), implementing 

the kinematical suppression of gluon radiation by the HQ mass in the "dead cone" . 

A direct study of the relative magnitude of coUisional (elastic) and radiative 
pQCD HQ E-loss in the GP has been undertaken in Ref.ESl. For the elastic E-loss 
of a parton with color Casimir constant, Cr, the leading logarithm expression in 
an ideal QGP with Nf effective quark fiavors at temperature T, 

'^^'^ CniralT^ (l + ^\ f{v) In(Bc), (61) 



dx " ' V 6 



has been used. In an ultrarelativistic gas of massless partons the jet-velocity function 
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Fig. 15. (Color online) Loft panel: average relative E-loss, AE/E, for u, c and ,b quarks as a 
function of jet energy, E, in a longitudinally (Bjorken) expanding QGP, with fixed path length 
L = 5fm, initial gluon rapiditiy density dNg/dy = 1000 and fixed Os = 0.3; the gluon (light-quark) 
mass is set to fj-n/V^ (A'-d/2), the c(6)-quark mass to rric = 1.2(4.75) GeV (solid lines: radiatve 
E-loss, dashed bands: elastic E-loss in two schemes as discussed in the text). Right panel: parton 
nuclear modification factor, Ra_a=Rq, for gluons, u-, c-, and 6-quarks as a function of pT for a 
fixed path length and dNg/dy = 1000 (dashed lines: radiative E-loss, solid linos: radiative-|-clastic 
E-loss). 

is given by 



while estimates for are taken from Refs. [Bjp^Sl, [TGp^, and [BTpEi. The 
different values for Be obtained in these models are considered as reflecting theo- 
retical uncertainties. The r adiat ive E-loss with in the D GLV reaction-operator ap- 
proach is calculated in 

Ref.lMI based on Refs! 121 | 125 | r^.^^ jgf^ ^^^^^i pjg ^ 
compares pQCD radiative and coUisional E-loss for various quark flavors (masses) 
at high PT > 5 GeV in a gluon plasma (GP) with T ~ 240 MeV. For light and 
charm (bottom) quarks the elastic E-loss is comparable to the radiative one up 
to pt — 10(20) GeV, and still significant above. The right panel of Fig. [15] reiter- 
ates that, within pQCD, coUisional E-loss is an essential component in calculating 
the suppression of light-parton and especially HQ spectra at RHIC. Note that the 
relative importance of coUisional E-loss is expected to increase if non-perturbative 
effects become relevant (which predominantly figure toward lower pr), or if the GP 
is replaced by a QGP. 

2.7. D Mesons in the Hadronic Phase 




(62) 



To complete the discussion of open charm in QCD matter we briefly address medium 
modifications of charm hadrons in hadronic matter. Pertinent studies may be di- 
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vidcd into calculations for cold nuclear matter as well as for hot meson matter. 

Early studies of D-mesons in cold nuclear matter focused on possible mass 
shifts due to scalar and vector mean fields acting on the light-quark content of 
the meso At normal nuclear matter density ^jat = = 0.16fm ^, attractive 
mass shifts of up to — lOOMeV have been reported for and mesons (where 
both mean fields contribute with the same sign) while the mass change of the 
and turned out to be small due to a cancellation of the mean fields. Similar 
findings have been reported in QCD sum rule calculation^^ where the (isospin- 
averaged) D-meson mass is reduced by about — 50MeV, mostly as a consequence of 
the reduction in the light-quark condensate. Rather different results are obtained 
in microscopic calculations of ZJ-meson selfenergies (or spectral f unctions] based 
on coupled channel T-matrices for DN scattering in nuclear mattci^lSHIES]^ These 
calculations incorporate hadronic many-body effects, most notably DN excitations 
into charm-baryon resonances not too far from the DN threshold, e.g., A^(2593) 
and Sc(2625), as well as charm exchange into ttAc and ttSc channels. In Ref : - ^U ep- 
arable meson-baryon interactions have been employed with parameters constrained 
to dynamically generate the Ac (2593) state. Since the in- medium D-meson spec- 
tral function figures back into the T-matrix, one is facing a sclfconsistency problem 
(much like for the heavy- light quark T-matrix discussed in Sec. 12. 3. 2|) . Selfconsistent 
calculations including nuclcon Pauli blocking and dressing of intermediate pion and 
nucleon propagators result in D-meson spectral functions with a significant broad- 
ening of up to r^i ~ 100 MeV but a rather s mall shift of the peak position of about 
— 10 MeV (for qj^ = g ^ 0.16 fm"'^). In Ref! ^^^ [ a somewhat stronger coupling of 
DN to the Ac(2593) results in a stronger collective Z?A^^^ Ac(2593) mode (about 
250 MeV below the free Z?-meson mass) and a pertinent level repulsion which pushes 
up the "elementary" D-peak by ^30 MeV. Also in this calculation the broadening 
is significant, by about ^50 MeV. The D^ was found to be rather little affected, 
neither in mass nor in width. Investigations in the selfconsistent cou pled-c hannel 
framework have been extended to a nucleon gas at finite temperaturJI^ with a 
more complete treatment of DN scattering, cf. left and middle panels of Fig. [TBI 
The thermal motion of nucleons implies that a larger kinematic regime in the eenter- 
of-mass (cm) energies in the scattering amplitude is probed (compared to r=0). 
For the real parts this leads to a further averaging of the positive and negative 
parts of the amplitude, while the imaginary parts are negative definite (some loss 
of interaction strength may occur in channels with resonances close to threshold). 
More quantitatively, at normal nuclear matter density, the resulting mass shifts are 
10-20 MeV at T = (attractive for D and repulsive for D), decreasing to about 
half (or less) at T = 150 MeV. On the other hand, the Z?-meson width is around 
100 MeV at both zero and finite temperature, while the D width is small at T = 
but increases to about 30 MeV at T = 150 MeV. 

Medium modifications of D-mesons in a hot pion gas have been studied in 
Ref.™. 

The main idea in this work is to implement the recently discovered scalar 
1)0(2310) and axialvector D']^(2430) states as chiral partners of the pseudoscalar 
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Fig. 16. (Color online) Real and imaginary parts of D-meson selfenergies in dense and/or hot 
hadronic matter. Left and middle panels: selfconsistent coupled-channel calculations for in-medium 
D+ (left) and Z)~ (right) potentials based on DN scattering in a hot nucleon gas as a function 
of temperaturJSni for two different nuclear densities (and for two different inputs for the DN 
scattering length); the real (upper panels) and imaginary (lower panels) part of the potentials 
are defined in terms of the on-shell solfenergy via ?/((?) = S£)((.i;™(g), (j')/2(.i;™ with a;™ the 
quasiparticle energy; note that the width is given by F = — 2ImC/. Right panel: width (upper 
lines) and mass shift (lower lines) of D (top) andD*(2010) (bottom) mesons in a hot pion gas 
based on resonant scattering via chiral partneriSD note that the isospin symmetric pion gas 
implies equal effects on D-mesons (-D+, -D") and anti-D mesons (-D~, D''). 



D and vector £)*(2010) mesons, respectively. Their large widths of 200-400 MeV 
are primarily attributed to S'-wavc pion decays into D and D* . In a thermal pion 
gas, Do(2310) and £'']^(2430) therefore act as strong resonances in Dtt and D*7r 
scattering, which have been treated in Breit-Wigner approximation. In addition, D- 
wave resonances, 151(2420) and £'2(2460), have been accounted for. The resulting 
collisional widths of D and D* reach up to about ~40-60MeV for temperatures 
around T = 175 MeV, while the mass shifts are attractive up to —20 MeV. It 
can be expected that (e.g., in a selfconsistent calculation) the inclusion of medium 
effects on the resonances (e.g., chiral partners are expected to approach degeneracy 
towards chiral restoration) will lead to a reduction in the mass shift (not so much 
for the widths). Above the critical temperature, it is then natural to switch to a 
quark-based description, i.e., c~\- q scattering, much like in the effective resonance 
model discussed in Sec. 12.3.11 

When combining the effects of pion and nucleon scattering on D-mesons in hot 
hadronic matter, their total width at temperatures around Tc adds up to V^^i^ — 
180 MeV) ~ 150 MeV. This is only by about a factor of ~2 smaller than what was 
found for c-quarks at T = 1.1 in the T- matrix approach for c-q scattering, cf. right 
panel of Fig. [TTJ Since it can be expected that other excited hadrons contribute to 
i3-meson rescattering (albeit with less strength) , D-meson transport properties may 
not be much different from those of c-quarks in the QGP, at least at temperatures 
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close to Tc- It is therefore of considerable interest to employ in-medium _D-meson 
T-matrices to evaluate heavy-Savor transport coefficients in hadronic matter. 

3. Heavy-Quark Observables in Relativistic Heavy-Ion Collisions 

One of the main motivations for the vigorous theoretical studies of HQ diffusion in 
the QGP is the possibility of utilizing HQ observables in ultrarelativistic heavy-ion 
collisions as a quantitative probe of the matter produced in these reactions. If the 
latter reaches approximate local thermal equilibrium, such applications can be per- 
formed by solving the Fokker- Planck equation for a heavy quark diffusing within a 
collectively expanding background medium with space-time dependent temperature 
and flow field (applicable for "sufficiently slow" charm and bottom quarks). This is 
typically realized with a Monte-Carlo simulation using a test-particle ansatz for an 
equivalent stochastic Langevin equation. In such a formulation, a direct relation be- 
tween the input in terms of a (temperature-dependent) HQ diffusion coefficient and 
the modifications of HQ spectra in the evolution can be established. In addition, 
the Langevin formulation admits an efficient implementation of the dissipation- 
fluctuation relation for relativistic kinematics. 

Alternatively, the modiflcations of HQ spectra in URHICs have been evaluated 
by implementi ng test parti cles into numerical transport simulations of the back- 
ground 

niediunPlIlSlIMl. 

This, in principle, accounts for non-equilibrium effects 
in the medium evolution (which could be particularly relevant for high-pt particles 
in the bulk), but the connection to the diffusion concept becomes less direct (HQ 
cross sections need to be evaluated in an equilibrium medium to extract "equiva- 
lent" diffusion coefficients). However, when analyzing theoretical predictions of HQ 
spectra we also compare to results of transport models for the bulk evolution. 

We start our presentation in this Section by briefly outlining how the Fokker- 
Planck equation is implemented into numerical simulations based on a relativistic 
Langevin process (Sec. 13. This is followed by a discussion of different models 
for the background medium in relativistic heavy-ion collisions, where we focus on 
thermal models including hydrodynamics and expanding fireballs (Sec. 13. 2p . The 
main task here is to provide realistic benchmarks for the conversion of transport co- 
efficients into modifications of HQ spectra. This furthermore requires the definition 
of controlled initial conditions for the HQ spectra (Sec. 13. 3p . usually taken from 
p-p collisions, possibly augmented by nuclear effects (in particular a "Cronin" pt- 
broadening). Available Langevin simulations combining different inputs are quanti- 
tatively compared at the level of the final-state HQ spectra resulting from the QGP 
and "mixed" phases in central and scmiccntral Au-Au collisions at RHIC fSec. l3.4() . 
Even though HQ spectra arc not observable, they provide the cleanest theoretical 
level of comparison, before further processing through hadronization, electron de- 
cay and charm/bottom composition occurs. The latter three steps are necessary to 
enable comparisons to currently available electron data (Sec. 13. 5|) . and thus arrive 
at an empirical estimate of the HQ diffusion coefficient characterizing the QCD 
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medium produced at RHIC. In a more speculative step, the extracted HQ trans- 
port coefficient may be used to schematically estimate the ratio of shear viscosity 
to entropy density (Sec. 13.6"]) . which has recently received considerable attention in 
connection with viscous hydrodynamic simulations at RHIC. 



3.1. Relativistic Langevin Simulations 

The Fokker-Planck equation, introduced in Sec. 12.11 is equivalent to an ordinary 
stochastic differential equation. Neglecting mean-field effects of the medium, the 
force acting on the heavy particle is divided into a "deterministic" part, describing 
its average interactions with the light particles in the medium (friction or drag), 
and a "stochastic" part, taking into account fluctuations around the average on 
the level of the standard deviation. Thus the relativistic equations of motion for a 
heavy quark become a coupled set of stochastic differential equations, which for an 
isotropic medium can be written in the form 

dxj = ^dt , 

E (63) 

dpj = -Tpjdt + VdtCjkPk , 

where E = (mg +p^)^/^, and T and Cjk are functions of {t,x,p) with j,k=l,2,3; 
they are related to the transport coefficients A and B (discussed in the previous 
section) below. F and Cjk describe the deterministic friction (drag) force and the 
stochastic force in terms of independent Gaussian-normal distributed random vari- 
ables p = {pi,p2,P3), 

P(p)=(^)^xp(4) , (64) 

In the limit dt 0, the covariance of the fluctuating force is thus given by 

(^F^^\t)Fl:''\t')) = C,iCki6{t-t') . (65) 

However, with these specifications the stochastic process is not yet uniquely defined, 
but depends on the spe cific choice of the momentum argument of the covariance 
matrix, Cjk, in Eq. (IMlP^, i.e., the definition of the stochastic integral. Usual 
schemes arc given by the pre-point Ito , the mid-point Stratonovic-Fisk, and the 
post-point Ito (or Hanggi-Klimontoviclll^^ interpretation of the stochastic inte- 
gral. We can summarize all these realizations of the stochastic process by specifying 
the actual momentum argument in the covariance matrix by 

Cjk~^Cjk{t,x,p + ^dp) , (66) 

where ^ ~ 0,1/2,1 corresponds to the pre-point Ito, the mid-point Stratonovic, 
and the post-point Ito realizations, respectively. The equation for the corresponding 
phase-space distribution function can be found by calculating the average change of 



August 1, 2009 6:56 WSPC/INSTRUCTION FILE hq-qgp4 



Heavy quarks in the quark- gluon plasma 41 

an arbitrary phase-space function, g{x,p), with time. According to Eq. (|63p . with 
the specification Eq. (|66p of the stochastic process, we find 



(67) 



(,g(. + d.,p + dp) - ,g(.,p)) = (|^| + ^ {-Tp, + e^Q. 

, 1 0^9 



2 dpjdpk 



CjiCki 




Here, the arguments of both, F and Cjk, have to be taken at {t,x,p) since the 
corrections are of the neglected order, 0{dt^^^). In the derivation of this equation 
the statistical properties of the random variables pi, implied by Eq. (|64p . 



(p,)=0, {p,Pk)=S,k, (68) 



have been used. It follows that the average of an arbitrary phase-space function is 
by definition given by the phase-space distribution function for the heavy particle 
(in our context a heavy quark), fQ{t,x,p), e.g., 

d_ 

dt 

After integrations by parts, and since Eq. ((^7|) holds for any function g, one finally 
arrives at the Fokker-Planck equation. 



iaix^p)) ^ J d^x J d^pg{x,p)^fQ{t,x,p) 



(69) 



dfg _^ Pldfq ^ _d_ 
dt E dxj dpj 



The drag term, i.e., the first term on the right-hand side of this equation, depends 
on the definition of the stochastic integral in terms of the parameter ^. Comparison 
with Eq. (jlOp shows that, independent of the choice of ^, the covariancc matrix is 
related to the diffusion matrix by 

Qfc = V2B'oP^i + V^iPI , (71) 

while the friction force is given by 

Tp, ^ Ap, ^ ^Qk^ ■ (72) 
dpi 

Numerical investigations have shown that the drag and diffusion coefficients in- 
ferred from microscopic models according to Eqs. p5|) in general do not warrant 
a good agreement of the long-time limit of the solution to the Fokker-Planck evo- 
lution with the relativistic equilibrium Jiittner-Boltzmann distribution (where the 
temperature is given by the background medium) . The problem is with the longitu- 
dinal diffusion coefficient, Bi, which induces an overestimate of the corresponding 
fluctuating forces. Thus, one typically adjusts the drag coe fficient by c hoosing Bi 
in Eq. ([70|) to satisfy the asymptotic equilibration conditiorff^IIISHEH] 

To find the dissipation-fluctuation relation, imposed by the equilibration condi- 
tion, we first study the heavy quark's motion in a heat bath in thermal equilibrium 
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in its rest frame. Then the momentum distribution of the heavy quarks should 
become a Jiittner-Boltzmann distribution, 

rQ'^{p)^cxp{-E/T) (73) 

with the temperature, T, imposed by the heat bath. For a Langevin process with 
Bo = Bi= D, i.e., 

Cjk = ^2D{E)5jk , (74) 

where the diffusion coefficient has been written as a function of the heavy quark's 
energy, _E, the cquihbration condition for a given parameter ^ in Eq. (|66|) is obtained 
by using Eqs. ^ and Eq. ^ in Eq. dTO]): 

A{E)ET-D{E)+T{l-S)D'{E)^Q . (75) 

Since the drag and diffusion coefficients are usually given numerically, the most 
convenient update rule for the Langevin process is achieved by setting ^ = 1, i.e., 
using the post-point Ito (Hanggi-Klimontovich) rule for the stochastic integral in 
Eq. (|63[) and imposing the simple relativistic dissipation-fluctuation relation, 

D ^ AET . (76) 

This guarantees the proper approach of the heavy quark's phase-space distribution 
to the appropriate equilibrium distribution with the temperature imposed by the 
heat bath. 

For the more general form of the cov ariancc matrix, Eq. (|7ip , the post-point Ito 
value, C = 1, has been chosen in Rcf.l^, and the longitudinal diffusion coefficient is 
set to 

Bi = AET , (77) 

while the drag coefficient A as well as the transverse diffusion coefficient, Bp, are 
used as given by Eq. (|15[) for the various microscopic models for HQ scattering in 
the QGP. Com pari ng to Eq. ((75l) . one finds that this is equivalent to the strategy 
followed in ^ei¥^ of using the prepoint-Ito rule, ^ = 0, but to adjust the drag 
coefficient according to the dissipation-fluctuation relation Eq. ((75|) . 



3.2. Background Medium in Heavy-Ion Collisions 

For HQ transport coefficients computed in an equilibrium QGP, the natural and 
consistent framework to describe the evolving medium in heavy-ion collisions are hy- 
drodynamic simulations, formulated in the same (thermodynamic) variables. This 
choice is further rendered attractive by the success of ideal hydrodynamics in de- 
scribing bulk observables at RHIC, i n particular pt spectra and elliptic flow of 
the most abundant species of hadronJI^E^IISIIS] r^Yie agreement with meson and 
baryon spectra typically extends to pr — 2-3 GeV, respectively. At the parton level, 
this converts into a momentum of pf ~ 1 GeV, which approximately coincides with 
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the "levcling-off" of the experimcntaUy observed V2{pt) (at higher momenta hy- 
drodynamics overestimates the elhptic flow). On the one hand, this appears as a 
rather smah momentum in view of the ambition of describing HQ spectra out to, 
say, ~ 5 GeV. However, one should realize that (a) more than 90% of the bulk 
matter is comprised of light partons with momenta below pt ~ 1 GeV, and (b) 
the velocity of a pf = 5 GeV charm quark (with rric ~ 1.5 GeV) is very similar 
to a pt = 1 GeV light quark (with rriq = 0.3 GeV). This suggests that most of 
the interactions of a pt = 5 GeV charm quark actually occur with soft light par- 
tons (which are wel l des cribed by a hydrodynamic bulk). This has been verified by 
explicit calculation^!^ and is, after all, a prerequisite for the applicability of the 
Fokker-Planck approach (i.e., small momentum transfer per collision). On the other 
hand, one may be concerned that the overestimate of the experimentally observed 
elliptic flow at intermediate and high pt within hydrodynamics may exaggerate the 
HQ elliptic flow in Langevin simulations. This is, however, not necessarily the case, 
since the transfer of V2 from the bulk to the heavy quark critically depends on the 
light-parton phase space density (the drag coefficient is proportional to it); since 
the hydrodynamic spectra fall significantly below the experimental ones at higher 
Pt, the phase space density of the hydrodynamic component is relatively small. It 
is therefore not clear whether the (small) fraction of thcrmalized particles at high 
Pt implies an overestimate of the total V2', this may be judged more quantitatively 
by comparing to transport calculations. 

In parallel to hydrodynamic descriptions of the bulk medium, expanding fireball 
models have been employed. These are simplified (and thus convenient) parameter- 
izations of a full hydrodynamic calculation in terms of an expanding volume and 
spatial flow-velocity field, but otherwise based on similar principles and variables. 
E.g., entropy conservation is used to convert a given volume into a temperature 
via an underlying equation of state (EoS). The reliability of a fireball model cru- 
cially hinges on a realistic choice of the parameters, in connection with a proper 
description of the final state in terms of particle production and collective flow. In 
principle, a fireball model offers some additional flexibility in varying the evolu- 
tion, which may provide useful checks of the sensitivity to specific aspects of the 
expansion. 

Several key parameters of thermal medium evolution models employed in HQ 
Langevin simulations are compiled in Tab.[T] The starting point of both hydro and 
fireball models are the initial conditions of the thermal medium, characterized by 
a formation time when the medium is first assumed to be (locally) equilibrated. At 
RHIC, typical formation times are estimated to be in the range of tq ~ 0.3-1 fm/c. 
With a total entropy fixed to reproduce the measured rapidity density of hadrons 
at a given centrality, e.g. at impact parameter 6 ~ 7 fm/c, these formation times 
translate into average initial temperatures of Tq ~ 250-350 MeV. If the entropy 
density scales as s (xT^, one can roughly compare the initial conditions in different 
approaches using S = sqVo and Vq oc tq. E.g., an initial T™^'' = 265 MeV based 
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M 1 




Arlxl 




To [fm/c] 


1.0 


0.33 


0.6 


0.6 


To [McV] 


265 (max) 


340 (avg) 


250 (avg) 


330 (max) 
260 (avg) 


Te [MeV] 


165 


180 


170 


165 


initial 


wounded 


isotropic 


lin. comb, of 


lin. comb, of 


spatial 


nuclcon 




iVcoii and iVpart 


TVcoii and A^part 


b [fm] 


6.5 


7 


5.5 


7 


bulk-i'2 


5 


5.5% 


3% 


4.75% 


TpB [fm/c] 


-9 


-5 


-9 


-7 


QGP-EoS 


massless 


massless 


massless 


massless 




{Nf = 3) 


(Nf = 2.5) 


(Nf = 3) 


{Nf = 3) 


HQ Int. 


pQCD HTL 


pQCD+rcso 
pQCD+T-mat 


AdS/CFT 


pQCD run. as 


TOc,6 [GeV] 


1.4 


1.5, 4.5 


1.5, 4.8 


1.5, 5.1 



Table 1. Survey of parameters figuring into hydrodynamic and fireball evolutions employed in the 
Langevin simulati ons of HQ spectr a for s emicentr al Au-Au collisi ons at RHIC, corresponding to 
Refs. [MTpiftol [HGRpMEU [AHHfllSEI] and [GAjSZfm 

on To = Ifm/c increases by a factor of 3^/"^ upon decreasing To = 0.33 fm/c, re- 
sulting in T^'^^ ~ 382 MeV; if the number of light flavors in the EoS is reduced. 
To increases as well; e.g., Tq^* = 260 McV based on To = 0.6fm/c and Nf~3 (as 
in Ref. increases to Ti^"^ = 260 MeV (0.6/0.33)^/3 (47 5/42,25)^3 ^ 330 MeV 
for To = 0.33 fm/c and Nf=2.5 (reasonably consistent with Ref. cf. Tab. [1]). 
The QGP-dominatcd evolution lasts for about 2-4 fm/c, followed by a mixed phase 
of similar duration at a critical temperature Tc ~ 165-180 MeV. The effects of a 
continuous (cross-over) transition, as well as of the hadronic phase, have received 
little attention thus far, but are not expected to leave a large imprint on HQ observ- 
ables. After all, the cross-over transition found in IQCD exhibits a marked change 
in energy density over a rather narrow temperature interval. A more important 
aspect is the consistency between the EoS used to extract the temperature of the 
bulk evolution and the corresponding degrees of freedom figuring into the calcu- 
la tion of the HQ transport coefficients. In hydrodynamical backgrounds used thus 
faJ2SIllSEll the evolution is described with a 2-|-l-dimensional boost-invariant sim- 
ulation with an ideal massless-gas EoS. The initial state is typically determined by 
distributing the entropy in the transverse plane according to the wounded nuclcon 
density. Unfortunately, the impact parameters in current Langevin simulations vary 
somewhat, which is particularly critical for the magnitude of the subsequently de- 
veloped elliptic fiow. The value of the critical temperature has some influence on 
the QGP lifetime (lower temperatures leading to larger QGP duration), as does 
the hadron-gas EoS (more hadronic states imply a larger entropy density at Tc and 
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thus a reduced duration of the mixed phase). The termination point of the evolu- 
tion (beginning, middle or end of mixed phase) is rather significant, especially for 
the HQ V2 which needs about 5fm/c to build up most of its magnitude. 

In Refs .122145] ^j^g medium is parameterized, guided by the detailed hydrody- 
namic calculations of 

Ref.I142], 

as a homogeneous thermal elliptic "fire cylinder" of 
volume V{t). The QGP temperature is determined via the QGP entropy density, 
s, under the assumption of isentropic expansion (total 5* = const), 

^ = 4 = |^T3(16 + 10.5iV,). (78) 

The thus inferred temperature is used in Eq. (|15|) to compute the friction coeffi- 
cients, A, and transverse diffusion coefficient, Bqi with the longitudinal diffusion 
coefRcient fixed by the dissipation- fluctuation relation, Eq. (j77p . In the mixed phase 
at Tc = 180 MeV the QGP drag and diffusion coefficients are scaled by a factor 
oc Q^/'^ to account for the reduction in parton densities (rather than using hadronic 
calculations). Special care has to be taken in the parameterization of the elliptic 
flow in noncentral Au-Au collisions: the contours of constant flow velocity are taken 
as confocal ellipses in the transverse plane with the pertinent transverse flow set 
consistently in perpendicular direction. The time evolution of the surface velocity of 
the semi-axes of the elliptic fire cylin der p arameterizes the corresponding results of 
the hydrodynamic calculations in in particular the time-dependence of the 

elliptic-flow parameter, for the light quarks. The parameters are adjusted such 

is) 

that the average surface velocity reaches v\_ = 0.5c and the anisotropy parameter 
V2 = 5.5% at the end of the mixed phase. Finally, the velocity field is specified by 
scaling the boundary velocity linearly with distance fro m th e center of the fireball, 
again in accordance with the hydrodynamic calculatiorP^H. 

3.3. Initial Conditions and Hadronization 

The Langevin simulations of HQ diffusion in the QGP require initial conditions 
for charm- and bottom-quark phase-space distributions. For the spatial part of the 
initial distribution in the transverse plane all calculations adopt binary-collision 
scaling following from a Glauber model, reflecting a hard process for the primordial 
production mechanism. Furthermore, all calculations thus far focus on a limited 
rapidity window around midrapidity, where the longitudinal distributio n is a ssumed 
to be uniform in space-time rapidity. As for the initial HQ pt spectra, Ref.^ employs 
a fit to a leading-order parton-model calculation from the CompHEP package^^, 

d^d^.^^^'^-y^WTA^' (79) 

with a = 3.5 and A = 1.849 GeV. In Refs.lSHHSl, PYTHIA resuhs for c-quark 
spectra have been tuned to reproduce available D-meson spectra in d-Au collisions 
at RHIC (assuming 5- function fragmentation, cf. left panel of Fig.fTT]). The pertinent 
semileptonic single-electron decay spectra approximately account for p-p and d-Au 
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PjlGeV\ PjlGeY] 



Fig. 17. (Color online) Fits of D- and D*-mcson px spectra in 200j4GcV d-Au collisions at RHIC 
with a modified PYTHIA simulation (left pane l) and the corresponding non-photonic single- 
electron pt spectra in p-p and d-Au coUisionsi^^ The missing yield of high-py electrons is fitted 
with the analogous B-meson decay spectra, thus fixing the bottom-charm ratio at Cj^g/crcc — 
4.9 • 10-3. 



spectra up to pt = 4 GcV; the missing part at higher pT is then supplemented by 
B-mcson contributions. This procedure results in a crossing of the D- and B-meson 
decay electrons a.t pr — 5 GeV and a cross-section ratio of (Ji,i,/(Jcc — 4 .9 ; l O""^ (see 
right panel of Fig. [T7|) . which is within the range of pQCD prediction j^^^ i 

With initial conditions and bulk medium evolution in place, one can evolve HQ 
phase-space distributions through the QGP (and mixed phase) of a heavy-ion colli- 
sion. The final HQ spectra, however, require further processing before comparisons 
to observables can be made. First, one has to address the hadronization of the HQ 
spectra into charm and bottom hadrons {D, D* , Ac etc.). Two basic mechanisms 
have been widely considered in hadronic collisions, i.e., fragmentation of an individ- 
ual quark and recombination with an extra quark from the environment. In general, 
the former is mostly applicable for high-energy partons while the latter requires a 
sufficient overlap of the mesonic wave function with the phase-space density of 
surrounding quarks and is therefore more relevant toward lower momentum. 

The fragmentation of a quark is implemented by applying the factorization 
theorem of At large transverse momenta, the production process of a par- 

ton occurs on a short time scale, Tprod — l/pt, while hadronization occurs at the 
considerably larger time scale Thad — 1/Aqcd- Thus the production cross section 
for a hadron can be factorized into an elementary parton-production cross section 
(hard process) and a phenomenological universal transition probability distribution, 
Dh/i{z), for a parton i of momentum pi to convert into a hadron with momentum 
fraction z = Phi Pi < 1- For light quarks and gluons the fragmentation functions, 
Dh/i, are rather broad distributions around z ~ 0.5, but for heavy quarks they 
become rather sharply peaked toward z = 1 and are sometimes even approximated 
by a i5-function, D[z) = 5{1 — z). 

The mechanism of recombination of a produced quark with other quarks or an- 
tiquarks in its environment (e.g., the valence quarks of the colliding hadrons) has 
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first been introduced in tlie late 1970's to explain flavor asymmetries in tt and K 
meson production in hadronic collisions at forward rapiditie j'^^^ i In particular, the 
recombination i dea has b een rather successful in describing flavor asymmetries in 



lisions, quark coalescence models, applied at the hadronization transition, provide a 
simple and intuitive explanation for the obse rved con stituent-quark number scaling 
(CQNS) of the elliptic flow of light hadronsil^HilSOl 

and the (unexpectedly) large 
baryon-to-meson ratios (e.g., p/n ~ 1 or {K + K) / {AK^) ~ 1.3 in central 200AGeV 
Au-A u collisions at RH IC) at intermediate transverse momenta (2 GeV ^ Pt S 
5 GeVj ^^^ l ^^^ l -'^^^ l ^^'^ l. CQNS refers to a scaling property of the hadronic elliptic 
flow, V2,h{PT), in terms of a universal function V2^q{pT/n) = V2,h{PT) /n, where n 
denotes the number of constituent quarks in a given hadron. h. CQNS naturally 
emerges from the recombination of approximately comoving quarks and antiquarks 
in a collectively flowing medium. Thus, within this picture, V2,q{pt) is interpreted 
as a universal elliptic flow of the quarks with transverse momentum pt at the mo- 
ment of hadronization (typically assumed to be the quark- hadron transition at Tc). 
It can be expected that the phenomenologically very successfu l coalescence con- 
cept also applies in the HQ sector of heavy-ion 

collisionfPUM]. 

Note that, unlike 

quark fragmentation, quark recombination adds momentum and elliptic flow to the 
produced hadron (through the quark picked up from the environment). 

At this point it might be instructive to reiterate a conceptual connection be- 
tween the quark coalescence model and the idea of resonance correlations in the 
QGP. The latter were found to be an efficient mechanism for arriving at a small 
HQ diffusion constant, both within the effective resonance model (Sec. 12.3. ip and 
within the T-matrix approach fSec. l2.3T^ . Especially in the T-matrix approach, the 
resonance correlations were found to strengthen toward the expected hadronization 
transition, and thus provide a natural emergence of heavy-light quark coalescence 
at Tc- These ideas have recently been implemented in a resonance-based descrip- 



tion of the coalescence process in kinetic theory based on a Boltzmann equatio n ' I 
This approach improves instantaneous coalescence formulations in that it respects 
energy conservation and establishes a well-defined equilibrium limit in the coales- 
cence process (i.e., the thermal distribution for the formed meson). Subsequently, 
resonance-recombination has been combined with "realistic" quark phase-space dis- 
tributions, as generated in relativistic Langevin simulation j-^^^ i In particular, it 
was found that CQNS could be recovered under the inclusion of space-momentum 
correlation s in th e quark phase-space distributions. 

In Refs.^2HIlll HQ spectra at RHIC have been hadronized in a combined coales- 
cence plus fragmentation scheme. For the hadronization of, e.g., charm quarks into 
D mesons one obtains the _D-meson spectra as 



the charm sector I I, even close to midrapidity. In the context of heavy-ion col- 






dy d'^pT dy d'^pT dy d'^pr 



For the first term on the right-hand side, the quark-coalescence model of Rem 
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has been employed, where the px spectrum of a Z? meson fohows from a convolution 
of light anti-quark and charm-quark phase-space distributions, fq,c, as 

J Ty^coal f n ■ r]rr C 

^ ^'""J hrJ d^9/^('7>-)/.(P,,-.)/c(Pc,r-.) . (81) 

Here, p = Pq + Pc the £)-meson momentum and go a combinatorial factor 
accounting for color-neutrality and spin averaging. The D-meson Wigner function, 
/£)(g, x), is assumed as a double Gaussian in relative momentum p^ — p^ and size, 
Tc — Tq, and da is the hyper-surface element 4-vector of the hadronization volume. 
The charm-quark distribution corresponds to the Langevin output at the end of 
the mixed phase of the fireball model, while the light-quark distributions are taken 
from pr eviou s applications of the coalescence model to light-hadron observables 
at This represents a parameter-free conversion of HQ distributions into 

heavy- meson spectra (note that the final state of the expanding fireball modePSIhas 
been matched to the parameterization of collecti ve ve locity and elliptic flow for the 
light-quark distributions in the coalescence modeff^H) . The coalescence mechanism 
does not exhaust all heavy quarks in the hadronization process, especially toward 
higher pt (where the light-quark phase-space density becomes small). Therefore, the 
remaining heavy quarks are hadronized using fragmentation, which for simplicity 
is treated in (5-function approximation (as has been done in connection with the 
initial conditions). The formation of baryons containing heavy quarks (e.g., Ac) 
has been neglected since it has been found to give only small contributions, i.e., 
Ac/D <^ 1. Quant itative refinements should, however, include these processes, see, 
e.g., Refs.™6Q]. 

Finally, the comparison to electron spectra requires to compute semileptonic de- 
cays of heavy-flavor hadrons. Thus far, these have been approximated in three-body 
kinematics, e.g., D evK . An important finding in this context is that t he result - 
ing electron V2(j)t) traces the one of the parent meson rather accurateljff^^^^^, 
implying that electron spectra essentially carry the full information on the heavy- 
meson V2 ■ In the Pt spectra, the decay electrons appear a t roughly half of the mo- 
mentum of the parent meson. It has also been pointed oulf^^^^^^^that Ac baryons 
have a significantly smaller branching fraction into electrons (about 4-5%) than D 
mesons (7% and 17% for neutral and charged D's, respectively). Thus, in case of a 
large Ac/D enhancement, a net electron "loss" could mimic a stronger suppression 
than actually present at the HQ level. In f act, even variations in the neutral to 
charged chemistry from p-p to A- A collisionJl^ could be quantitatively relevant. 



3.4. Model Comparisons of Heavy- Quark Spectra at RHIC 

We are now in position to conduct quantitative comparisons of diffusion calculations 
using transport simulations for HQ spectra in 200 AGcV Au-Au collisions at RHIC. 
We focus on Langevin simulations but also allude to Boltzmann transport models. 
The modifications of the initial spectra are routinely quantified in terms of the 
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nuclear modification factor, Raa, and elliptic-flow parameter, V2-i defined by 



RAAipub) = 

V2{pt;b) = 



.,,div^ (82) 
/d'/'d^;f|^cos(20) 

respectively; dNQ^{b)/dpt denotes the HQ pt spectrum in an A- A collision at im- 
pact parameter, b, which is scaled by the spectrum dNj^^ /dpt from p-p collisions 
times the number of binary nucleon-nucleon collisions, Nco\\(b) (to account for the 
same number of heavy quarks). Thus, any deviation of Raa from one indicates 
nuclear effects (from the QGP but possibly also in the nuclear initial conditions or 
from the pre-equilibrium stages). The elliptic-flow parameter, V2{pt), is the second 
Fourier coefficient in the expansion of the (final) momentum distributions in the 
azimuthal angle, (f), relative to the reaction plane {x-z plane) of the nuclear colli- 
sion. At midrapidity. where the "directed" flow (vi) is expected to vanish, the V2 
coefficient is the leading source of azimuthal asymmetries. A non-zero V2 is only 
expected to occur in noncentral A- A collisions due to an "almond" -shaped nuclear 
overlap zone (with a long (short) axis in y (x) direction). Typical sources for a non- 
zero elliptic fiow arc a path-length difference for absorption of particles traversing 
the reaction zone or an asymmetry in the collective (hydrodynamic) flow due to 
stronger pressure gradients across the short axis. Both effects convert the initial spa- 
tial anisotropy, ^2, in a positive momentum anisotropy in the particle pt-spectra. 
While the former mechanism is usually associated with high-pt particles (typical 
leading to a W2 < 5%), the latter is driven by collective expansion due to thermal 
pressure mostly applicable to low-pt particles (with significantly larger V2 values, in 
excess of 5%). Since in the Langevin simulations heavy quarks are assumed to be 
exclusively produced in primordial N-N collisions (i.e., their number is conserved 
subsequently), the HQ Raa can be simply calculated as the ratio of the HQ pt dis- 
tribution function at the moment of hadronization to the initial distribution (taken 
from p-p collisions) , 

n ./Q(ihad,Pt) 

Raa = —r-r, ^ i 

fQ(to,Pt) ^g^^ 

^ N /d0/Q(thad,Pt,0)cos(2(?!)) 

V2{pt) = -TY. N , 

while the V2 is computed using its definition given above. 

The next five figures (|18ll22p encompass calculations of Raa and V2 in semicen- 
tral Au-Au collisions at RHIC for the following approaches: 

(i) Fig. [TH] [MTpSl displays Langevin simulations for c quarks (with the prc-point 
Ito realization of the stochastic integral) using a hydrodynamic evolution for 
6 = 6.5 fm; the HQ drag and diffusion coefficients are based on LO hard-thermal 
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Fig. 18. (Color online) Nuclear modification factor (left panel) and elliptic flow (right panel) of 
charm quarks as a function of transverse momentum in semic ent ral (6 = 6.5 fm) Au-Au collisions 
using a hydrodynamic evolution of the bulk medium at RHI(j2S]. The calculations are performed 
for HTL-improved LO-pQCD scattering with variable strong coupling and fixed Debye-screening 
mass in t-channel gluon-exchange scattering (/!£> = 1.5 T). The relation of the spatial diffusion 
coefficient, Ds (denoted D in the figure legend), to the strong coupling constant, Os, is given by 
the approximate relation 2-KTDa ~ 6(0.5/os)^- 



loop scattering matrix elements with variable a,, but fixed Debye screening 
mass. 

(ii) Fig. [H [HGRpSI displays Langcvin simulations for c and b quarks (with the 
post-point Ito (Hanggi-Klimontovich) realization) using a thermal fireball ex- 
pansion for 6 = 7 fm; the HQ drag and diffusion coefficients arc based on 
the effective resonance-|-pQCD modcI-^^J for variable resonance width (coupling 
strength) an d a., = 0.4 in the pQCD part. 

(iii) Fig. [20] [AHHpJ^ displays Langcvin simulations for c quarks (with the pre-point 
Ito realization) using a hydrodynamic expansion for b = 5.5 fm; the HQ drag 
and diffusion coefficients are based on the strong-coupling limit with AdS / CFT 
corr espon dence with a variable coupling strength estimated from matching to 



QCEffoa. 

(iv) Fig. ED [hmgr|3SI 

displays Langcvin simulations as under (ii) but with HQ 
transport coefficients based on the T-matrix+pQCD approach for two IQCD- 
based input poten tials. 

(v) Fig. [22] [Mol] I displays Boltzmann transport simulations using a covariant 
transport model for b — S fm; the HQ interactions are modeled by schematic 
LO pQCD cross sections, including upscaling by "if factors" . 

Before going into details, let us try to extract generic features of the calculations. 
In all cases there is a definite correlation between a reduction in RAA{pt > 3GeV) 
and an increase in V2{pt), i.e., both features are coupled to an increase in interaction 
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Fig. 19. (Color online) The HQ Raa (left panel) and V2 (right panel) for scmicentral {b = 7fm) 
Au-Au collisions at RHIC within the effective resonance + pQCD model compared to results 
from LO pQCD elastic scattering only with as = 0.4 and corresponding Debye-screening mass 
= v'4-irOsT'. 

strength (decrease in the spatial HQ diffusion coeffieient). Furthermore, the V2{pt) 
shows a typical, essentially linear, increase reminiscent of a quasi-thermal regime 
followed by a saturation characteristic for the transition to a kinetic regime. In 
all Langevin calculations the saturation for charm quarks occurs at about pt = 2- 
3 GeV. For the largest interaction strength considered [Ds — l/(27rT)), the left 
panels of Figs. [18] and [20] even suggest a turnover of V2 (at this point one should 
recall that all calculations displayed in this section utilize elastic scattering only 
which is expected to receive appreciable corrections at high pt due to radiative 
processes). On the other hand, for pt — 2-3 GeV the nuclear modification factor 
is still significantly falling, leveling off only at larger pt ~ 5-6 GeV. As expected, 
bottom quarks exhibit much reduced effects for comparable diffusion constants due 
to their factor ^3 larger mass (see Fig. [T9l and lower panels in Fig. [2T|) . 

Next, we attempt more quantitative comparisons. Some re prese ntative numbers 
for the resulting Raa and V2 values are compiled in Tab. [ jf First we com- 
pare the LO-pQCD calculations for HQ diffusion in the hydrodynamic and fireball 
backgrounds corresponding to Figs. [T8l and [T9l respectively; for a comparable spa- 
tial diffusion coefficient, Ds — 24-30/ (27rr), both calculations show a maximal V2 of 
about 2% and a Raa {pt=5 GeV) ~ 0.7 (recaU the smaller 6 = 6.5 fm in [MT] vs. 7 fm 
in [HGR] which may lead to somewhat smaller V2, and the smaller Tq = 265 MeV 
[MT] vs. 340 MeV in [HGR] which entails somewhat less suppression). For the 
[AHH] hydro calculation with an AdS/CFT- motivated ansatz for the HQ friction 
constant, 



a diffusion constant of Dg = 21/( 2ttT) leads to similar results (note that Tab. [2] 
contains results for 6 = 7.1 while Fig. [20] is calculated for b ~ 5.5 fm). Let 

us now turn to stronger coupling, still focusing on the three Langevin simulations 




(84) 
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Fig. 20. (Color online) Raa (left) and V2 (right) of charm quarks resulting from hydrodynamic 
simulations of b = 5.5 fm Au-Au (y'sivjv = 200 GeV) collisions using AdS/CFT- motivated charm- 
quark diffusion constants with variable strength parameter, ""^HSl which corresponds to the con- 
stant C in Eq. | [84t . 

in Figs. \TE[ [TOl and which all utilize friction coefficients with a similar temper- 
ature dependence, essentially 7 oc (recall right panel of Fig. [75]) . correspond- 
ing to an approximately constant spatial diffusion constant times temperature, 
Ds{27tT) ~ const (recall left panel of Fig. [H]). For = G/{2ttT), all calcula- 
tions are again in semi-/quantitative agreement, w ith a maximum V2 of 4-6% and 
RAA{pt=5GeV) ~ 0.25-0.3. The 4% value for Ref.ED [AHH] will increase some- 
what if the hydro evolution is run to the end of the mixed phase rather than 
terminated in the middle of the mixed phase (this is supp orte d by the discussion 
in Sec. 13. 5. We also note that in the fireball model of 

Ref.|2Hl 

[HGR] the inclusive 



Model [Ref.] 


Ds{2TrT) 


b [fm] 




RAA{pt=5 GeV) 


hydro + LO-pQCEi^ 


24 


6.5 


1.5% 


0.7 


hydro + L0-pQCeI2S1 


6 


6.5 


5% 


0.25 


fireball + L0-pQCeI2H1 


-30 


7 


2% 


0.65 


fireball + reso+LO-pQCrl^ 


-6 


7 


6% 


0.3 


hydro + "AdS/CFT" (|84|'-^^^ 


21 


7.1 


1.5-2% 


-0.7 


hydro + "AdS/CFT" 1^^^ 


2tt 


7.1 


4% 


-0.3 


transport + LO pQCEtii^ 


-30 


8 


-2% 


-0.65 


transport + LO pQCeM^ 


-7 


8 


10% 


-0.4 



Table 2. Overview of model approaches (1^' column) and input parameters (2"^^ column: spatial 
charm-quark diffusion coefficient, S''^ column: nuclear impact parameter) for Langcvin simulations 
of charm-quark spectra in Au-Au collisions at RHIC; selected values for the resulting elliptic flow 
^j^max ^ V2{pt = 5 GeV)) and nuclear modification factor are quoted in columns 4 and 5. The 
last two rows represent charm-quark transport calculations in a transport model for the bulk. 
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Fig. 21. (Color online) Charm- (top row) and bottom-quark (bottom row) spectra in Au- 
A\i(^snn = 200 GeV) collisions at RHIC using Langevin simulations for HQ diffusion in an ex- 
panding fireball model. The results using HQ difl^usion based on LO-pQCD and reso nance-l-pQCD 
approaches (see also Fig. I19I I arc compared to the T-matrix+pQCD calculationJS] (the pertinent 
uncertainty band reflect s differ ent input potentials based on parameterizations of IQCD HQ free 
energies as given m Refs.Enm, The left (right) panels show the Raa ("2) for central (semiccntral) 
collisions. 



(pt-integratcd) V2 at the end of the mixed phase was adjusted to the experimen- 
tally observed light-particle 112 — 5.5-6 % at an impact parameter of & = 7fm, i.e., it 
presumably includes an extra 20% of bulk-W2 compared to the hydrodynamic calcu- 
lation^. Such an amount is typically built up in the subsequent hadronic phase of 
hydrodynamic evolutions and thus not present in pertinent HQ simulations within 
a QGP (-1- mixed) phase. 

We now make some comments specific to individual calculations. The Langevin 
calculations using HQ T-matrix interactions (supplemented with pQCD scattering 
off gluons) shown in Fig. [21] are rather close to the effective resonance model, even 
though they do not involve tunable parameters. However, they are still beset with 
substantial uncertainty, as indicated by the use of two different input potentials 
(in addition, the use of the free energy, Fi, instead of the internal energy, Ui, as 
potential significantly reduces the effects). One also notices that the V2 at low pt 
is very similar to the resonance model while the suppression at high pt is some- 
what less pronounced. This is so since the T-matrix transport coefficients (a) fall 

'This adjustment ensures compatibility of the fireball freezeout with the coalescence modeffSU, 
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off stronger with three-momentum (the resonant correlations arc close to the Q-q 
threshold) , and (b) decrease with increasing temperature (resonance melting) . The 
latter combines with the facts that the suppression is primarily built up in the very 
early stages (where the T-matrix is less strong) while the bulk V2 takes a few fm/c 
to build up (at which point the T-matrix has become stronger). Furthermore, the 
T-matrix calculations lead to stronger medium effects on b quarks than the effec- 
tive resonance model; this reflects the stronger binding due to the mass effect in 
the T-matrix calculation. 

A principal limitation of the Langevin approach is the treatment of fluctua- 
tions which are by definition implemented in Gaussian approximation. The latter 
arises due to enforcing the dissipation-fluctuation relation (mandatory to ensure 
the HQ distributions to approach equilibrium) which tends to underestimate the 
momentum fluctuations especially at high momentum, compared to a full transport 
calculation. This leads to an overestimate of the quenching effect at high pt even for 
the same average energy loss. One may assess these limitations mor e quanti tatively 
by comparing to Boltzmann simulations including partonic phaseJ^^^EH^ an ex- 
ample of which is displayed in Fig.[22]for charm quarks in 6 = 8fm Au-Au collisions 
at RHIC. The baseline LO-pQCD calculations indicated by the crosses in Fig. [22l 
labeled by "1.33 mb", may be compared to the fireball-Langevin simulations rep- 
resented by the blue lines in Fig. \T9\ In both calculations the underlying elastic 
parton-HQ cross sections correspond to a strong couphng constant of as — 0.4. 
The quenching and elliptic flow come out quite similar in both calculations at least 
up to Pt ~ 5 GeV, especially when accounting for the slightly different centrality. 
E.g., in the Boltzmann treatment, the Raa for Pt = ^ GeV charm quarks is about 
0.6-0.7 with a W2 of a few percent. For a four-fold increase of the cross section (which 
would roughly correspond to a reduction of Ds{2ttT) from ~30 to ^7), one finds 
RAA{pt = 5 GeV) ~ 0.4 and a maximum V2 of close to 10%. While the latter value 
is somewhat larger than the Langevin predictions, the agreement is not too bad. 

Finally, Fig. [53] shows results from an exploratory calculation in the Langevin 
approach where HQ drag and diffusion coefficients from elastic scattering in the 
effective resonance model (cf. Sec. I2.3.ip are conibined with induced gluon radia- 
tion in the DGLV E-loss formalism (cf. Sec. I2.6[j ^^'^ [ One of the uncertainties in 
this calculation is the extrapolation of the radiative E-loss into the low-momentum 
regime, where it still contributes rather substantially; e.g., the elliptic flow of charm 
quarks is increased over elastic pQCD+resonance model by ca. 40%, and even more 
(ca. 100%) relative to pQCD elastic scattering only. Another limitation is the above 
mentioned caveat in Langevin theory of underestimating the (E-loss) fluctuations 
implying an overestimate of the quenching at high pt. This c an also be seen when 
comparing to the pQCD radiative E-loss calculationEpIIS^EH]^ where the gluon radi- 
ation is treated microscopically within an opacity expansion. A consistent merging 
of radiative and elastic processes in HQ transport thus remains a challenging task. 
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Fig. 22. (Color online) Raa (left) and V2 (right) of charm qua rks re sulting from covariant transport 
simulations of fe = 8 fm Au-Au (y^sjvjv = 200 GeV) coUisionilMl 
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Fig. 23. (Color online) The nuclear modificat ion factor, Raa (left panel ) . and ellipt ic flow, t>2 
(right panel), for heavy quarks from coUisional^Sl and radiative E-losJEU^ cf. Ref.CSSI 



3.5. Heavy-Meson and Electron Observables 

To compare to observables, the HQ spectra discussed in the preceding section need 
to be converted into spectra of color-neutral final-state particles. At the minimal 
level, this requires hadronization into charm and bottom mesons and baryons. Thus, 
a measurement of identified HQ hadrons constitutes the most direct way to make 
contact with theoretical predictions. Currently, the richest source of information on 
HQ spectra in Au-Au collisions at RHIC are single-electron (e*) spectra, which, 
after the subtraction of sources coupling to a photon ("photonic sources"), are 
associated with semileptonic decays of HQ hadrons. As discussed in Sec. 13.31 the 
decay electrons largely preserve the modifications of the parent hadron spectra, 
albeit shifted in pt (by roughly a factor of ^^2). The more severe complication is the 
composition of the spectra, most notably the partition into charm and bottom 
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parentJI^. Since the heavier bottom quarks are, in general, less affected by the 
medium, their contribution significantly influences the resulting spectra. Unless 
otherwise stated, the calculations discussed below include a "realistic" input for the 
charm/bottom partition, i.e., either in terms of pQCD predictions ior p-p spectra or 
via empirical estimates from D-meson and electron spectra in p-p and d-An. Within 
the current theoretical and experimental uncertainties, both procedures agree, with 
an expected crossing of charm ad bottom electrons at pt — 3-6 GeV in p-p collisions 
at RHIC energy. 

Almost all of the approaches for computing HQ diffusion and/or energy loss 
introduced in Sec. [2] have been applied to data at RHIC. We organize the 
following discussion into (mainly perturbative) E-loss calculations (usually applied 
within a static geometry of the nuclear reaction zone) as well as perturbative and 
nonpertubative diffusion calculations using Langevin simulations for an expanding 
medium. 



3.5.1. Energy-Loss Calculations 

Radiative energy loss (E-loss) of high-energy partons in the QGP is believed to be 
the prevalent mechanism in the suppression of light hadrons with high pt > 6 GeV. 
It turns out that the application of this picture to the HQ sector (Sec. 12. 6p cannot 
account for the observed suppression in the non-photonic spectra. 

In the DGLV formalism, the high-pT suppression due to radiative E-loss of c 
and b quarks falls short of the data by about a factor of 3, cf. left p anel of Fig. 12'jp^. 
This led the authors to consider elastic E-loss (see also Refs.l^^^^i^^ which was found 
to be comparable to the radiative one out to the highest electron px measured thus 
far (^10 GeV). Their combined effect still underestimates the measured suppression 
by about a factor of ~2 for pt > 4 GeV. Similar findings were reported within the 
BDMPS approach: for a transport coefficient of q = 14GeV^/fns|, the spectra 
cannot be reproduced either, unless an unrealistic ass ump tion of neglecting the 
bottom contribution is ma de, cf. middle panel of Fig. I24P^. 

Both E-loss calculationJSSESlare performed for a static (time-averaged) medium 
of gluons, with fragmentation as the sole mechanism for hadronization. This is 
expected to be a good approximation at high px- Processes leading to an energy 
gain in the spectra, e.g., due to drag effects or coalescence with a light quark, are 
not included. Such processes lead to an increase in the pT of the final-state hadron 
and thus to an increase in the electron Raa, which would augment the discrepancy 
with data. The neglect of the diffusive term becomes particularly apparent in the 
elliptic ffow. In the E-loss treatment the only source of an azimuthal asymmetry 
in the pt spectra in non-central Au-Au collisions is the spatial geometry of the 
overlap zone: particles traveling along the short axis are less likely to be absorbed 

^This value is a factor 5-10 larger than pr edicted by pQCD, and at the upper limit of being 
compatible with light hadron suppressiorl22l. 
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Fig. 24. (Color online) PQCD E-loss cal culations fo r heavy quarks compared to single-electron 
observables in Au-Au collisions at Left panel: nuclear modification factor, 

^AA' ™ central collisions in the DGLV approacli^SJ (upper band: radiative E-loss only; lower 
band: radiative plus elastic E-loss) for a gluon plasma with rapidity density dN^ /dy=1000\ the 
bands represent the uncertainty in the leading-logarithm approximation of the elastic part (as 
described in Sec. I2.6| |. Middle and right panel: R\j^ in central and ii| in 0-80% central collisions 
within the BDMPS approach using a transport coefficient of g=14 GeV dashed and dash- 

dotted curves represent the results for c- and fe-quark contributions separately, while the solid 
curve is the combined result with a band indicating the pQCD uncertainty in the charm/bottom 
partitioJiSl. 



than those moving along the long axis of the approximately elliptic reaction zone. 
The positive V2 generated by this mechanism amounts to up to a few percent and 
significantly falls short of the observed electron V2, see right panel in Fig. [24l 

As an alternative mechanism, the collisional dissociation of D and B mesons 
from HQ f ragm entation in the QGP (Sec. l2.3."3|) has been implemented into an E-loss 
calculatiorF-4ij. A rather striking prediction of this calculation is that the shorter 
formation time of B mesons leads to stronger suppression than for D mesons above 
hadron momenta of pt — 15 GeV at RHIC, cf. left panels in Fig. [551 This turns out 
to be an important ingredient in the successful reproduction of the suppression 
data as shown in the right panel of Fig. 

3.5.2. Langevin Simulations 

The importance of clastic scattering for HQ diffusion and E-lo ss has been empha- 
sized, prior to quantitative measurements of spectra, in Refs.^^SUfi]^ albeit within 
rather different realizations of the underlying HQ interaction (recall Sees. l2.3Tl and 
I2.2.2[ respectively). 

The HQ spectra of the effective resonance + LO-pQCD modePS (using Langevin 

''The compatibility with the radiative picture of high-pt light-hadron suppression is presumably 
little affected since the Lorentz-dilated formation times of light quarks largely result in hadroniza- 
tion beyond the QGP lifetime. 
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Fig. 25. (Color online) Energy-loss calculations em ploy ing the coUisional dissociation mechanism 
for HQ fragmentation into heavy mesons in the QGP^^. Left panels: nuclear modification factor for 
D and B mesons in central 200 AGeV Au-Au and Cu-Cu at RHIC (upper panels) for gluon-rapidity 
densities of dJV^/djy = 1175 and 350, respectively, and in central 5.5 ATeV Pb-Pb collisions at LHC 
for dNS /dy = 2000 and 3500 (lower panels). Right panel : Ra a for e± in central 200 AGeV Au-Au 
collisions at RHIC compared to PHENIX and STAR datJSM the yellow bands reflect theoretical 
uncertainties due to an impact parameter expansion of the heavy-light quark interaction. 

simulations within an expanding fireball, recall Fig. ll9p have been converted into 
spectra utilizing a combined c oalescen ce/fragmentation scheme at Tc followed by 
heavy-meson three-body decayi^HMl. The predic ted spectra an d elH ptic flow 
show approximate agreement with 2005 PHENIX ESUMl ^nd STAPllSi data, see 
upper left and right panel of Fig. [26l respectively. Compared to the results for LO 
pQCD interactions only (blue lines), the resonance interactions (red bands) turn 
out to be instrumental in generating the required suppression and elliptic flow (see 
upper panels of Fig. [26]). LO-pQCD scattering alone, even with a strong couphng 
of as=0.4, does not produce sufficient coupling to the bulk medium to suppress the 
primordial quark spectra, nor to drag the heavy quarks along with the collective flow 
of the expanding fireball. The effect of heavy-light quark coalescence is illustrated 
by a calculation where only fragmentation is used as a hadronization mechanism 
(lower panels in Fig. [26]). In this case, the shape of the Raa and the magnitude 
of the V2 are not well reproduced. Coalescence processes add both momentum and 
V2 to the meson (and thus to the e^) spectra, i.e., the suppression becomes smaller. 
It is furthermore instructive to compare the LO-pQCD results with fragmentation 
only (lower left panel in Fig. [26]) to the pQCD E-loss calculations, especially to 
the elastic DGLV results where 0^=0.3 has been used (left panel of Fig. [24]) . The 
suppression level in the pertinent electron Raa is quite comparable for a rather large 
range in pT- The increasing trend in the Langevin calculations for pt ^ ^ GeV is 
presumably due to the dominant 5-quark contribution (which is barely suppressed 
even in the resonance model up to pr ^ 5 GeV, see left of Fig. [T9]) . Let us also 
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Fig. 26. (Color online) Nuclear modification factor (left panels) and elliptic flow (right panels) of 
spectra from heavy-flavor decays in Au-Au collisions at RHIC, as computed within HQ Langevin 
simu lations in a thermal fireball employing the resonance-|-pQCD model for HQ transport in the 
The upper panels include effects from heavy-light quark coalescence while the latter 
are switched off in the lower panels. The bands represent the full results for D- and _B-meson 
resonance widths of F = 0.4-0.75 GcV; the blue lines are obtained for LO-pQ CD interact ions only, 
while the purple lines neglect bottom contributions. The data are from Refs! 33|166 | 16 3 



estimate the impact of radiative contributions on the resonance model. Within 
DGLV the electron suppression due to radiative E-loss alone amounts to about 0.6- 
0.8 for pt — 4-10 GeV. Upon multiplying the Raa for the rcsonance-|-pQCD model 
in the upper left panel of Fig. [26] with this factor, the result would be compatible 
with current RHIC data. 

The PHENIX collaboration has conducted a comprehensive comparis on of their 
2006 e'^ datcP2l to theoretical calculations predicting both Raa and ^ j28 | 36 | 17 1J^ 
cf. left panel of Fig. [27l The interpretation reiterates some of the main points made 
above: (i) the missing drag in (radiative) E-loss calculations entails a substantial 
underprcdiction of the V2] (ii) Langevin calculations using elastic scattering require 
rather small HQ diffusion coefficients, Ds{2TrT) ~ 4-6, to be compatible with the 
observed level of suppression and elliptic flow, and, (iii) quark coalescence improves 
the simultaneous description of these two observables. 

The heavy-light quark T-matrix approach, based on input potentials estimated 
from thermal lattice-QCD, has been applied within the same Langevin-fireball -I- 



August 1, 2009 6:56 WSPC/INSTRUCTION FILE hq-qgp4 



60 R. Rapp, H. van Hees 




P,[Gev;c] p^[GeV] 



Fig. 27. (Color online) Left panel: Comparison of PHENIX datJ22J to theoretical predictions 
based on (i) Langevin simulations with the resonance+LO-pQCD interactions plus quark coales- 
cence (bandspH, (ii) Lan gevin simulations with (upscaled) HTL-improvcd LO-pQCD inte ractions 
(dash-dotted and dotted lines j23, and (iii) radiative E-loss calculations (dashed line s |551. For fur- 
ther illustration, PHE NIX da ta for tt^ sup pression and elliptic flow arc 

showjnni. Right panel: 

Comparison of and STAH^ data to Langevi n-fir eball simulations for HQ trans- 

port using T-matrices with IQCD-based heavy-quark potentialJ3Sl an uncertainty due to different 
IQCD internal energies is indicated by the solid and dash-dotted curve; the dashed lines are 
obtained without the effects of heavy-light quark coalescence at Tc ■ 



coalescence/fragmentation scheme as the effective resonance modeP^. The perti- 
nent spectra (cf. rig ht piancl o f Fig [27]) exhibit a comparable level of agreement 
with current RHIC datiSlIMlSi 

as in the resonance model, with a similar uncer- 
tainty due to different extractions of the HQ internal energy. Although the T-matrix 
calculations involve, in principle, no tunable parameters, the inherent theoretical 
uncertainties are appreciable (e.g., in the definition of the in-medium potential in 
terms of internal or free energy). Let us, however, recall a rather general feature of 
the T-matrix approach which was visible already at the level of the HQ spectra in 
Fig. 1211 the weak temperature dependence of, e.g., the friction coefficient implies 
that the HQ coupling to the medium remains rather strong in the later QGP and 
mixed phase stages of the evolution. Since the bulk V2 is largest in these later stages, 
while the suppression largely occurs in the first l-2fm/o^, the T-matrix interac- 
tions generate relatively more V2 than suppression compared to, e.g., the resonance 
model (or, alternatively: for the same V2, the suppression in the T-matrix approach 
is smaller). This traces back to the increasing color-Debye screening with increasing 
temperature, which leads to a gradual melting of the resonance correlations and a 
marked increase of the spatial diffusion constant, Ds/(2ttT), with temperature (re- 
call Fig. I14|) . Such a temperature dependence appears to improve the consistency 
in the simultaneous description of the Raa and V2, but more precise data are 
needed to scrutinize this feature (including d-A collisions to quantify the Cronin 
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Fig. 28. (Colo r on line') Boltzmann-transport model results for electron Raa (upper panels) and V2 
(lower panels |S3 in a hvdrodvna mic ev olution using different versions of LO-pQCD HQ interac- 
tions, compared to PHENIX datel^SESlin 200AGeV Au-Au collisions at RHIC. The curves in the 
left panels are computed with fixed Os (SttT) at given temperature, conventional infrared regulator 
(/i|j = r/i|j with r=re=l) and large iC-factors; the curves in the right panels are computed with a 
running a^, reduced IR regulator (r=0.2) and reduced K factors of 1.5-2 (or 2-3), as represented 
by the bands, ttrmax or ttrmin indicate frcezcout at the beginning or end of the mixed phase, 
respectively. 



effect, which could increase the Raa without noticably affecting U2). 

Finally, we reproduce in Fig. [2H1 selected results of the Boltzmann-transport 
approach of HQ diffusiorPS with the background medium described by the hy- 
drodynamical model of Ref.El hadronization is treated in a combined coales- 
cence-|-fragmentation approach similar to the one in Refs.^^^SSl. The HQ inter- 
actions in the QGP are implemented via the elastic pQCD scattering amplitudes 
described in Sec. 12.2.31 The left panels in Fig. [55] refer to a model with fixed cou- 
pling constant, as(27rT), at given temperature and standard screening mass (r=l). 
A large K factor of A' = 12 is needed to simultaneously reproduce the electron 
Raa V2 data. The final elliptic flow is found to be rather sensitive to the late 
QGP stages of the evolution, favoring hadronization at the end of the mixed phase 
(i.e., a t a small transition energy-density, Ctrans); this is consistent with the findings 
of Ref.lSSl and the above discussion of the T-matrix interaction. It thus corrob- 
orates that anisotropic matter-flow can only be transferred to the heavy quarks 
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if the latter is sufficiently large, while E-loss (reflected in high-p^ suppression) is 
mostly effective when the fireball density is high. The simulations find little impact 
of the initial-state Cronin effect on the final heavy-quark V2, but the suppression 
is somewhat reduced primarily for momenta of pt — 1-3 GeV. The right panel 
of Fig. [28] illustrates that similar results can be achieved with smaller /C-factors, 
i^=1.5-2, if the pQCD cross sections are augmented by a running coupling, a^it) 
[t: 4-momcntum transfer in the elastic scattering process), and a small infrared 
regulator, jl^ = with r=0.2, in t-channel gluon-exchange scattering. 



3.6. Viscosity? 

In this section we utilize the quantitative estimates for the HQ diffusion coefficient 
as extracted from current RHIC data to obtain a rough estimate of the ratio of 
shear viscosity to entropy density, 77/s, in the QGP. This quantity has received 
considerable attention recently since (a) it allows to quantify deviations from the 
predictions of ideal fluid dynamics for observables like the elliptic flow, and (b) 
conformal field theories in the strong coupling limit are conjectured to set a universal 
lower bound for any liquid, given by 77/s = l/(47r|I3, referred to as KSS bound. In 
the following we will bracket the estimates derived from HQ observables by using 
relations of Ds and 77/s in the weak- and strong-coupling limit (the latter assumed 
to be given by the AdS/CFT correspondence). 

Following the discussion in Sec. 12.41 the strong coupling limit in the AdS/CFT 
framework results in a (spatial) HQ diffusion constant of ~ 1/ {2ttT). Combining 
this with the lower bound of the viscosity-to-entropy-density quoted above, one 
obtains 

- = . (85) 

s 2 

For a weakly coupled dilute gas, an estimate for a relation between Dg and 77 can 
be inferred starting from the calcu lation of the shear viscosity from kinetic theory 
for an ultrarelativistic gaJI^^IIZSl^ 

4 

77 « —n (p) Atr , (86) 
15 

where n denotes the particle density, {p) the average momentum of the gas particles 
and Atr its transport-mean free path. With the estimate, n (p) ~ e, for the energy 
density and the corresponding equation of state, Ts — e + P — 4/3e, one arrives at 

1 , ^ 

- ^ -TXu- ■ 87 
s 5 

Finally, equating the transport mean-free path Atr to the mean-free time Ttr and 
taking into account the delay due to the mass effect of the heavy quark on the 
thermalization time, tq w Tt^T/niQ, one finds 

^ « ItD, . (88) 
s 5 
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Fig. 29. (Color online) The ratio of shear viscosity to entropy density, rj/s. Left panel: schematic 
estimates using charm-quark diffusion constants based on (a) schematic LO pQCD elastic scatter- 
ing (cis=0.4) in the weakly interacting limit, Eq. I I88I I (dashed line), (b) pQCD elastic scattering 
with running coupling constant and small IR regulator (band enclosed by dash-dotted lines using 
the weak- and strong-coupling limits), (c) the effective resonance -|- pQCD model in the strong- 
coupling limit, Eq. I)86|l (band enclosed by long-dashed lines for r=0.4-0.75 GeV), and (d) the 
lattice-QCD potential based T-matrix approach augmented by pQCD scattering off gluons (band 
enclosed by solid lines constructed from th e wea k- and strong-coupling limits). Right panel: lat- 
tice QCD com putations in a gluon plasmJSll compared to results inferred from perturbation 
theorJmS].^ 



In comparison to the "strong-coupling" estimate within AdS/CFT, Eq. the 
shear viscosity appears to be underestimated when the kinetic theory for a dilute 
gas is applied to liquids. These estimates are now applied to several of the HQ 
diffusion calculations discussed above, see the left panel of Fig. [29l Since ry/s cx 
Ds{2ttT), the main features of Fig.[T4]are transmitted to 7]/s, in particular the weak 
temperature dependence of the LO-pQCD calculations and the effective resonance 
model. Of course, the absolute values of these calculations differ considerably. A 
different behavior is only found for the T-matrix -|-pQCD model, which suggests a 
transition from a strongly coupled regime close to Tc to relatively weak coupling 
above ^2Tc. In fact, the uncertainty band has been constructed as follows: for the 
lower limit; the weak-coupling estimate Eq. (|88p is used; for the upper limit, the 
strong-coupling limit estimate, Eq. ([86|) . at r=:0.2GeV is linearly interpolated with 
the LO-pQCD weak-coupling limit at r=0.4GeV (the strong-coupling estimate for 
T-matrix-|-pQCD overshoots the LO-pQCD result at this temperature). As for the 
spatial diffusion constant, the increase oi-q/s with temperature is related to color- 
Debye screening of the IQCD-based potentials which entails a gradual melting of the 
dynamically generated resonances in the heavy-light quark T-matrix. It is tempting 
to interpret the decrease of rj/s when approaching Tc from above as a precursor- 
phenomenon of hadronization and thus connected to the phase transition itself. It 
remains to be seen whether a similar mechanism is operative in the light-quark 
and/or gluon sector (three-body interactions are unlikely to produce this due to 
the decrease in particle density when approaching Tc from above). Such a behavior 
is rather general in that it has been observed around phase-transition points for a 
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large variety of substances, see, e.g., the discussion in Refs.^^^^^^^. 

Fi nally we show in the right panel of Fig. [29] a quenched IQCD computation 
of r;/Jl^. The error bars are appreciable but the results tend to favor rj/s values 
which are below LO-pQCD calculations. The specific pQCD result inclu ded i n this 
plot employs a next-to- leading logarithm calculation for the shear viscositjff^ and a 
self-consistent hard-thermal-loop calculation for the entropy densitjEH. It is rather 
close to the schematic LO calculation (using as=QA) in the left panel of Fig. [29l 



4. Heavy Quarkonia in Medium 

In recent years it has become increasingly evident that observables in the heavy- 
quarkonium and open heavy-flavor sectors are intimately connected. In the orig- 
inal picture of charmonium suppression as a probe of color-screening in hot and 
dense QCD matteJI^ there are no obvious such connections. Several recent de- 
velopments have changed this situation. Thermal lattice QCD calculations find 
that charmonium correlation functions are remarkably stable up to temperatures 
of ^2Tc or higher, suggestive for the survival of the ground state (77c, J/V") well 
into the QGP. This interpretation is supported by probabilistic extractions of the 
pertinent quarkonium spectral functions. It implies that quarkonia can not only 
dissociate but also regenerate in the It immediately follows that the yield 

and spectra of regenerated quarkonia are, in principle, sensitive to the abundance 
and momentum spectra of open-charm states in the system. E.g., for a fixed total 
charm number in the system, a softening of the heavy-quark spectra is expected to 
increase c-c overlap in phase space and thus enhance the probability for charmo- 
nium formation. At the same time, elliptic flow of charm quarks will imprint itself 
on regenerated charmonia. Furthermore, HQ interactions with light quarks (and 
possibly gluons) may be closely related to the interaction (or potential) between 
two heavy quarks. E.g., the T-matrix approach discussed in the previous section 
is directly based on potentials which are extracted from the HQ free energy com- 
puted in lattice QCD. As we argued there, this approach to evaluate HQ diffusion 
has several attractive features, both theoretically (it may provide maximal inter- 
action strength in the vicinity of Tc) and phcnomenologically (it describes current 
HQ observables at RHIC fairly well). 

In the remainder of this section we address several aspects of quarkonia in 
medium and in heavy-ion collisions with a focus on connections to the open heavy- 
flavor sector. More extensive reviews on quarkonia in medium have recently been 
given m Refs.m which we do not attempt to reproduce here. In Sec. l4.1l we give 
a brief review of theoretical issues in the understanding of in-medium quarkonium 
spectral properties, in terms of thermal lattice QCD results for correlation and spec- 
tral functions and their interpretation using effective potential models (Sec. I4.LT|) . 
The latter are employing input potentials extracted from heavy-quark free ener- 



'Note that higher dissociation rates also imply higher formation rates. 



August 1, 2009 6:56 WSPC/INSTRUCTION FILE hq-qgp4 



Heavy quarks in the quark- gluon plasma 65 



gics computed in lattice QCD. thus enabling, in principle, an internal consistency 
check, provided a suitable potential can be defined. While color screening is a key 
medium effect in the potentials (governing the binding energy of the bound states), 
a quantitative assessment of spectral functions requires the inclusion of finite- width 
effects induced by dissociation reactions and possibly elastic scattering (Sec. 14.1^ . 
In Sec. 14.21 we elaborate on recent developments in describing heavy-quarkonium 
production in heavy- ion collisions. The main focus is on transport models which 
track the dissociation and regeneration of charmonia (and bottomonia) through the 
QGP, mixed and hadronic phases fSec. I4.2?T|) . complemented by a brief discussion 
of initial conditions as affected by cold-nuclear-matter (CNM) effects (shadowing, 
Cronin effect and nuclear absorption). This is followed by an assessment of the 
current status of charmonium phenomenology at SPS and RHIC. 



4.1. Spectral Properties of Quarkonia in the QGP 
4.1.1. Lattice QCD and Potential Models 

iQiri 

The phenomenological Cornell potential^ for the interaction between two heavy- 
quark (color-) charges in the color-singlet channel, 

yQQ{r;T = Q) = -~+ar , (89) 

has been very successful in reproducing the vacuum spectroscopy of charmonium 
and bottomonium bound states. It consists of a (color-) Coulomb plus a (hnear) 
confining part with a very limited number of parameters, i.e., a strong coupling 
constant, as, and string tension, a (in addition, an effective HQ mass, toq, needs 
to be specified). Subsequent developments have put this framework on a more 
rigorous footing by showing that (a) the potential desc ription can be recovered as 
a low-energy effective theory of QCD with heavy 

quarkJSSm^ and (b) lattice QCD 
computations of the color-singlet heavy-quark free energy, Fi , have found exc ellent 
agreement with the f unct ional form (and parameters) of the Cornell potentiaff^. 

Early calculationJISI! of Q-Q bound-state properties in the QGP have supple- 
mented the Cornell potential by a phenomenological ansatz for color screening of 
both the Coulomb and confining parts, 

tiniT) V ^ J ir 

The key quantity carrying the temperature dependence is the Dcbye screening mass, 
liD (oc gT in thermal pQCD). Already at that time the possibility was established 
that ground-state charmonia (and even more so bottomonia) can survive until tem- 
peratures (well) above Tc- More recently, quantitative IQCD computations of the 
finite-temperature color-singlet free energy of a HQ pair, Fi(r;T), have become 
available, see, e.g.. Fig. [30l The results nicely illustrate the color-screening effect 
and its gradual penetration to smaller distances. When inserting the in-medium 
free energy as an improved estimate of the finite-temperature HQ potential into a 



^Oq(^; T) = 7-?^ (1 - - ^ e---(^)'^ . (90) 
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Fig. 30. ( Colo r online) Free e nergy of a static color-singlet HQ pair as computed in lattice QCD 
for Nf = ^23 (left) and Arj=ilS2l flavors (right). The critical temperature is Tc = 193(170) MeV 
for the Nf = 3(2) calculation, and the string tension typically amounts to (t^/^ ~ 420 MeV = 
1/(0.45 fm). 



Schrodingcr equation, t he "m elting" temperature of the J/^p Xc) was found to 
be just above (below) tJI^SI. 

Further progress in thermal IQCD came with the eomputation of heavy quarko- 
nium correlation functions, 

Ga(r,r)-(0„(r,r)jt(0,O))) (91) 

(also referred to as temporal correlators), as a function of Euclidean time, r. j„ rep- 
resent the creation/annihilation operators of a hadronic current of given quantum 
numbers, a. In the pseudoscalar and vector charmonium channels (corresponding to 
c-c iS- waves with 77c and J/ip states, respectively), the Euclidean correlators were 
found to exhibit a surprisingly weak temperature dependence up to ~2Tc, even 
at large r, suggestive for rather stable bound states. The temporal correlators are 
related to the physical spectral function, (7aiE,p',T), via 

00 

(r, p; T) = J dE a„ {E, p; T) K{E, r; T) (92) 


with a thermal integral kernel 

cosh[i;(r - i/2r)] 
sinh[£;/2T] • 

Eq. ([nil) implies that the extraction of the spectral function from the Euclidean 
correlator requires a nontrivial integral inversion. Especially at finite T, where pe- 
riodic boundary conditions limit the information on Ga{T,p] T) to a finite interval, 
< r < 1/T, and for a finite number of r points, the unambiguous inversion 
to obtain aa{E,p]T) becomes an ill-defined problem. However, using probabilistic 
methods (in particular the so-called maximum entropy method (ME M)), a sta- 
tistical reconstruction of aa{E,p;T) is possible and has been appliecP^lSl^ xhe 
approximate constancy of the temporal correlators lead to spectral functions with 
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Fig. 31. HQ free energy in the color-singlet channel computed in thermal A'^ 
and corresponding numerically extracted internal energy (right 



2 lattice QCD 



rather stable ground-state peaks corroborating the notion of surviving ground states 
well above Tc. 

To resolve the apparent discrepancy with the low dissociation temperature found 
in the potential model discussed above, it has been suggested to employ as potential 
the internal rather than the free energy, which are related via 



F{r-T)^U{r-T)-TS{r-T) 



(94) 



Especially in the color-singiet channel, the (positive) entropy contribution rises sig- 
nificantly with Q-Q separation, r, thus producing "deeper" potentials (cf. Fig. W\\ 
entailing stronger binding. Consequently, pertinent evaluations of quarkonium 
spectra lead to larger dissociation temperatures, which seemingly agree bet- 
ter with t he IQCD spectral func tions. These assertions have been made more 
quantitativj5M58|59|60|61|62|63| 

by employing potential models to calculate in- 
medium spectral functions, perform the straightforward integral in Eq. (|92p and 
compare to the rather precise temporal correlators from IQCD. It is important 
to realize that the Eu clid ean correlators involve the pertinent spectral function 
at all energies, the in-medium bound-state spectrum obtained from a 

Schrodinger equation (using either a screened Cornell potential or IQCD internal 
energies) has been combined with a perturbative ansatz for the Q-Q continuum 
abo ve thre shold. No agreement with IQCD correlators could be established. In 
jigfg 157158] importance of rescattering effects for the interacting Q-Q continuum 
was emphasi zed and implemented into the calculations of the Q-Q spectral func- 
tions. In continuum correlations were impl eme nted via Gamov resonance 
states in Breit-Wigner approximation, while in a thermodynamic T-matrix 
approach was employed, 

d^fc 



(27r)= 



Gqq{E; k) T^{E) [1 - /q«) - fgiu)^)] 



(95) 



exactly as introduced in the context of HQ diffusion in Sec. 12.3.21 recall Eq. (j46|) . 
The T-matrix approach enables a consistent treatment of bound and continuum 
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internal HQ free energies extracted from A^^ = 3 thermal lattice 

qceIS!. 

For numerical purposes, 

a small charm-quark width of 7c = 20 MeV has been implemented into the intermediate QQ 
propagator in the scattering equation l|46| l. 



states on equal footing, as well as the implementation of medium effects (selfener- 
gies) into the intermediate two-particle propagator, Gqq, recall Eq. (|48|) . Pertinent 
results for S- and P-wavc charmonium spectral functions, using the internal en- 
ergy extracted from the Nj = 3 free energjl^Sl (left panel of Fig. [5D|) . are shown in 
Fig. [2H for a constant (T-independent) charm-quark mass of nic = 1.7 GeV. One 
clearly recognizes the reduction in binding energy as a result of color screening by 
the upward moving bound-state peak position with increasing temperature (the 
cc threshold is fixed at 2mc = 3.4 GeV). Also note that nonperturbative rescat- 
tering effects close to and above threshold induce a substantial enhancement in 
the Q-Q spectral function over the non-interacting continuum (indicated by the 
red long-dashed lines), an effect which is of prime importance in the calculation 
of HQ diffusion in q-Q scattering as well. When applied to the calculation of Eu- 
clidean correlatorgSSI in Eq. (|92|) , the upward shift of low-energy strength due to 
the moving bound states in the S'-wave spectral function shown in the left panel 
of Fig. [32] entails a suppression of G(r) with temperature which disagrees with the 
weak temperature dependence found in IQCD. 

Another important ingredient to understand the behavior of the correlators is 
the temperature dependence of the HQ mass. Schematically, the, say, Jf-ip bound- 
state mass may be written as 

mj/i^ = 2"i* - es • (96) 

This illustrates that a small (large) binding energy, es, can be compensated by 
a small (large) effective quark mass in a way that the ground-state mass stays 
approximately constant. Indeed, when interpreting the asymptotic value of the in- 
medium potential as an effective HQ mass correction, 

m*=m2 + Amc , ArUc = X{r = oo;T)/2 , (97) 

with X = U or F (or an appropriate combination thereof), the use of U im- 
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E_._^_ [GeV] 

Fig. 33. (Color online) S-wave cliarmonium spectral functions {rjc or J/V') computed employ ing: 

(i) a screened Cornell potential within nonrelativistic Green's function ap proa ch and 

(ii) an IQCD-based = 3 internal energ Pl within a T-matrix approacM Both calculations 
account for in-medium charm-quark masses, and approximately reproduce the weak temperature 
dependence of the temporal correlators computed in IQCD up to temperatures of at least 2Tc 
(see, e.g., inset in the left panel). However, in calculation (i) the bound-state has disappeared at 
temperatures below 1.5Tc, while in calculation (ii) it is visible up to ~2.5Tc. 



plies strong binding with large effective quark masses while the use of F leads to 
weak binding with small to* (recall from Fig. [31] that the "{/-potential" is deeper 
than the "i^-potential" but features a larger asymptotic value for r — > cxd). Conse- 
quently, it has been found that reasonable agreement with IQCD correlators can be 
achieved with different spectral functions, covering a rather large range of dissocia- 
tion temperatures, e. g., sligh tly above Tc using screened Cornell potentials similar 
to the free en ergy using a linear combination of free and internal 

gj^gj.gji57163|^ or up to ~2.5Tc when using internal energieJ^^, cf. Fig. [33l 

To resolve this redundancy, it will be necessary to develop independent means 
of determining the in-medium quark mass and the appropriate quantity to be iden- 
tified with the HQ potential. First estimates of the HQ mass from thermal IQCD 
have been obtained by approximating the HQ number susceptibility within a quasi- 
particle model with effective quark massS^. The results suggest a rather moderate 
temperature variation of the latter, w hich devia tes significantly from the pcrturba- 
tive predictions up to T~3Tc. In Refs.^^ ^^ l ^"^ ! hard-thermal-loop and HQ effective 
theory techniques have been applied to derive the leading terms in a perturbativc 
and HQ mass expansion of a finite temperature potential. An interesting finding of 
these investigations is that the potential develops an imaginary part in the medium 
which arises from the Landau damping of the exchanged gluons, representing a 
decay channel of the HQ bound state. A more general discussion of the in-medium 



■■We do not address here the issue of so-called zero-mode contributions to quarkonium correla- 
tors, which arise on the lattice due to the periodic boundary conditions in temporal directioriS^. 
These contributions are essential to obtain quantitative agreement with the IQCD correlators in 
all mesonic quantum- number channels except the pseudoscalar one (»?c); they arc rather straight- 
forward to implement in quasi-particle approximation. 
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decay width of heavy quarkonia. which plays a central role for phenomenology in 
heavy-ion collisions, is the subject of the following section. 

The impact of finite-width effects on charmonium correlators has been studied 
within the T-matrix approach in RefP, by implementing an imaginary part into 
the charm-quark propagators. A broadening of charmonium spectral functions leads 
to an enhancement of the temporal correlators (due to additional strength at lower 
energies), which, however, is only a few percent for a charmonium width on the order 
of r$~100MeV. On the one hand, such a value for the width is phenomenologically 
significant, as it implies that about 60% of the charmonia decay within a time of 
2 fm/c. On the other hand, for larger widths, their impact on the correlators should 
be accounted for in quantitative comparisons to IQCD "data" . 

4.1.2. Dissociation Widths 

The spectral width of a quarkonium state propagating through matter can, in prin- 
ciple, receive contributions from elastic and inelastic reactions with the medium 
particles. Elastic scattering affects the momentum distribution of the quarkonium 
while inelastic interactions change its abundance (via dissociation or formation). 
More formally, the quarkonium acquires a complex selfenergy which can be ex- 
pressed via the in-medium scattering amplitude, A4^i, folded over the (thermal) 
distribution, fi, of the medium particles. 



The real part of characterizes in-medium changes of the quarkonium mass 
while the imaginary part determines its width, r$(£^)=-2 Im I]^(£'j^. Most of the 
attention thus far has been directed to the inelastic reactions (rather than elastic 
scattering). Using the optical theorem to relate the imaginary part of the forward 
scattering amplitude to the cross section, one arrives at the well-known expression 



where Wrci denotes the relative velocity of the incoming particles and s=(p + fc)^ the 
squared center-of-mass energy of the ^'-i collision. The first evaluation of t he inelas- 
tic quarkonium reaction cross section with gluons was conducted in 
Employing Coulomb wave functions for the quarkonium bound state, the leading- 
order process, ^ + <? ^ Q + Q, is the analog of photo-dissociation of hydrogen, see 
left panel of Fig. [34l For an S'-wave ^ bound state with binding energy eb, the 




(98) 




(99) 



'^In addition, mass and widtii changes arc induced at the Q-Q level via in-medium effects on the 
Q-Q potential and direct i]/ — > Q-(-Q decays, respectively. These effects can be accounted for, e.g., 
in the underlying Q-Q T-matrix, Eq. I|95| l. 
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Fig. 34. Di agrams f or quarkonium-dissociation reacti ons via parton impact; left panel: gluo- 
dissociation CSSHHSl right panel: quasifroe dissociation H3. 



cross section is given as a function of incoming gluon energy, fcg, by 

1(1) 7!^) J, w.°-f\ ,i„„, 

3 V 3 / \eB J mf^ (fco/es)^ 

This expression has a rather pronounced maximum structure with the peak cross 
section reached for a gluon energy, fc™''^ = l^s- ^^'^ applicability of the gluo- 
dissociation formula should be reasonable for the free bottomonium ground state 
(ex — 1 GeV), but borderline for J/t/; {ej/^ ~ 0.6 GeV). Taken at face value 
for QGP temperatures of T = 300-400 MeV, where the typical thermal energy of 
(massless) gluons is around ko ~ I GeV, the convolution of the gluo-dissociation 
cross section with a thermal gluon distribution function in Eq. (j99p results in an 
inelastic J/ip width (lifetime) of Tj/^ ~ 150-400 MeV (tj/^ ~ 0.5-1.3 fm/c), see 
the dashed line in the right panel of Fig. 

The situation changes if the quarkonium binding energy decreases due to color- 
screening as discussed in the previous section (or for excited charmonia which are 
weakly bound even in vacuum). In this case, the peak of the gluo-dissociation cross 
section moves to smaller energies and becomes rather narrow; the loss of phase space 
can be basically understood by the fact that for a loosely bound ^' state, the absorp- 
tion of an on-shell gluon on an (almost) on-shell quark is kinematically impossible 
(suppressed). Consequently, with decreasing biiiding energy, the cross section has 
less overlap with the thermal gluon spectrum SZGHZT leading to a decreasing width 
with temperature (cf . dotted line in the left panel of Fig. [35]) . This unphysical be- 
havior signals the presence of other inelastic processes taking over. In Ref.H3 the 
so-called "quasifree" dissociation mechanism has been suggested, J/'4'+P c+c+p, 
where a thermal parton {p = g, q, q) scatters "quasi-elastically" off an individual 
heavy quark in the bound state (see right panel of Fig. [34]). The p-Q scattering 
amplitude has been evaluated in leading-order (LO) perturbation theory in- 
cluding thermal parton and Dcbyc masses and slightly modified kinematics due to 
the small but finite binding energy. While naively of next-to-leading order (NLO) 
in Qg compared to gluo-dissociation, the additional outgoing parton opens a large 
phase space rendering the quasifree process significantly more efficient for weakly 
bound states. Therefore, it readily applies to excited states as well, thus enabling 
a treatment of all charmonia on an equal footing (which is essential even for J/ip 
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Fig. 35. (Color online) Parton-induced dissociation widths for J/ip at rest in a QGP as a function 
of temperature. Left panel: gluo-dissociation width based on the cross section, Eq. IIIOOI I. with 
free (dashed line) and in-medium decreasing binding energy (dotted li ne'), as w ell as "quasifree 
dissociation" width (J/ip + p —> c + c + p) with in-medium decreasing e jSHSZl. Right panel: full 
NLO calcu lation for (a) quark- and (b) gluon-induced dissociation for different thermal parton 
masseiSHl. 



observables, since the latter receive significant feed-down contributions from Xc and 
tfj' states, see Sec. 14.2.1] below). For a coupling constant of as ~ 0.25, the quasifree 
dissociation rate reaches Fj/^ = 100-200 MeV for temperatures, T = 300-400 MeV 
(cf. soHd line in the left panel of Fig. [35]). These values could be substantially en- 
hanced if non-perturbative p-Q scattering mechanisms are operative, much like the 
ones discussed in Sec. 12.31 A complete NLO calculation for parton-induced char- 
monium destruction has recently been carried out in Ref.Q^S', including the effects 
of in-medium reduced binding energies. The right panels of Fig. [35] show the re- 
sults for quark- and gluon-induced breakup of the J/^p for as = 0.5 and different 
(fixed) thermal parton masses. For temperatures around 250 MeV, the sum of both 
contributions (Fj/^ ~ 350 MeV) is about a factor ^4 larger than the quasifree 
results (Fj/^ ~ 80 MeV) in the left panel, calculated for as = 0.2j|^. Thus, there 
is good agreement between these two calculations, since the rate is basically oc 
(for T — 250MeVthe T-dependent Dcbye-mass in fj,D = gT, amounts to 

~ 440 MeV). 

The three-momentum dependence of the dissociation rate of a moving quarko- 

niunij_r^_(V) , has been calculated for full NLO and the quasifree rates in 
Refe.[189ll90]^ 

respectively. In both calculations a weak increase of the rate with 
increasing three-momentum is found. Since the quasifree cross section is essentially 
constant, this increase is caused by the increasing flux of thermal partons encoun- 
tered by the moving bound state. A simil ar re sult has been obtained in a calculation 
employing the AdS/CFT correspondenci^^ (recall Sec. l2.4l for more details on this 
framework and its caveats). On the other hand, gluo-dissociation leads to a rather 
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pronounced decrease of the dissociation rate with increasing three-momentum, since 
the pertinent cross section is peaked at relatively low energies and falls off rapidly at 
large center-of-mass energies, s = (p + k)"^ . Of course, if gluo-dissociation becomes 
ineffective, its three-momentum dependence becomes immaterial. A decreasing p- 
dependence should also be expected in quasifree dissociation if non-pcrturbative 
(resonance-like) q-Q interactions are operative (recall Sec. 12. 3p . since the pertinent 
cross sections are concentrated at small ^/s as well. 

To utilize quarkonium observables as a probe of QGP formation in URHICs, it 
is mandatory to have good control over the modifications of quarkonia in hadronic 
matter, in particular their inelastic reaction rates. From current lattice QCD cal- 
culations it is very difficult to extract information on excited states (say, In 
addition, the results in the P-wave (scalar and axialvector) channels (corresponding 
to Xc.o a-nd Xc.i, respectively) are sensitive to the so-called zero- mode contributions 
which are not directly related to the bound-state properties (as briefly mentioned in 
a previous footnote). Potential models find that ij}' ^ Xc,o and Xc,i "melt" close to or 
even below Tc, suggesting substantial modifications in the hadronic phase. Even for 
the J /ip, hadronic dissociation may lead to significant suppression (in addition to 
suppressed feed down from i})' and Xc states). One main obstacle in a reliable assess- 
ment of these reactions is that low-energy reactions of the type h+J/tfj D+D+X 
constitute a nonperturbative problem with little experimental information available 
to constrain effective models. 

An initial estimate of quarkonium dissociation by light hadrons has been ob- 
tained by using the gluo-dissoci ation cross section, Eq. (jfOOp . convoluted over the 
gluon distribution inside hadronJI^. Since the latter is rather soft (fco ~ O.I GeV), 
the gluon energy is in general not sufficient to break up the J/ip, leading toa (low- 
energy) cross section of order o-™f/,^ — O.Imb. Quark-exchange reactionJI^, e.g., 
in meson-induced breakup, h + J/tp D + D (including excited D mesons in the 
final-state), are presumably more relevant. Effect ive quark models predict dissocia- 
tion cross sections of order 1-2 mb, see, e.g., Refs.^l^^lUS] An a lternative approach 
is to construct effective hadronic models, pioneered in Guiding principles 

are basic symmetries including gauge invariancc for vector mesons { J/ip^ p) as well 
as flavor symmetries, most notably SU(4) (albeit explicitly broken by the charm- 
quark mass) and chiral symmet ry which is operative in interactions with (pseudo- 
) Goldstone bosons (tt and j^j l97 | I98 | I99 | 200 | 20I | 20 2l -pj^^ ^^.^ uncertainty in 

these models remains a reliable determination of the cutoff scales figuring into the 
hadronic vertex form factors. With cutoff values of around I GeV, the agreement 
with quark models is quite reasonable; dissociation reactions induced by p mesons 
appear to be the most important channel. Their thermal density in hadronic matter 
(i.e., for temperatures of ~I80MeV) is not very large, so that the total J/ip width 
does not exceed a few MeV, and therefore is substantially smaller than in the QGP. 
E.g., for a total hadron density of Qh = 3^po and a thermally averaged cross section 
of (cjjj/^) = I mb (corresponding to significantly larger peak cross sections), a 
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rough estimate for the dissociation rate gives T^'J^^ = ('''^j/^ ^rci) 9h — 5 McV. 
An interesting possibihty to constrain efFe ctiye had ronic vertices in a rather model- 
independent way is to use QCD sum rulePM. Pertinent estimates yield, e.g., 
a thermally averaged 7rJ/-0 dissociation cross section of (cjj'jy^ Vi-ei) = 0.3 mb at 
T = 150 McV, in the same range as the above estimate. We finally remark that 
in- medium effects, e.g., modified spectral distributions of Z?- mesons, can increase 
the final-state phase space and lead to an appreciable increase of the dissociation 
ratePQ^l, This is especially pertinent to excited charmonia like the ip' , whose mass 
is close to the free DD threshold, so that a slight reduction (or b roadening) in the 
ZJ-meson mass can open the direct decay channel, ip' £)jj 205 | 



4.2. Quarkonium Production in Heavy-Ion Collisions 

Similar to the open heavy-flavor sector, a key objective (and challenge) in the 
quarkonium sector is to connect their equilibrium properties to obscrvables in 
heavy-ion collisions, and eventually deduce more general insights about basic prop- 
erties of QCD matter, e.g., color de-/ confinement and Dcbye screening. Since 
quarkonium states in heavy-ion collisions are even more rare than individual heavy 
quarks, it is suitable to adopt a transport treatment for their distribution functions 
in a realistic "background medium" whose evolution is not affected by the heavy 
quarks or quarkonia. The connection between observables extracted from the dis- 
tribution function (after its "transport" through the medium) and the equilibrium 
properties discussed in the previous section is given by the coefficients and equilib- 
rium limit of the transport equation, as elaborated in the following section, 14.2.11 
The current status in comparing various model implementations to charmonium 
data at SPS and RHIC will be assessed in Sec. 14.2.21 



4.2.1. Quarkonium Transport in Heavy-Ion Collisions 

The (classical) Boltzmann equation describing the time-evolution of the phase-space 
distribution function, /^(r, r;p), of an (on-shell) quarkonium state, ^' (with energy 
Po ~ ujp = {p^ + rn^y/'^), may be written as 

p^df,U{r,T]p) = -Up r>t(r,T;p) f^{r,T]p) + ujp (3^{r,T;p) (101) 

(a mean-field term has been neglected assuming that the real part of the ^ self- 
energy is small). When focusing on inelastic reactions, r$(r,r;p) represents the 
dissociation rate discussed in Sec. 14.1.21 above, which governs the loss term, i.e., 
the first term on the right-hand-side {rhs) of Eq. (|10ip . The (r,r) dependence of 
r$ typically converts into a temperature dependence via the fireball evolution of 
a heavy-ion reaction for given projectile/target (A/B), collision energy (-y/s) and 
impact parameter (6). The second term on the rhs of Eq. (jlOip is the gain term 
accounting for the formation of quarkonia. For a 2 ^ 2 process (as, e.g., realized 
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via the inverse of gluo-dissociation, Q + Q g + ^'), it takes the fornPSSl 

xM^^*(s) e[rdi..-r(r,T)] {2iTf 5^'\p + k - pq - pa) . (102) 



To ensure detailed balance, the cross section figuring into the formation probabil- 
ity, W^q{s) = cr^QQ^g^ ^i:c\'^^pq^pq^ li^s to be the same (up to a kinematic and 
statistical factor) as the one used in the dissociation rate, Eq. For reactions 

beyond 2 2 (such as the quasifree process, p + Q-\-Q^p-\-'^), the microscopic 
evaluation of the gain term becomes more involved. Note the explicit dependence 
on the HQ phase-space distribution functions, f'^ '^, in Eq. (|102[) . whose modifi- 
cations in heavy-ion reactions are the central theme in Sec. [3] of this review. The 
temperature-dependent step function in Eq. (|102[) signifies the limit set by the dis- 
sociation temperature, Tdiss, above which a well-defined ^' state no longer exists 
and thus formation reactions arc not meaningful. 

It is both instructive and useful for practical applications to simplify the gain 
term by integrating out its spatial and three-momentum dependence. This is pos- 
sible under the assumption of a homo gene ous medium and thermally equilibrated 

19071 

HQ distribution functions; one obtain a ' 

^ = -r*(iv*-7v;^) , (103) 

dr 

which now clearly exhibits detailed balance in terms of the approach to the equi- 
librium limit, N^, of the state The latter quantity is given by 

N'^"" = Vfb n^{m^-T, 7q) = d^^l J -0^f^(up; T) , (104) 

carrying the explicit dependence on the (in-medium) quarkonium mass, to^ (which, 
in turn, depends on a combination of Q-Q binding energy and in-medium HQ mass, 
cf. Eq. ([QH]) ): denotes the spin degeneracy of the 5* state and Vfb the (time- 
dependent) fireball volume. The appearance of a HQ fugacity, 7q = 7q, owes its 
origin to the (theoretically and experimentally supported) postulate that QQ pro- 
duction is restricted to hard A^-A^ collisions upon initial nuclear impact. The number 
Nqq = Nq = Nq of heavy anti-/quarks is then conserved in the subsequent fireball 
evolution (separately for charm and bottom), which is achieved by introducing jq 
at given fireball volume and temperature into the thermal densities of open and 
hidden HQ states, i.e., 

for cither charm {Q~c) or bottom {Q~b). The thermal open-charm (and -bottom) 
number, Aop, depends on whether one is evaluating it in terms of individual quark 
states, A"op = Vfb jQ2riQ{mQ,T) (appropriate for a (weakly interacting) QGP), or 
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in terms of hadronic states, Nop — Vfb7q '^ct'^(-^' A^s)- This is, in principle, a 
nontrivial issue, since both hadronic and partonic evaluations of A^op can be subject 
to corrections, see, e.g., Refs.^^nZESllOH] ^j^g hadronic phase one expects the spec- 
tral functions of D-mesons, Ac baryons, etc. to undergo significant medium effects, 
e.g., reduced masses and/or increased widths. In the partonic phase, especially close 
to Tc, it is not inconceivable that hadronic bound states are still present and thus 
an approximation with weakly interacting quasiquarks may not be an accurate one. 
Even within a quasiquark description, significant uncertainty is associated with the 
value of the HQ mass adopted in the calculation of A^^p and thus in the quantitative 
determination of 7q. The general trend is that for a given temperature, volume and 
Nqq, a smaller value for mg results in a larger value for rig and thus in a smaller 
value for 7q, which, in turn, reduces A^'' quadratically. The underlying physics is 
that of relative chemical equilibrium: for a given number of heavy anti-/quarks, the 
latter preferentially occupy the states of the lowest energy. In the simplest case, 
where a quasiquark description applies and the VE* mass is given by the expression 
([96|) . the number is essentially determined by its binding energy (larger eb im- 
plying larger N^). The gain term as written in Eq. (|102p is, strictly speaking, only 
applicable in the quasiquark approximation. If additional resonances are present in 
the medium (e.g., ZJ-meson resonances or cq diquark states), additional reaction 
channels would have to be included in a coupled rate-equation framework to ac- 
count for the competition of these resonances to harbor c quarks. In the simplified 
treatment given by Eq. (jl03p , this competition is included via the c-quark fugacity 
figuring into N^. 

A slightly differ ent v iew on regeneration and suppression processes in the QGP 
is advocated in Ref! "^^^ [ based on the strongly coupled nature of the QGP (sQGP) 
as produced at SPS and RHIC (i.e., at not too high temperatures). It is argued that 
a small charm-quark diffusion constant (cf. Sees. [2] and [3]) inhibits the separation 
of the produced c and c pair in the sQGP. In connection with the survival of J/ ip 
bound states well above Tc (as, e.g., in the right panel of Fig. [33]), this enhances 
the probability for a produced cc pair to bind into a charmonium state (relative 
to p-p collisions). In particular, this approach accounts for the possibility that the 
pairwise produced c and c quarks do not explore the entire fireball volume as usually 
assumed in equilibrium models. Such an effect has also been imple mented i n a more 
simplified manner in the thermal-rate equation approach of Refs.f^QZIIlQ] terms 
of a time dependent correlation volume. 

Let us briefly discuss the initial conditions for the quarkonium distribution 
functions. Starting point are measured quarkonium spectra in p-p collisions. In 
a heavy-ion collision, these are subject to modifications before the medium can be 
approximated with a thermal evolution. "Pre-equilibrium" effects may be distin- 
guished according to whether they occur before or after the hard QQ-production 
process takes place. The former include nuclear modifications of the parton dis- 
tribution functions generically denoted as "shadowing" , as well as pt broadening 
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(Cronin effect) attributed to a scattering of the projectile/target partons on their 
way through the target /projectile nucleus prior to the fusion reaction into Q-Q. 
In a random-walk picture, the accumulated transverse momentum is approximated 
by = agN{l), where (l) is an average nuclear path length of both gluons be- 
fore the hard scattering, and agpf parameterizes the transverse-momentum kick per 
path length in giuon-nucleon scattering. Both "pre-fusion" effects are in principle 
universal, i.e., not directly linked to the Q-Q-production process. In p-p collisions, 
a fraction of 1-2% of cc pairs (^0.1% of bb pairs) deve lop a correlation that leads to 
the formation of a charmonium (bottomonium) statJ^IHl nuclear collisions in- 
elastic collisions of the produced QQ pair with passing-by nucleons can destroy this 
correlation. This so-called nuclear absorption may be parameterized by an effective 
absorption-cross section, cr*^^!!! As is well-known, the finite (and different) forma- 
tion times of quarkonia imply that the QQ pair interacting with a nuclcon does, in 
general, not represent a fully formed quarkonium, but rather a pre-resonance state. 
A rnicroscopic description of nuclear absorption is therefore a rather challenging 
At a minimal level, finite formation times imply that the values for 
effective nuclear absorption cross sections should be expected to depend on collision 
energy {\/s), rapidity (y) and bound-state quantum numbers (since different binding 
energies imply different formation times). A careful measurement and systematic 
interpretation of quarkonium suppression in p-A collisions, where the formation of 
a thermal medium is not expected (at least at SPS and RHIC), is therefore an 
inev iteble prereq uisite for quantitative interpretations of heavy-ion data, see, e.g., 
Refs! I I I for recent work. Pre-equilibrium effects not only modify the mo- 
mentum dependence of the quarkonium distribution functions but also their spatial 
dependence. 

4.2.2. Quarkonium Phenomenology in Heavy-Ion Collisions 

As discussed in the Introduction, there is ample evidence for both chemical and 
thermal equilibration in the low-pt regime of (bulk) particle production in ultrarel- 
ativistic heavy- ion collisions. A well defined set of thermodynamic variables charac- 
terizing the temperature evolution and flow fields of the fireball greatly facilitates 
the comparison of independent calculations for quarkonium production and main- 
tains direct contact to their in-medium properties in equilibrated QCD matter. 
In this section we therefore focus on rate-equation approaches implemented into 
thermal background media. 

We recall that the experimental quarkonium yields usually include feed down 
contributions due to decays of higher resonances. E.g., for J/ip production in p-p 
collisions about 30% (10%) of the observed number arises from decays of Xc ("0') 

'Nuclear absorption typically occurs at a rather large 'J'-Af center-of-mass energy (comparable to 
the of primordial N-N collisions) and is therefore in a very different energy regime than the 
low-energy hadronic absorption cross section relevant for the later hadron-gas stage of the fireball 
evolution. 
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Fig. 36. (Color online) J/ip production at SPS i n the thermal rate-equation 

approacll2Q71190| (left 

panels) and the statistical hadronization model^S] (right panels). Upper panels display the cen- 
trality dependence of the inclusive J/ip yield (normalized to Drell-Yan production or the number 
of binary N-N collisions), lower panels the average squared J/ip transverse momentum. 



states I " I ' [ The standard assumption in heavy-ion coUisions is that primordial 
production fractions of excited states scale as in p-p collisions, but subsequent sup- 
pression (and/or regeneration) will change these ratios (due to different inelastic 
cross sections at all stages). This needs to be taken into account for realistic com- 
parisons to heavy-ion data (unless otherwise stated, it is included in the theoretical 
models discussed below). After thermal freezeout, the decay branchings are as- 
sumed to be as in vacuum (since the J/tp lifetime after freezeout (ca. 2000 fm/c) is 
about a factor of ~200 larger than the fireball lifetime, in-medium dilcpton decays 
contribute a small fraction to the spectrum observed in the detectors). 

Let us start by analyzing J/i/j pro duction in Pb-P b(^/s = 17.3 AGeV) collisions 
at SPS in the context of NA50 datJ^llinilinilHl^ cf. Fig [Ml The left panels 
display the outcome of thermal-rate equation calculationiSOHinQl^ where quasifree 
dissociation rates in the QGP (cf. left panel of Fig. [35]) and hadronic SU(4) cross 
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sections for mcson-induccd dissociation in the HG arc evolved over an expanding 
fireball model (adjusted to empirical information on hadron production and flow 
velocities). The prevalent effect is identified as suppression in the QGP, controlled by 
an effective strong coupling constant, ~ 0.25, in the quasifree rate. This value is 
adjusted to reproduce the suppression level in central collisions (where the average 
initial temperature amounts to about To ~ 210 MeV). Regeneration is a rather 
small effect, based on a p-p open-charm cross section of = 5.5 /ib distributed 
over two fireballs (the covered rapidity window amounts to Aw = 3.6, resulting in a 
rapidity density of dtTgc/dy ~ 1.53 /ib. A recent compilation^^ of charm production 
at fixed-target energies finds a total cross s ectio n of dec — 3.6-5.2 /xb; with an 
experimental rapidity width of around Ay = ^223]^ ^^j^g resulting rapidity density is 
approximately ducc/Ay — (2.2 ± 0.5) ^h). In addition, a correction for incomplete 
charm-quark thermalization has been implemente d ^ via a kinetic relaxation time 
{t^'^) reducing the equilibrium J/ip number. The lower left panel of Fig. [36l suggests 
that the centrality dependence of the average J/ip^s transverse-momentum squared, 
(Pf) . is l argely governed by the Cronin effect as extracted from experimental p-A 
dataPl^. The q uasifree charmonium dissociation rates, which increase with three- 
momentunffS^ lead to a slight suppression of (pj) for central collisions. 

T he ri ght panels in Fig. [36] are calculated within the statistical hadronization 
modePH, assuming that all primordial charmonia arc suppressed and produc- 
tion entirely occurs at the critical temperature for hadronization based on rela- 
tive chemical equilibrium of open- and hidden-charm hadrons (with Ncs fixed as 
in Eq. (jTUS])). This also implies that the charm-quark momentum distributions are 
kinetically equilibrated. With a rapidity density for the p-p charm cross section of 
dtTgc/dy = 5.7 lib the NA50 data can be reproduced reasonably well. This input 
cc number is larger by a factor of ~4 compared to the input in the left panels, 
which accounts for most of the difference to the regeneration yield in the rate- 
equation calculation for central collisions (remaining discrepancies are largely due 
to the c-quark relaxation correction which becomes more pronounced toward more 
peripheral collisions jf]. The interpretation of the J/-0's average pt is also rather 
different, in that it entirely stems from a thermal source (with moderate collective 
flow) in the vicinity of Tc (the resulting (p^) is quite consistent with the regeneration 
component in the lower left panel of Fig. [36l). 

NA50 has also measured ip' 

productioiP41225l223 

Usi ng p -A collisions, the 

extracted nuclear absorption cross section has been update (i222lto at, = (7.7 ± 
0.9) mb, which is significantly larger than for J/jp, aiLt = (4.2 ± 0.5) mb. In Pb-Pb 
collisions {^/s = 17.3 AGeV), the ratio ^'/{•JJ'^)} ^ suppressed substantially already 
in rather peri pheral collisions, cf. Fig. [g^ ^4 | 225| \Yithin the thermal rate-equation 
approacn^^, this behavior cannot be explained by inelastic reactions in the QGP 
alone, since very little (if any) QGP is formed in peripheral Pb-Pb collisions at 

™The charm ensemble at SPS is in the canonical limit, 
O.SAfop in Eq. I|105|l . and thus oc Nce- 
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Fig. 37. rcolor online) NA50 data for the tp'/{J/ii) ratio in Pb-Pb(v^ = 17.3 AG eV) 
collisionJ221Ii2Sl compared to calculations withi n the thermal rate-equation 

approaclJ2nil226| (igf^ 

panel) and the statistical hadronization model^Sl (right panel). In the left panel, the dashed 
(solid) curve is obtained without (with) the inclusion of hadronic medium effects (modeled 
by reduced D-meson masses); both calculations utilize a nuclear absorption-cross section of 

^nuc — ^^nuc — 

4.4 mb, while in the third calculation (dash-dotted line) the value has been up- 
dated to 7.9 mb (with the associated band indicating a ±0.6 mb uncertainty). In the right panel, the 
statistical model yield involves either the full hadronization volume (dash-dotted line) or excludes 
a dilute nuclear surface ("corona"), where neither suppression nor regeneration is operative. 



SPS. However, hadronic dissociation of the ip' can account for the su ppression 
pattern, but only if in-medium effects are included (cf. left panel of Fig. 
specifically a reduction of the DD threshold which accelerates ip' suppression due to 
the opening of the direct decay mode, tJj' DD (a similar effect can result from a 
broadening of the D-meson spectral functio ns as discussed in Sec. l2.7p . The updated 
(larger) ip' nuclear absorption cross scctiorP^ZI also plays a significant role in the 
quantitative description of the low-centrality data. The statistical hadronization 
model predicts a flat i/j' / {J 14') ratio, basically given by the ratio of thermal densities 
at the hadronization temperature. The shape and magnitude of the calculated ratio 
is rather consistent with the NA50 data for central and semiccntral collisions where 
hadronization from a QGP can be expected to be applicable; deviations occur for 
more peripheral centralities. Thus, the ip' /{J/^) ratio does not provide a clear 
discrimination of regeneration- and suppression-dominated scenarios at SPS. 

One of the controversies in the interpretation of the NA50 data has been 
whether they feature any "sharp" drop in their centrality dependence, e.g., around 
Et ~ 35 GeV (or iVpart ^ 120) in the upper left (right) panel of Fig. Ei Such a 
drop has been associated with a threshold behavior for QG P for mation resulting in 
an abrupt "melting" of the Xc states due to color screeninJ229| (recall that Xc feed 
down presumably makes up ~30% of the inclusive J/ip yield). The investigation of 
this question was one of the main objectives of the successor experiment of NA50, 
NA60, w here J/^p production in a medium size nuclear system (In- In) has been 
measurecP^. Fig. [38] compares the NA60 J/tJj data as a function of centrality to 
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Fig. 38. (Color online) NA60 data for the centrality dependence of J/i/) production in In-In (v^ = 
17.3 AGeM) coUision j^^^ l, compa red t o theoretical predictions based on (from top to bott om): 
(a) the threshold-melting scenariJSSSl (dashed line), (b) th e the rmal rate-equation approaclflSZ] 
(dash-dotted line) and (c) the comover suppression approacl^SQ] (dotted line) . Nuclear absorption 
effects have been divided out of the data and calculations based on measured suppression in p-A 
collisions. 



three theoretical 

prcdictioniMHZESni, 

all of which reproduce the NA50 d ata rea- 
sonably well. The predictions of the thermal rate-equation 

approadPSZUHZI 

roughly 

reproduce the onset and magnitude of the suppr ession (except for the most central 
data points) ; the threshold-melting scenaric)229| misplaces the onset of the suppres- 
sion (which in the data is b elow iVpart = 100, contrary to the Pb-Pb system) and 
the comover calculatiorPSOl overpredicts the suppression. The leveling-ofF (or even 
increasing trend) of the data for A^part > 150 is somewhat unexpected and deserves 
further study. 

The ther mal r ate-equation f rame work has been used to predict J /tp produc- 
tion at RHld2n^ Wit h upd ateJ^^Q^ for the experimental input (a smaller nuclear 
absorption cross sectiorPSS ajuj g, larger J /il) number in p-p coUisionJ^SI which 
figures into the denominator of the nuclear modification factor and leads to a rel- 
ative reduction of the regeneration yield) , an approximate agreement with current 
PHENIX data on the centrality dependence of inclusive J /if: production and pt 
spectra emerges, see left panels of Fig. |39] (the underlying charm cross section in 
p-p, (JcE = 5 70 /ib, tr anslates into dugc/dy ~ 100 /ib, consistent with PHENIX 
measurementi^SSESllj rpj^g main features of this interpretation are an about equal 
share of (suppressed) primordial and regenerated J/ip's in central Au-Au collisions 
(where the average initial temperature is about Tq = 370 MeV), as well as a sig- 



°The underlying fireball model is the same as used in the open heavy-flavor sector in connection 
with Fies. [T9l 1211 etc., with average initial temperatures of Tq ~ 340-370 MeV for scmi-/ccntral 
Au-Au collisions. 
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Fig. 39. (Color onlin e') J/ip production in 200 AGeV Au-Au collisions at RHIC in the transport 
approach of Rof.E22l (left panels, using an expanding thermal fireball with reaction rates based 
on the quasifree dissociation process) and Ref.'^'^^ (right panels, using a hydrodynamic evolution 
with reaction rates based on the gluo-dissociation process). The upper (lower) panels display the 
inclusive J/ip Raa (avera ge p ^) as a function of nucleon-participant number. The data are from 
the PHENIX coUaboratioJiSl. 



nificant reduction of the average due to secondary production, as compared to 
primordial production with an estimated Cronin effect (the latter is not yet ac- 
curately determined from p- A data). Consequently, the regeneration component is 
concentrated in the low-pt regime of the spectra. The dependence on the kinetic re- 
laxation time for c quarks is rather moderate while the inclusion of the momentum 
dependence in the quasifree dissociation ratc-'^^'^ has little effect. 

The right p anels of Fig. [39] show the results of the rate-equation approach 
of Refs.^2Q6l23T]^ where gluo-dissociation rates in the QGP are convoluted over 
a 2-|-l-dimensional hydrodynamic evolution (emp loying a charm cross section of 
dcTgc/dy = 120 ^b, in line with PHENIX datipSS^; reactions in the hadronic phase 
are neglected. It is very encouraging that the results are in good agreement with 
the ones in the left panelJ^^ which are obtained with similar physics input but 
in a different realization (e.g., fireball vs. hydro but with comparable initial tem- 
peratures (To ~ 350 MeV in central collisions) and charm cross section, etc.). The 
centrality dependenc e of the inclusive J/^p yield shows a slight step-structure in 
the upper right paneP^, induced by different dissociation temperatures for J/tJj 
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(Tdiss = 320 MeV) and Xc V'' (Tdiss = 7c = 165 McV) above which the suppression 
is a ssum ed to be practically instantaneous (similar findings have been reporte d in 
Ref. 12351). 

The (j p?) o f the regenerated component is somewhat smaller in 

Ref.l23l] 

compared to Ref.™, since in the former it is computed via continuous regeneration 
based on Eq. (|102p . while in the latter the pt spectra of the regenerated component 
are approximated with a thermal blast wave. 

Recent experimental data on J /ip production in Cu-Cu collision s indicat e that 
the nuclear modification factor tends to increase at high pt > 5 Ge\^^36|237] g^^j^ 
a trend is not present in QGP suppression calculations based on t he a uasifree disso- 
ciation rate which increases with the three-mor nentum of the However, the 
inclusion of charmonium-formation time (which reduce the dis- 
sociation rate due to time dilation in the development of the hadronic wave packet) 



and the contributions from bottom feed-dowrPSl can lead to an increase of Rjl^ 
at high p^^. 

An initially promising signature to discriminate suppression and c-c coalescence 
mechanisms is the elliptic flow of charmonia (which also provides a close connection 
to, and thus consistency check with, the collective flow of open charm). If only sup- 
pression mechanisms arc operative, the azimuthal asymmetry of the charmonium 
momentum distributions entirely develops from the path length differences caused 
by the long vs. the short axis of the almond shaped nuclear overlap zone. As for 
open-charm (recall middle panel of F ig. [24|) . this effect is rather small, generating 
a maximal V2{pt) of up to 2-3/^243| ^j^g other hand, for c-c coal escence, the 
charmonium bound state inherits up to twice the c-quark ti jl55 | 15 7j g^^. ^j^g time 
of formation, especially if the c-c binding energy is small (in that case little V2 
is carried away by an outgoing light parton). This effect is further maximized if 
the coalescence occurs late in the evolution, e.g., at the hadronization transition 
(where most of the elliptic flow is believed to have built up). However, according 
to the above discussion, the charmonium regeneration yield is mostly concentrated 
at rather low pt < 3 GeV, where the magnitude of the c-quark W2(pt/2) is not 
very large (this is a consequence of the mass ordering of V2, whose rise in pt is 
shifted to larger values for heavier particles), recall, e.g., the right panel of Fig. [TOl 
Consequently, within the transport models displayed in Fig. I39[ the net V2ipt) of 
regenerated and primordial J/^j's combined does not exceed 2-3%'^'-'^ ^"^^ "^^"^ and 
would therefore be difficult to discriminate from primordial production only. An 
interesting question, which thus far has received little attention, concerns elastic 
interactions of charmonia in the medium and whether they could contribute to 
their V2 in heavy-ion collisions. Elastic interactions should become more relevant 
as the binding energy of the charmonium increases, rendering them more compact 
objects which are less likely to break up. Interestingly, the NA60 collaboration has 
reported a rather large inclusive {pt > 0.5 GeV) elliptic flow of V2 = ( 6.8 ± 4) % 
for J/tjj's in semicentral In-In (^i = 17.3 AGeV) coUisions at the SP^^. This 
observation will be difficult to explain based on dissociation reactions alone. 
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Fig. 40. (Color online) Rapidity dependence of J/'i/i production in Au-Au [y^ = 200 AGeV) 
collisions at RHIC. Left and middle panel: results of the statistical hadronization modeJZZH for 
the Jhh rapidity density in central (left) and semicentral (middle) collisions, compared to PHENIX 
datEp2J, The central lines (shade d err or bands) correspond to (the uncertainty in) pQCD charm 
cross sections, a^^ = 256]^] /jJHS] Right panel: the ratio of ^?,4A(^part) for forward to mid- 
rapidity J/ip yields in the thermal rate-equation approaclJ2S]^ compared to PHENIX datjS]. 

Next, we address the rapidity dependence of J/ip production at RHIC, where 
PHENIX measurements in the dielectron channel at central rapidity, \y\ < 0.35, 
and in the dimuon channel, at |y| = 1.2-2.2, indicate a maximum in RAA{y) around 
y = for central and semicentral collisions. This trend can be nicely reproduced 
by the statistical hadronization model as a consequence of the underlying c-quark 
distributions (cf. left and middle panel in Fig. l40|) which are expected to be narrower 
than for bulk particle production, implying higher c-quark densities and thus larger 
charmonium yields at central rapidity. This dependence is more difficult to explain 
in the thermal rate-equation approacheJ?^ where only about 50% (or less) of the 
J/ip's originate from c-c regeneration. In addition, the charm ensemble at RHIC is 
not yet fully in the grand canonical limiiElfor which h/Io 1 in Eq. ()105p . and thus 
the sensitivity on the c-quark density (fugacity) is less pronounced. Furthermore, 
the thermal suppression of the primordial component exhibits an opposite trend, 
being slightly less suppressed at forward y due to reduced light-particle production. 
It is quite conceivable that cold-nuclcar-matter effects imprint significant rapidity 
dependencies on the prim ordial com ponen t (e.g ., stronger shadowing and/or nuclear 
absorption at forward 

In Rem, e.g. , stronger nuclear absorption of 
J ftp's at RHIC has been found as a consequence of different production mechanisms 
which probe different kinematics in the nuclear parton distribution functions. 

Recent calculations of charmonium production within microsc opic trari sport 
models for the bulk-medium evolution can be found, e.g., in Refs .2232111. The 
results are generally quite reminiscent of the rate-equation calculations discussed 
above. In particular, the description of RHIC data requires the inclusion of regen- 
eration interactions. This is also true for the so-called comover approach, which has 




°In part, this is due to a finite correlation volume introduced for cc quarks in the approach of 
RefsMMl, 
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been extended to include charmonium formation reactions in 

Re« 

The interplay of suppression and regeneration should lead to interesting conse- 
quences for the excitation function of charmonium production. The approximate 
degeneracy of J/i/' suppression by about a factor of ~3 in central A- A collision s at 
both SPS and RHIC has been anticipated in the two-component model of Ref!^^^ 
with a rather flat behavior for ^ = 17-200 AGeV. This degeneracy is expected to 
be lifted at higher (LHC, = 5500 AGeV) and lower (FAIR, ^ = 8AGeV) 
collision energies. At LHC, the statistical hadronization model predicts the in- 
clusive J/ip Raa in central Pb-Pb collisions to recover the level in p-p colli- 
sions, i.e., RAA{Npait = 350) 1, based on a p-p open-charm cross section of 
darr Jdy = 640 /ib. At F AIR e nergies, on the other hand, statistical production is 
smalP^, while transporlPSDl and rate-equation approaches predict about a factor 
of two suppression, mostly dominated by nuclear absorption. The effective nuclear 
absorption cross section will thus be an essential quantity to be determined in p-A 
reactions at FAIR. 

Finally, let us turn to bottomonium production, which adds several new aspects 
compared to charmonium production: (i) the binding energies of bottomonium 
states are larger by about a factor of ~2 which opens a wider window to study 
their dependence on color screening (due to larger dissociation temperatures) and 
makes them more robust in the hadronic phase; (ii) at given collision energy, the 
number of bb pairs is substanti ally srn aller than the number of cc ones (e.g., by 
about a factor of ^^200 at (iii) bottom-quarks are less susceptible to 

changes in their momentum distributions due to their factor '--^S larger mass (as 
discussed in Sees. [2] and [3] of this review). The latter 2 points suggest that regen- 
eration processes play less of a rolc|p|. Early analyses of T production in heavy-ion 
collisions have focused on the pt-dependence of suppression scenarios where in- 
stantaneous dissociation above a cr itical energy density has b een combined with 
formation-time effects, both at LHd -^4 | 251 | 252r ^^d RHld2521. The opposite limit 
of secondary production alone has been evaluated in the statistical hadronization 
modc|221l. The thermal rate-equation approach, Eq. (|103p . has been applied to T 
production in Ref.l^, in analogy to the charmonium sector as displayed in the left 
panels of Figs. [36land[39l The time evolution of the T{1S) yield in central Au-Au 
collisions at RHIC has been calculated for the following two scenarios (as shown in 
Fig. l4ip : in the flrst one (left panel), reduced in-medium binding energies (acco rding 
to solutions of a Schrodinger equation with a color-screened Cornell potcntiaff^l^ 
are combined with quasifree dissociation (and formation) reactions; in the second 
one, the gluo-dissociation process is applied to the T(IS') with vacuum binding en- 
ergy (assuming rub = 5.28 GeV in connection with eb ~ 1-1 GeV). One finds that 



PRowever, care has to be taken in deducing that this renders T regeneration irrelevant at RHIC, 
since (a) the bottom ensemble is in the canonical limit, and (b) the regeneration yield needs to 
be compared to the primordial yield: in p-p collisions the ratio T /{bb) ~ 0.1% is about a factor of 
10 smaller than in the charm sector where J/i/'/(cc) ~ 1%. 
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Fig. 41. (Color o nline ) Time dependence of T(1S) production in central Au-Au(y^ = 200 AGeV) 
collisions at RHidSS] using the quasifree dissociation-cross section with in-medium reduced binding 
energies (left panel) and the gluo-dissociation cross section with the free binding energy (assumed 
to be = 1.1 GeV; left panel). The time evolution of the background medium is based on an 
expanding fireball as in the left panels of Fig. 1391 



color-screening enables a ^40% suppression of the T(IS') within the first 1-2 fm/c 
(where the average medium temperature is above 200 MeV), with insignificant con- 
tributions from regeneration since is too small. On the other hand, with its 
vacuum binding energy, the T(15') is basically unaffected at RHIC. This suggests 
a very promising sensitivity of T{1S) production to color-screening at RHIC (of 
course, the observed T(IS') yield contains feed down contributions, amounting to 
ca. 50% in p-p collisions; this underlines again the importance of measuring the 
excited T' and Xb states) . Note that A''^'' is quite different in the two scenarios even 
though rriT = 2mf, — is assumed to be constant in both cases. The reason is the 
difference in relative chemical equilibrium: a larger binding energy implies a larger 
rrib, which leads to a larger via Eq. (jlOSp . i.e., it is thermodynamically more 
favorable to allocate b and b quarks in an T (relative to a smaller with smaller 
mt and thus smaller 7^). 

5. Conclusions 

Heavy-quark physics is an increasingly useful and adopted tool in the theoretical 
analysis of hot and dense QCD matter and its study in ultrarelativistic collisions 
of heavy nuclei. Key reasons for this development lie in a combination of exciting 
new data becoming available not only for quarkonia but also for open heavy-flavor 
observables, together with attractive features of charm and bottom quarks from 
a theoretical point of view. These features are, of course, rooted in the large scale 
introduced by the heavy-quark (HQ) mass, which enables the use of expansion tech- 
niques, most notably HQ effective theories and Brownian motion for HQ diffusion 
in a QGP fluid. Furthermore, the production of heavy quarks in nuclear reactions 
is presumably restricted to primary N-N collisions which renders HQ spectra a cal- 
ibrated probe of the medium over the entire range of transverse momentum. The 
latter provides a unique opportunity for a comprehensive investigation of the QCD 
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medium at all scales, ranging from diffusion physics in the low-pt limit to a "stan- 
dard" hard probe at pt S> toq. Heavy quarks thus connect transport coefficients 
in the QGP and observables in ultrarelativistic heavy-ion collisions in the arguably 
most direct way. Finally, relations between the open and hidden heavy-flavor sectors 
promise valuable mutual constraints, both theoretically and phenomenologically. In 
this review we have largely focused on aspects of soft physics for HQ propagation 
and binding in the QGP. 

Interactions of slowly moving heavy quarks in a QGP are dominated by clas- 
tic scattering off thermal partons. A perturbative expansion of these interactions, 
specifically for the HQ diffusion coefficient, shows poor convergence for coupling 
constants believed to be relevant for a QGP as formed in heavy-ion collisions. Sev- 
eral options of amending the perturbative treatment have been suggested, e.g., a 
reduced screening mass or running coupling constant at low momentum transfer. 
While increasing the interaction strength, they inevitably face the problem of little 
control over higher order "corrections" . Nonperturbative approaches have been put 
forward which can, in principle, overcome this problem by a (partial) resummation 
of large contributions. E.g., a potential-based T-matrix approach characterized by 
a scattering equation becomes particularly promising if the input interaction can 
be specified in a model-independent way, i.e., from thermal lattice QCD. Currently, 
open questions remain as to the validity of the potential approach at finite tem- 
perature, and a suitable definition of the potential from the HQ free energy. Here, 
a close connection between the open and hidden heavy-flavor sectors emerges via 
the same low-energy interaction operative for HQ diffusion and quarkonium bound 
states. Qualitatively, one finds that, if ground-state quarkonia survive until tem- 
peratures well above T^, potential scattering of heavy quarks in the QGP builds up 
resonance-like correlations which are instrumental in obtaining a small HQ diffusion 
coefficient close to Tc- As an alternative nonperturbative approach, HQ diffusion has 
been estimated in the strong-coupling limit of conformal field theory (CFT) using 
a conjectured correspondence to string theory in Anti-de-Sitter (AdS) space. With 
CFT parameters adapted to resemble QCD, the resulting HQ diffusion constant 
is comparable to the T-matrix approach close to Tc, but is approximately con- 
stant with increasing temperature while the T-matrix interaction and ultimately 
approaches pQCD estimates. 

Quantitative applications of the Brownian-motion framework to HQ observables 
at RHIC critically hinge on a reliable description of the background medium evo- 
lution. The latter specifies the ambient conditions in the (approximately) thermal 
bath including its temperature and collectivity, whose magnitudes directly impact 
the nuclear modification and elliptic flow of HQ spectra. A survey of available 
calculations indicates that the translation of a given HQ diffusion coefficient into 
suppression and elliptic flow of HQ spectra is currently at the ~50% accuracy level. 
This needs to be further scrutinized and improved. In line with the theoretical 
expectations for HQ diffusion in a QGP at T=l-2Tc, the current data call for sig- 
nificantly stronger interactions than provided by LO pQCD. A simultaneous and 
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consistent evaluation of spectra (RAAiPt)) and elliptic flow {v2{pt)) is pivotal to this 
conclusion. A non-negligible role is played by the hadronization process: heavy-light 
quark coalescence processes seem to improve the experimentally observed correla- 
tion of a rather large V2 and a moderately suppressed Raa in the electron spectra at 
moderate pfp<3 GeV. More robust conclusions will require a better knowledge of the 
Cronin effect figuring into the initial conditions for the HQ spectra. A theoretically 
appealing aspect of quark coalescence is its close relation to resonance correlations 
in the QGP, which could be at the origin of the nonperturbative interaction strength 
in HQ diffusion. From the experimental side, important discrimination power will 
come with an explicit measurement of Z?-mesons, to explicitly separate the bottom 
contribution (present in the electron spectra). Theoretical investigations should take 
advantage of the opportunity to predict angular correlation measurements which 
will become feasible soon. We believe that charm data at SPS energy would consti- 
tute a valuable complement to RHIC data, which could help in deciding how much 
of the bulk flow (and suppression) can be imparted on charm quarks in a medium 
at smaller temperatures. 

In the heavy-quarkonium sector we have started with a brief synopsis of cur- 
rent applications of potential models in medium. At this stage, the comparison of 
calculated spectral functions to Euclidean correlation functions computed in lattice 
QCD suggests that scenarios with either strong binding and rather large in-medium 
HQ mass, or weak binding and smaller HQ mass, are both viable. On the contrary, 
inelastic reaction rates are rather sensitive to the binding energy and thus to the 
strength of color-Debye screening, especially for bottomonia. This translates into 
a promising discrimination power of bottomonium suppression measurements at 
RHIC and LHC. The phenomenology of charmonia is presumably more involved; 
e.g., kinetic rate-equation calculations for J/ip production in central Au-Au colli- 
sions at RHIC indicate that the number of surviving primordial J/i/j's is comparable 
to the number of secondary produced ones via c-c coalescence in the QGP (and/or at 
hadronization) . Such an interpretation is consistent with J/ijj px spectra where the 
coalescence yield is concentrated at low pT, thus reducing the average p^ compared 
to primordial production. Kinetic approaches furthermore suggest that regeneration 
is subleading at the SPS, and that the observed suppression occurs at energy den- 
sities above the critical one. Deeper insights will follow when advanced theoretical 
models meet future precision data, including rapidity dependencies, elliptic flow, 
excited charmonia and much needed constraints from d-Au collisions at RHIC (to 
pin down cold-nuclear matter effects). An extended excitation function via a RHIC 
energy scan, LHC and FAIR will further disentangle suppression and regeneration 
effects. We emphasize again that the presence of regeneration mechanisms would 
imply valuable connections to the open heavy-flavor sector, as coalescing heavy 
quarks necessarily imprint their kinematics on HQ bound states. 

In summary, we believe that in-medium HQ physics will continue as a chal- 
lenging but rewarding forefront research field for many years to come, with ample 
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opportunities for surprises, insights and progress. 
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